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Abstract 

We lay down the foundations of the theory of Poisson vertex algebras aimed at its 
applications to integrability of Hamiltonian partial differential equations. Such an equation 
is called integrable if it can be included in an infinite hierarchy of compatible Hamiltonian 
equations, which admit an infinite sequence of linearly independent integrals of motion in 
involution. The construction of a hierarchy and its integrals of motion is achieved by making 
use of the so called Lenard scheme. We find simple conditions which guarantee that the 
scheme produces an infinite sequence of closed 1-forms ujj , j £ Z + , of the variational complex 
ft. If these forms are exact, i.e. uij are variational derivatives of some local functionals 
Jhj, then the latter are integrals of motion in involution of the hierarchy formed by the 
corresponding Hamiltonian vector fields. We show that the complex f2 is exact, provided 
that the algebra of functions V is "normal"; in particular, for arbitrary V, any closed form 
in Q becomes exact if we add to V a finite number of antiderivatives. We demonstrate 
on the examples of the KdV, HD and CNW hierarchies how the Lenard scheme works. 
We also discover a new integrable hierarchy, which we call the CNW hierarchy of HD type. 
Developing the ideas of Dorfman, we extend the Lenard scheme to arbitrary Dirac structures, 
and demonstrate its applicability on the examples of the NLS, pKdV and KN hierarchies. 

Keywords and phrases: evolution equation, evolutionary vector field, local functional, integral 
of motion, integrable hierarchy, normal algebra of differential functions, Lie conformal algebra, Poisson 
vertex algebra, compatible A-brackets, Lenard scheme, Beltrami A-bracket, variational derivative, Frcchct 
derivative, variational complex, Dirac structure, compatible Dirac structures. 
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Introduction. 



An evolution equation is a system of partial differential equations of the form 
(0.1) — - =Pi(u,u',u",...), ie I = {!,...,£} (I may be infinite) , 

where ui = Ui(t,x) are functions on a 1-dimensional manifold M, depending on time t and 
x G M, and Pi are differentiable functions in u = (iti)iej and a finite number of its derivatives 

/ _ ldui\ II _ (GPui\ 

One usually views UijU'^u'- , . . . as generators of the algebra of polynomials 
(0.2) R = C[u\ n) \i el, neZ + ] , 

equipped with the derivation d, defined by d(u^) = u[ n+1 \ n G Z+. An algebra of differential 
functions V is an extension of the algebra R, endowed with commuting derivations - ^ , i G 



du 

M 

du 

for all but finitely many i G / and n G Z + . Then d extends to a derivation of V by 



/, n G Z+, extending the usual partial derivatives in ii, and such that, given / G V, ^ = 



(0.3) 8= y n ; (n+1) ^-. 

An element of V/dV is called a /oca/ functional, and the image of / G V in V/<9V is de- 
noted by J* /, the reason being that in the variational calculus local functionals are of the form 
J M f(u, u', . . .) dx, and taking V/dV provides the universal space for which integration by parts 
holds. A local functional //is called an integral of motion of (0.1), or is said to be conserved 
by the evolution equation (0.1), if its evolution in time along (0.1) is constant, namely 

<°- 4 > /I - °- 

The element /, defined modulo dV, is then called a conserved density. By the chain rule, 
equation (0.4) is equivalent to 

(0.5) /*>(/)= 0, 

where, for P = (Pj)j g / G V , Xp is the following derivation of V: 

(0.6) A' 

Note that [9, Xp] = 0. A derivation of the algebra V with this property is called an evolutionary 
vector field. It is easy to see that all of these are of the form (0.6) for some P G V . 

Here and further V denotes the space of all I x 1 column vectors with entries in V. Also, 
V®^ will be used to denote the subspace of i x 1 column vectors with only finitely many non-zero 
entries. Though, in this paper we do not consider any example with infinite I, we still use this 
distinction as a book-keeping device. 
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We have the usual pairing 

v ®i x y£ Vj defined by FP = £\ gJ FjPj. Note that integrating 

by parts transforms (0.5) to 



ou 



"here & = (*) € V « and 



(0.7) V _ ^ ,_ ar 



is the variational derivative of / by Uj. 

Given a sequence of commuting evolutionary vector fields, i.e. a sequence P n E V^, n = 
0, 1, 2, . . . , such that the corresponding evolutionary vector fields commute (compatibility con- 
dition): 

(0.8) [X pm ,X P n] = , for all m,n G Z + , 

one considers a hierarchy of evolution equations 

(0.9) ^ = P n (u, u>y,...), neZ + , 

dt n 

where jj^- = (gr*) ie j £ V* and U{ = Ui(ti,t2, ■ ■ ■)■ If the evolutionary vector fields Xpn, n £ Z + , 
span an infinite-dimensional vector space, then each of the equations of the hierarchy (0.9) (and 
the whole hierarchy) is called integrable. For example, the hierarchy of linear evolution equations 

(0.10) ^ = ut ] , iG/,nGZ + , 

ut n 

is integrable. Note that the compatibility condition (0.8) can be equivalently expressed by the 
following identities: 

dP m dP n 

-n~ = ' rn,neZ + . 
utn aT m 

Note also that the problem of classification of integrable evolution equations is equivalent to the 
problem of classification of maximal infinite-dimensional abelian subalgebras of the Lie algebra 
of evolutionary vector fields. 

Remark 0.1. The above setup can be generalized to the case of any finite-dimensional manifold 
M, replacing d by partial derivatives di = and V/dV by V /'Y^ i {diV). However, we shall be 
concerned only with the 1-dimensional case. 

A special case of an evolution equation is a system of equations of the form 

where H (d) = (Hij (u, u', . . . ; <9)) i - eJ is an I x I matrix, whose entries are finite order differential 
operators in d with coefficients in V, and J h € V/dV is a local functional (note that, since, in 
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view of (0.3), -J- o d = Ofor all i £ I, the RHS is well defined). In this case one considers the 
"local" bracket "[FT] G I): 

(0.12) {ui(x),Uj(y)} = H ji (u(y),u'(y),...;d y )S(x-y), x,y G M , 

where 5(x — y) is the <5-function: j M f(x)5(x — y) dx = f{y). This bracket extends, by the 
Leibniz rule and bilinearity, to arbitrary /, g G V: 

(0.13) {f(x),g(y)}= £ u,(y)} . 

Using bilinearity and integration by parts, equation (0.13) suggests a way to define a bracket 
on the space of local functionals V/dV. Indeed we have 



{I M f(.x)dx, f M g(y)dy} = V f [ ^{ X 5 ° V {ui(x),Uj(y)}dxdy 

£ / 5 ^H j My),u\y),... i d y) 8 Mdy, 



where, for the last identity, we performed integration by parts. Hence, given two local func- 
tionals j f,Jg G V/dV, we define their bracket associated to the operator H(d), as the local 
functional 

Likewise, given a local functional J f G V/dV and a function 5 G V, we can define their bracket, 
which will now be an element of V, by the formula 

(0.15) {ff,g} = «',... 

i,jei du j 1 
nez+ 

Of course, by integration by parts, we get that the brackets (0.14) and (0.15) are compatible in 
the sense that 

I {If, 9} = UfJg}- 

We can then rewrite the evolution equation (0.11), using the above notation, in the Hamiltonian 
form: 

(0.16) ^ = {Jh,u}. 

Here and further, for u = (uj)ier £ V> {JX - "} stands for Uj})^ ^ G V £ . 

The bracket (0.13) is called a Poisson bracket if the bracket (0.14) on V/dV satisfies the 
Lie algebra axioms (this is the basic integrability conditions of the evolution equation (0.11)). 
The skew-commutativity of the bracket (0.14) simply means that the differential operator H{d) 
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is skew-adjoint. The Jacobi identity is more involved, but can be easily understood using the 
language of A-brackets, which, we believe, greatly simplifies the theory. The idea is to apply 
the Fourier transform F(x,y) i— > J M dxe x ( x ~ y ^ F(x, y) to both sides of (0.13). Denoting 

{f(y)xg(y)} = f M e x ^{f(x), g ( y )}dx, 

a straightforward calculation gives the following important formula [DK]: 

(0.17) {f x g}= £ J^id + Xnuta+wUi-d-Xr-^, 

i,jei du j du i 

where {v,i\Uj} = Hji(u, u', . . . ; A), and {uig +x Uj}^ means that d is moved to the right. 

In terms of the A-bracket (0.17), the bracket (0.14) on the space of local functionals V/dV 
becomes 



(0-18) {Jfjg} = J{hg} 



A=0 ' 



Likewise the bracket (0.15) between a local functional Jf € V/dV and a function g € V can be 
expressed as 

{If, 9} = {f\g}\ x=0 , 
so that the evolution equation (0.11) becomes 

(0.19) ^ = {hxu} \ x=Q . 

It is easy to see that the A-bracket (0.17) satisfies the following sesquilinearity properties 
(0-20) {df x g} = -X{f x g} , {f x dg} = (8 + X){f x g} , 

and the following left and right Leibniz rules: 

(0.21) {f x gh} = {f x g}h + g{f x h} , {fg x h} = {f x+d h}^g + {g x+9 h}^f . 

Furthermore, the skew-commutativity of the bracket (0.14) is equivalent to 

(0-22) {f x g} = -{g-g-xf} , 

where now d is moved to the left, and the Jacobi identity is equivalent to 

(0-23) {fx{9»h}} - {gMxh}} = {{fxg}x+„h} . 

A commutative associative unital differential algebra V, endowed with a A-bracket V® V — > 
C[A] <8> V, denoted by a (8) b 1— ► {a x b}, is called a Poisson vertex algebra (PVA) if it satisfies all 
the identities (0.20)-(0.23). If the A-bracket (0.17) defines a PVA structure on V, we say that 
H(d) = \Hij(u, u', . . . ; 5)) . . j in (0.12) is a Hamiltonian operator, and that equation (0.11) 
(or, equivalently, (0.19)) is a system of Hamiltonian equations associated to the Hamiltonian 
operator H{&). 
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It follows from (0.17) that 



{fh, •} = {h X -}\ x=0 = X H(d)5h/5u for fh G V/dV . 

It is easy to check that, provided that Jacobi identity (0.23) holds, the map f f t— > X H (g)Sf/8u 
defines a homorphism of the Lie algebra V/dV of local functionals (with the bracket (0.18)) to 
the Lie algebra of evolutionary vector fields. Recall that, by definition, a local functional f f is 
an integral of motion of the Hamiltonian equation (0.16) if 

(0-24) jft ={Ih ' If} = - 

By the above observation this implies that the corresponding evolutionary vector fields X H ^sh/s u 
and X H ^sf/Su commute: 

(0-25) [Xu(d)5h/Sui X H(d)Sf/Su] = • 

In fact, we will see in Section 1.6, that in many cases this is a necessary and sufficient condition 
for the commutativity relation (0.24). 

The basic problem in the theory of Hamiltonian equations is to establish "integrability" . 
A system of Hamiltonian equations du/dt = {Jh,u} = H(d)5h/5u, is said to be integrable 
if, first, fh lies in an infinite-dimensional abelian subalgebra of the Lie algebra (V/dV, { , }), 
in other words, if there exists an infinite sequence of linearly independent local functionals 
fho = fh, fh±, fh2, ■ ■ ■ G V/dV, commuting with respect to the Lie bracket (0.18) on V/dV, 
and, second, if the evolutionary vector fields X H ^sh n span an infinite-dimensional vector 

space. By the above observations, if such a sequence fh n exists, the corresponding Hamiltonian 
equations 

^ = {Jhn,u} = H(d) 6 -^ , n G Z + , 

form a hierarchy of compatible Hamiltonian equations. Thus, if the center of the Lie algebra 
V/dV with bracket (0.18) is finite-dimensional, then we get a hierarchy of integrable evolution 
equations. (Note that the hierarchy (0.10) of all linear evolution equations is not Hamiltonian, 
but its part with n odd is.) 

The problem of classification of integrable Hamiltonian equations consists of two parts. The 
first one is to classify all A-brackets on V, making it a PVA. The second one is to classify all 
maximal infinite-dimensional abelian subalgebras in the Lie algebra V/dV with bracket (0.18). 

What apparently complicates the problem is the possibility of having the same evolution 
equation (0.19) in two different Hamiltonian forms: 

(°- 26 ) ft = U h i\ u } \\=o = U h 0\ U }l\\=0> 

where {• \ i = 0, 1, are two linearly independent A-brackets on V and fhi, i = 0, 1, are some 
local functionals. However, this disadvantage happens to be the main source of integrability in 
the Hamiltonian approach, as can be seen by utilizing the so called Lenard scheme [M], [O], 
[D]. Under some mild conditions, this scheme works, provided that the two A-brackets form a 
bi-Hamiltonian pair, meaning that any their linear combination makes V a PVA. Namely, the 
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scheme generates a sequence of local functionals Jh n , n G Z+, extending the given first two 
terms, such that for each n G Z + we have 



(0.27) 



{Jh n+lx u 



} 



0lA=0 



{Jh nx u} 



1 1 A=0 ' 



In this case all Jh n G V/dV, n G Z + , pairwise commute with respect to both brackets 
{• , i = 0, 1, on V/9V, and hence they are integrals of motion in involution for the evo- 
lution equation (0.26). 

In Section 2 we provide some conditions, which garantee that the Lenard scheme works. As 
applications, we discuss in detail the example of the KdV and the dispersionless KdV hierarchies, 
based on the Gardner-Faddeev-Zakharov and the Virasoro-Magri PVA. In particular, we prove 
that the polynomials in u form a maximal abelian subspace of integrals of motion for the 
dispersionless KdV, and derive from this that the infinite-dimensional space of integrals of 
motion for the KdV hierarchy, obtained via the Lenard scheme, is maximal abelian as well. In 
fact, due to the presence of a parameter in the Virasoro-Magri A-bracket (the central charge), 
we have a triple of compatible A-brackets. This allows one to choose another bi-Hamiltonian 
pair, which leads to the so called HD integrable hierarchy, which we discuss as well. We also 
discuss the example of the coupled non-linear wave (CNW) system of Ito, based on the PVA 
corresponding to the Virasoro Lie algebra acting on functions on the circle. Here again, due to 
the presence of a parameter in the A-brackets, we can make another choice of a bi-Hamiltonian 
pair, leading us to a new hierarchy, which we call the CNW hierarchy of HD type. 

Another important class of Hamiltonian equations is provided by symplectic operators. Let 
S(d) = (Sij(u,u' , . . . ; 3)) . - eJ be a matrix differential operator with finitely many non-zero 
entries. The operator S(d) is called symplectic if it is skew-adjoint and it satisfies the following 
analogue of the Jacobi identity (i,j, k G I): 



Here {■ \-}b is called the Beltrami A-bracket, which we define by {ui\Uj} = 5ij and extended 
to V ® V -> C[A] tg> V using (0.20) and (0.21). 

In the symplectic case there is again a Lie algebra bracket, similar to (0.18), but it is defined 
only on the following subspace of V/dV of Hamiltonian functionals: 



(0.28) 



{u i xS kj (^)} B -{u jl j,S k i(X)} B + {S j i(X)x+nU k } B = 0. 



It is given by the following formula (cf. (0.15)): 



(0.29) 




where — = S(d)P . 



Su 



Consider an evolution equation of the form 



(0.30) 



dt 



where P G V e is such that S(d)P = ^ for some Jh G J~s- A local functional jf G J~s is 
called an integral of motion of the evolution equation (0.30) if {Jh,J f}s = 0, and equation 
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(0.30) is called integrable if it admits infinitely many linearly independent commuting integrals 
of motion Jh n ,n G Z+, commuting with jh = fho, and the corresponding evolutionary vector 
fields Xpn commute and span an infinite-dimensional vector space. Thus, classification of 
integrable Hamiltonian equations associated to the symplectic operator S(d) reduces to the 
classification of infinite-dimensional maximal abelian subspaces of the Lie algebra J~s- Based 
on the above definitions, the theory proceeds in the same way as in the Hamiltonian case, 
including the Lenard scheme. Following Dorfman [D], we establis in Section 4 integrability of 
the potential KdV equation and the Krichever-Novikov equation, using this scheme. 

In fact, there is a more general setup, in terms of Dirac structures, introduced by Dorf- 
man [D], which embraces both the Hamiltonian and the symplectic setup. A Dirac struc- 
ture involves two matrix differential operators H(d) = (fly (it, u', . . . ; d)^. and S(d) = 

(Sij(u, u', . . . ; 9)) i . g p and the corresponding Hamiltonian equations are of the form ^ = P , 

where P G V e is such that 

(0.31) S(d)P = H(d)^, 

for some JhG. Fh Si the corresponding Lie algebra of Hamiltonian junctionals, defined by 

(0.32) T H ,s = {// € V/8V 

with the Lie algebra bracket 
(0.33) {Jfj9}s = where fl(5)^ = 5(a)P. 

Under a suitable integrability condition, expressed in terms of the so-called Courant-Dorfman 
product, (0.33) is well defined and it satisfies the Lie algebra axioms. 

We develop further the theory of the Lenard scheme for Dirac structures, initiated by 
Dorfman [D], and, using this, complete the proof of integrability of the non-linear Schrodinger 
system, sketched in [D]. 

Applying the Lenard scheme at the level of generality of an arbitrary algebra of differen- 
tial functions V requires understanding of the exactness of the variational complex f2(V) = 
(Bkez+ nk ( v ) at A; = and 1. We prove in Section 3 that, adding finitely many antiderivatives 
to V, one can make any closed fc-cocycle exact for any k > 1, and, after adding a constant, for 
k = (which seems to be a new result). 

The contents of the paper is as follows. In Sections 1.1 and 1.2 we introduce the notion 
of an algebra V of differential functions and its normality property. The main results here are 
Propositions 1.5 and 1.9. They provide for a normal V (in fact any V after adding a finite 
number of antiderivatives) algorithms for computing, for a given / 6 V such that g = 0, 
an element g 6 V such that / = <9<7+const., and for a given F S V®^ which is closed, i.e. 
such that Dp{d) = D F (d), where Dp{d) is the Frechet derivative (defined in Section 1.2), an 
element / G V such that x£ = F. These results are important for the proof of integrability of 
Hamiltonian equations, discussed in Sections 2 and 4. 

In Sections 1.3 and 1.7 we give several equivalent definitions of a Poisson vertex algebra 
(PVA) and we explain how to introduce a structure of a PVA in an algebra V of differential 



H(d)f G S(d)V £ \ , 
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functions or its quotient (Theorems 1.15 and 1.21, Propositions 1.16 and 1.37). In particular, 
we show that an equivalent notion is that of a Hamiltonian operator. In Section 1.5 we show 
how to construct Hamiltonian equations and their integrals of motion in the PVA language. In 
Section 1.6 we explain how to compute the center of a PVA. In Section 1.4 we introduce the 
Beltrami A-bracket (which, unlike the Poisson A-bracket, is commutative, rather than skew- 
commutative), and interpret the basic operators of the variational calculus, like the variational 
derivative and the Frechet derivative, in terms of this A-bracket. Other applications of the 
Beltrami A-bracket appear in Section 3. 

In Section 2.1 we introduce the notion of compatible PVA structures on an algebra V of 
differential functions and develop in this language the well-known Lenard scheme of integrability 
of Hamiltonian evolution equations [M], [O], [D]. The new results here are Propositions 2.9, 
2.10 and 2.13, and Corollary 2.12, which provide sufficient conditions for the Lenard scheme to 
work. 

In Section 2.2, using the Lenard scheme, we discuss in detail the integrability of the KdV 
hierarchy, which includes the classical KdV equation 

— = 3un' + cv!" , ceC. 
dt 

The new result here is Theorem 2.15 on the maximality of the sequence of integrals of motion of 
the KdV hierarchy. We discuss the integrability of the HD hierarchy, which includes the Harry 
Dim equation 

^ = 03(„-V2) +a d{u- 1 ' 2 ) , aeC, 

in Section 2.3, and that of the CNW hierarchy, which includes the CNW system of Ito (c £ C): 

i % = cv"' + 3uu' + vv' 
X I = d(uv) 

in Section 2.4. In Section 2.5 we prove integrability of the CNW hierarchy of HD type, which 
seems to be new. The simplest non-trivial equation of this hierarchy is (c G C): 

X § = 

In Sections 3.1 and 3.2 we introduce, following Gelfand and Dorfman [GD2], the variational 
complex il(V) as the reduction of the de Rham g-complex f2(V) by the action of d, i.e. O(V) = 
£l(V)/dQ(V). Our main new result on exactness of the complexes Q(V) and O(V) provided 
that V is normal, is Theorem 3.2, proved in Section 3.3. Another version of the exactness 
theorem, which is Theorem 3.5, is well known (see e.g. [D], [Di]). However it is not always 
applicable, whereas Theorem 3.2 works for any V, provided that we add to V a finite number 
of antiderivatives. 

In Sections 3.4 and 3.5, following [DSK], we link the variational complex Q(V) to the Lie 
conformal algebra cohomology, developed in [BKV], via the Beltrami A-bracket, which leads 
to an explicit construction of f2(V). As an application, we give a classification of symplectic 
differential operators and find simple formulas for the well-known Sokolov and Dorfman sym- 
plectic operators [S], [D], [W]. In Section 3.6 we explore an alternative way of understanding 
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the variational complex, via the projection operators Vk of Cl onto a subspace complementary 

to on k . 

In Sections 4.1-4.3 we define the Dirac structures and the corresponding evolution equations. 
In Sections 4.4 and 4.5 we introduce the notion of compatible Dirac structures and discuss the 
corresponding Lenard scheme of integrability. The exposition of these sections follows closely 
the book of Dorfman [D], except for Proposition 4.16 and Corollary 4.17, guaranteeing that the 
Lenard scheme works, which seem to be new results. 

In Section 4.6 we prove integrability of the non-linear Schrodinger (NLS) system 

f 3f = v" + 2v(u 2 + v 2 ) 

\ f = -u" - 2u(u 2 + v 2 ) ' 

using a compatible pair of Dirac structures and the results of Section 4.5. 

In Section 4.7 we prove that the graph of a Hamiltonian differential operator is a Dirac 
structure (see Theorem 4.21), and a similar result on the bi-Hamltonian structure vs. a pair of 
compatible Dirac structures (see Proposition 4.24), which allows us to derive the key property 
of a bi-Hamiltonian structure, stated in Theorem 2.7 in Section 2.1. 

Likewise, in Section 4.8 we relate the symplectic differential operators to Dirac structures 
and prove Dorfman's criteria of compatibility of a pair of symplectic operators [D]. In Sections 
4.9 and 4.10 we derive the integrability of the potential KdV equation (c € C) 

f t = 3(u') 2 + cu'» , 

and of the Krichiver-Nivikov equation 

du „, 3 {u"f 
~dt ~ U 2 u' ' 

The latter uses the compatibility of the Sokolov and Dorfman symplectic operators. 

In this paper we are considering only the translation invariant case, when functions / E V 
do not depend explicitly on x. The general case amounts to adding J| in the definition (0.3) 
of d. 

Many important integrable hierarchies, like the KP, the Toda lattice, and the Drinfeld- 
Sokolov, are not treated in this paper. We are planning to fill in these gaps in future publications. 



1 Poisson vertex algebras and Hamiltonian operators. 

1.1 Algebras of differential functions. By a differential algebra we shall mean a unital 
commutative associative algebra over C with a derivation d. Recall that dl = 0. One of the 
most important examples is the algebra of differential polynomials in i variables 

R = £[uf ] \ i e {1, ...,£} = I, n e Z+] , 

where d is the derivation of the algebra R, defined by d(iq ) = u\ n+1 \ 

Definition 1.1. An algebra of differential functions V in a set of variables {ui}i^i is an ex- 
tension of the algebra of polynomials R = C[u[ \i € I,n G Z+], endowed with linear maps 
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— j^r : V — ► V, for all i G / and n G Z + , which are commuting derivations of the product in 
V, extending the usual partial derivatives in R, and such that, given / £ V, — f^y/ = for all 

but finitely many i G I and n G Z+. Unless otherwise specified, we shall assume that V is an 
integral domain. 

Typical examples of algebras of differential functions that we will consider are: the algebra 
of differential polynomials itself, R = C[u[ n ^ \i G I,n G Z+], any localization of it by some 
element f & R, or, more generally, by some multiplicative subset S C R, such as the whole field 
of fractions Q = C{u[^ \ i G I, n G Z+), or any algebraic extension of the algebra R or of the 
field Q obtained by adding a solution of certain polynomial equation. An example of the latter 
type, which we will consider in Section 2.3, is V = ^-[y/u^ 1 ,u" , . . obtained by starting 
from the algebra of differential polynomials R = C[u^ \ n £ Z+], adding the square root of the 
element u, and localizing by ^Ju. 

On any algebra of differential functions V we have a well defined derivation d : V — > V, 
extending the usual derivation of R (defined above), given by 



(1-1) d = Y, 



(n+l)_9_ 



Moreover we have the usual commutation rule 

d 



(1.2) 



9 8 



where the RHS is considered to be zero if n = 0. These commutation relations imply the 
following lemma, which will be used throughout the paper. 

Lemma 1.2. Let D^z) = E„ e z+ * n £isy- Then f or ever y /l ( A ) = E^Lo^A™ G C[A]®V and 
/ G V £/ie following identity holds: 

(1.3) Di(z)(h(d)f) = {Di(z){h(d)))f + h(z + d)(Di(z)f), 

where Di(z)(h(d)j is the differential operator obtained by applying D{(z) to the coefficients of 
h(d). 

Proof. In the case h(X) = A, equation (1.3) means (1.2). It follows that 
Di(z) od n = {z + d) n o Di(z) for every n G Z+ . 
The general case follows. □ 

As before, we denote by V® r C V r the subspace of all F = (i^)j g j with finitely many 
non-zero entries (r may be infinite), and introduce a pairing V r x V® r — > V/9V 

(1.4) (P,F)^JP-F, 
where, as before, J" denotes the canonical map V — ► V/dV. 
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Recall the definition (0.7) of the variational derivative ^ : V — > V®^. It follows immediately 
from (1.2) that 

(1.5) — od = o, iei, 

dUi 

i.e. dV C Ker 4-. In fact, we are going to prove that, under some assumptions on V, apart 
from constant functions, there is nothing else in Ker j-. 

First, we make some preliminary observations. Given / G V, we say that it has differential 

order n, and we write ord (/) = n, if 9 f n) / for some i E I and d / m) = for all j E I and 

du) du- 

m > n. In other words, the space of functions of differential order at most n is 

(1.6) V n := {/€V| ^yy = for all jel,m> n) . 

Clearly ord (/) = n if and only if / € V n \V n -i- We also denote by C C V the space of constant 
functions /, namely such that ®f m) = for all i G / and m G Z+, or, equivalently, such that 

du\ 

df = (this equivalence is immediate by (1.1)). We will also let C(uf) C Vo the subspace of 

functions / depending only on Uj, namely such that 9 f n) = unless j = i and n = 0. 

d Uj 

We can refine filtration (1.6) as follows. For i G I and n G Z + we let 

(1.7) V n ,i := {/ G V n = for all j > i\ C V„ . 

In particular, V n / = V n . We also let V n fi = V n -\ for n > 1, and Vo,o = C. 

It is clear by the definition (1.1) of d that if / = dg ^ and 5 G V n) j\V ni i_i, then 
/ G Vn+i^YVji+iji-i) an d in fact / has the form 

(1.8) f = J2 h j uf +1) +Y:h j u < f ) +r, 

j<i j>i 

where hj G V nj j for all j G /, r G V n j, and /ij 7^ 0. In particular, it immediately follows that 

(1.9) CndV = 0. 

Proposition 1.3. (a) The pairing (1.4) is non- degenerate, namely JP ■ F = /or every 
F G V® r ijf and on/y ifP = 0. 

(b) Moreover, let I be an ideal o/V® r containing a non-zero element F which is not a zero 
divisor. If j P ■ F = for every F in the ideal I, then P = 0. 

Proof. We can assume, without loss of generality, that r = 1. Suppose that P 7^ is such that 
j P ■ F = for every F G V. In this case, letting i 71 = 1, we have that P G <9V has the form 

(1.8). But then it is easy to see that ?4™ P does not have this form, so that f u\ n+ P 7^ 0. 
This proves (a). The argument for (b) is the same, by replacing P by PF. □ 
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We let g be the Lie algebra of all derivations of V of the form 
(1.10) X= Yl hn-^j , K n eV. 

By (1.1), d is an element of g, and we denote by g® the centralizer of d in g. Elements 
are called evolutionary vector fields. For X G g d we have X(uj n ^) = X(d n Ui) = d n X(ui), so 
that X = nez + ^C^l ) g f^l * s com pl e tely determined by its values X(ui) = Pi, i G /. We 
thus have a vector space isomorphism — 5 s , given by 

(1-H) V'BP={P i ) ieI ~X P = Y ^)A) £ a - 

The ^-tuple P is called the characteristic of the evolutionary vector field Xp. 
1.2 Normal algebras of differential functions. 

Definition 1.4. We call an algebra of differential functions normal if — f^y(V n i) = V n i for all 



i G /, n G Z + . Given / G V n ,i, we denote by J du^ f G V n ,i a preimage of / under the map 



3 



^ ^y. This antiderivative is defined up to adding elements from V n) i_i. 

Clearly any algebra of differential functions can be embedded in a normal one. 

We will discuss in Section 3 the exactness of the variational complex f2. In this language, 
Propositions 1.5 and 1.9 below provide algorithms to find, for a "closed" element in V, a 
preimage of d (up to adding a constant), and, for a "closed" element in Q 1 , a preimage of j^. 

Proposition 1.5. LetV be a normal algebra of differential functions. Then 

(1.12) Ker — = C ® dV . 

ou 

In fact, given f G V n i such that % = 0, we have d f n] G V n -i i, and letting 

f r (n-1) df 

g n -i,i= / du\ — t^t eVn-i,,, 
J du\ 

we have f — dg n -n G V n j_i. This allows one to compute an element g G V such that f = 
dg+const. by induction on the pair (n,i) in the lexicographic order. 

Proof. Equality (1.12) is a special case of exactness of the variational complex, which will be 
established in Theorem 3.2. Hence, if / G Ker 4- D V n ,i, we have / = dg + c for some g G V n -i,i 
and c G C. By equation (1.8), ®f n) G V n _i i for every i G I. The rest immediately follows by 

du) 



d.f 

the definition (1.1) of d. □ 



Example 1.6. In general (1-12) does not hold. For example, in the non-normal algebra V = 
C[u , u', u", . . . ] the element ^- is in Ker ^, but not in dV. But in a normal extension of V 
this element can be written as d\ogu. 
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In order to state the next proposition we need to introduce the Frechet derivative. 

Definition 1.7. The Frechet derivative Dt of / G V is defined as the following differential 
operator from V e to V: 

D f {d)P = X P (f) = £ 

One can think of Df(d)P as the differential of the function f(u) in the usual sense, that is, 
up to higher order in P, 

/(« + P) -/(«)« £ y$pM = D f {S)p, 

jei,nez+ ® u i 

hence the name Frechet derivative. 

More generally, for any collection F = (f a )a&A °f elements of V (where A is an index set), 
the corresponding Frechet derivative is the map Dp : V e — ► V" 4 given by 



(1.13) (D F (P))a = D fa (d)P [ = Xp(f a )), aeA. 

Its adjoint with respect to the pairing (1.4) is the linear map D F {9) : V®" 4 — ► V®^, given by 

(1.14) (D F (d)G) t = £ (~ d r{^G a ), iel. 

On the other hand, if F has only finitely many non-zero entries, namely F G V®" 4 C V" 4 , it is 
clear from the definition (1.14) that D* F can be extended to a map V" 4 — > V . In particular, 
for F G V® £ , both D F (d) and D* F (d) are maps from 1/ to V® £ . 

The following lemma will be useful in computations with Frechet derivatives. 

Lemma 1.8. (a) Both the Frechet derivative Dp{d) and the adjoint operator D* F {d) are 
linear in F G V" 4 . 

(b) Let M be an Ax B matrix with entries in V. We have, for F 6 V s , G 6 V® A and P G V*, 

D MF (d)P = £ F p {D Ma0 {d)P) + M af3 {D Fp (d)P) 
= (D M (d)P) -F + M- {D F (d)P) , 

and 

D* MF (d)G= Y D* Ma ^d){FpG a )+ D*F (d)(M a pG a ) 

a£A,/3eB a£A,f3£B 

= D* M {d)(G T ® F) + Dp{d)(G T M) . 

(c) Let F G V A , G G V®- 4 and P G V e . We have 

D 9F (d)P = d{D F (d)P) , D* dF (d)G = -D* F (d){dG) . 
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Proof. Part (a) is obvious. For part (b) the first equation follows immediately from the definition 
(1.13) of the Frechet derivative, while the second equation can be derived taking the adjoint 
operator. Part (c) can be easily proved using equation (1.2). □ 

We have the following formula for the commutator of evolutionary vector fields in terms of 
Frechet derivatives: 

(I- 15 ) [ X Pi X 0\ = X Dq(d)P-D P {d)Q ■ 

Elements of the form |^ G V®^ are called exact. An element F G V®^ is called closed if its 
Frechet derivative is a self-adjoint differential operator: 

(1.16) D* F (d) = D F (d) . 

It is well known and not hard to check, applying Lemma 1.2 twice, that any exact element in 
V® £ is closed: 

(1.17) D*si{d) = Dsi(8) for every // G V/dV. 

If V is normal, the converse holds. In fact, the next proposition provides an algorithm for 
constructing / G V such that = F, provided that F is closed, by induction on the following 
filtration of V . Given a triple (n,i,j) with n G Z + and i,j G I, the subspace (V® ) n i . C V®^ 

consists of elements F G V®^ such that F^ G V n j for every k £ I and j is the maximal index 
such that Fj G V n)i \V^_i. We also let (V®\, j = (V®\_ Mj and (V®\. = (V®\,_ M . 



Proposition 1.9. Let V be a normal algebra of differential functions. Then F G V® is closed 
if and only if it is exact, i.e. F = ^ /or some / G V. in /act, if F E (V®^) n f . is closed, then 
there are two cases: 

OF 

(a) if n is even, we have j <i and - (t 3 t) G V n /2j, and letting 



(1-18) /n,y = (-l) n/2 / ^" /2) | *4 n/2) £§y e V B/2>i , 

ffrj i/n is odd, we have j < i and — A- G V( n _i)/2i, and letting 

(1-19) = (-D ( " +1)/2 / < / du"^-^ G V (n+1)/2 , , 

J J du) ' 



we have F - %^ G (V®*) 



n,i,j—l' 



(Note that f n ,i,j in (1-18) is defined up to adding an arbitrary element g G V n / 2j j such that 
— 1§/2) £ Vn/2.j-ii an d similarly for f n ,i,j i n (1-19)J Equations (1.18) and (1.19) allow one 

to compute an element f G V such that F = x£ by induction on the triple (n,i,j) in the 
lexicographic order. 
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Proof. The first statement of the proposition is a special case of the exactness of the complex 
of variational calculus, which will be established in Theorem 3.2. Hence there exists / G V, 
defined up to adding elements from dV, such that F = 

Lemma 1.10. (a) Up to adding elements from dV , f can be chosen such that f G Vmk\Vmk-i 
and T~Tm) ^ V m ,fc\V m _i ,k f or some me1 + and k G I. 

(h) If f is chosen as in part (a), let (n,i,j) be the maximal triple (with respect to the lexico- 
graphic order) such that 

(1.20) — , f— 7-r + for some p G Z+ . 

1 3 

Then, there are two possibilities: 

A. n = 2m, j < i < k; 

B. n = 2m — 1, i > j = k. 

In both cases the only p G for which (1.20) holds is p = m. 
(c) g G {V m ) n4J \(V m ) nid _ v where (n,i,j) is as in (b). 

Proof. For (a) it suffices to notice that, if / G V m fc\V m k-i is such that d / m) G V m _i then 

f-dj d 4 m_1) ^py e V m , fe -i, 

and after finitely many such steps we arrive at the / that we want. 

Let (n,i,j) be as in (b). Notice that, since condition (1.20) is symmetric in i and j, we 
necessarily have i > j. Since / G V m ,fc, it is clear that (n,i,j) < (2m,k,k). Moreover, since, by 
( a )> d{ g V m k\V m -i fej it easily follows that (n, i,j) > (2m — 1, k + 1, fc) in the lexicographic 

order. These two inequalities immediately imply that either A. or B. occurs. In order to prove 
(c) we consider the two cases A. and B. separately. In case A. condition (1.20) exactly means 
that: 

1. -^y G V m ,i-i for every h > j, so that M- G VWi-i; 

dv> h ' oa h 

2. -^y G V m .i\V m ,i-i, so that # G V 2m ,i\V2m,i-l I 

3. -^y G V m/ i for every h < j, so that 4£r G V 2 m,i- 

By definition, the above conditions imply that |^ G (V®^) 2m 4 .. Case B. is treated in a similar 
way. □ 
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Returning to the proof of Proposition 1.9, let F = where / G V is chosen as in Lemma 
1.10(a). We then have 

9F ' - 9 5f - 8 / (-9r^ + (-sr- 1 ^T + -) = (-ir *' 



where, for the last equality, we used the commutation relation (1.2) and Lemma 1.10(b)-(c). 
The rest follows easily. □ 

In Section 3 we shall prove exactness of the whole complex Q in two different contexts. 
First, under the assumption that V is normal, and second, under an assumption involving the 

degree evolutionary vector field: 

(1-21) A=^ U f^, 

l<i<i 

for which we are going to need the following commutation rules: 

o or c 

(1.22) 9A = A9, — = (A + l)— pr , — A = (A + 1) — . 

Using this, one finds a simpler formula (if applicable) for the preimages of the variational 
derivative, given by the following well-known result. 

Proposition 1.11. Let V be an algebra of differential functions and suppose that F G V® e is 
closed, namely (1.16) holds. Then, if f G A _1 (u • F), we have 

(1.23) - F t G Ker (A + 1) , for all i £ I , 

but 

where u • F = Y^iel u iFi- 
Proof. We have, for i £ I, 

±-{u.F) = £ { - dr ° {ujFj)=Fi+ £ {-dT(f^u 3 

B F 

= F *+ E — ^9 n Uj = (A+m. 

In the third equality we used the assumption that F is closed, i.e. it satisfies (1.16), while in 
the last equality we used definition (1.21) of the degree evolutionary vector field. If we then 
substitute u • F = Af in the LHS and use the last equation in (1.22), we immediately get 
(1.23). □ 

Example 1.12. Consider the closed element F = (u' 3 , —u'2, —u'i) G R% 2 21 where R is the algebra 
of differential polynomials in m, U3. Since Ker (A + 1) = in R, we can apply Proposition 
1.11 to find that / = A" 1 ^!?^ — U2u' 2 ' — "U3^'i) = ^( UlU 3 — U2u' 2 ' — u^u^) satisfies F = 
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On the other hand, consider the localization V by U4 of the algebra of differential polynomials 
in ui, U2, U3, U4 and the following closed element F G (V 4 )2,2,2: 

Fx = u^ 2 u' 3 , F2 = 2u^ 3 u' 2 u' 4 — u^ 2 u' 2 ' , F3 = 2u^ 3 uiu' 4 — u^ 2 Ui , F4 = —u^ 3 u\u' s — u^ 3 (u 2 ) 2 ) ■ 

In this case A(u ■ F) = 0, hence formula (1.23) is not applicable, but we still can use the 
algorithm provided by Proposition 1.9, to find / = /2,2,2 + /i,4,3i where ^2,2,2 = \ u a 1 ( ? 4) 2 an d 
/i,4,3 = u^ 2 uiu' 3 . 

1.3 Poisson vertex algebras. 

Definition 1.13. Let V be a C[<9]-module. A X-bracket on V is a C-linear map V®V — > C[A]<8>V, 
denoted by / (8> g 1— ► {f\g}, which is sesquilinear, namely one has (/, g, h G V) 

(1-24) {df x g} = -\{f x g} , {fxdg} = (d + A){/ A <?} . 

If, moreover, V is a commutative associative unital differential algebra with a derivation d, a 
X-bracket on V is defined to obey, in addition, the left Leibniz rule 

(1-25) {f x gh} = {f x g}h + {f x h}g , 

and the right Leibniz rule 

(1-26) {fg x h} = {f x+d h}^g + {gx +9 h}^f . 

One writes {fxg} = Sngz + ^r(/(n)5)) where the C-bilinear products f^g are called the 
n-th products on V (n G and /( n )P = for n sufficiently large. In (1.26) and further 

on, the arrow on the right means that A + d should be moved to the right (for example, 

{fx +9 h}^g = J2nez+U(n) h )^lp-9)- 

An important property of a A-bracket {. x •} is commutativity (resp. skew-commutativity): 

(1-27) {g x f} = ^{f-x-89} (resp. = - J,f^ 9 g}) • 

Here and further the arrow on the left means that —A — d is moved to the left, i.e. <__{/_ ,\_a<7} = 

Snez+ \\ (/(n)ff) = e9lx {f-xg}- in case there is no arrow, it is assumed to be to the left. 
Another important property of a A-bracket {• \ •} is the Jacobi identity: 

(1-28) {fx{g,.h}} - {g u {f x h}} = {{fxg}x+„h} . 

Definition 1.14. A Poisson vertex algebra (PVA) is a differential algebra V with a A-bracket 
{• A •} : V ® V — > C[A] (g) V (cf. Definition 1.13) satisfying skew-commutativity (1.27) and Jacobi 
identity (1.28). A Lie conformal algebra is a C[9]-module, endowed with a A-bracket, satisfying 
the same two properties [K]. 

We next want to explain how to extend an arbitrary "non-linear" A-bracket on a set of 
variables with value in some algebra V of differential functions, to a Poisson vertex 

algebra structure on V. 
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Theorem 1.15. Let V be an algebra of differential functions, which is an extension of the 
algebra of differential polynomials R = C[uj n ^ \i E I,n G Z+]. For each pair i,j G /, choose 
{u ix uj} G C[A] &> V. 

(a) Formula 

(1.29) {f x g}= Yl 7%)( A + 9 ) re KA+a%K(-A-5r-^. 

defines a X-bracket on V (c/. Definition 1-13), which extends the given X-brackets on the 
generators Ui, i G /. 

(b) The X-bracket (1.29) on V satisfies the commutativity (resp. skew-commutativity) condi- 
tion (1.27), provided that the same holds on generators: 

(1.30) {ui\Uj} = ±^-{uj_ x _ d Ui} , for all i,j G I . 

(c) Assuming that the skew-commutativity condition (1.30) holds, the X-bracket (1.29) satis- 
fies the Jacobi identity (1.28) (thus making V a PVA), provided that the Jacobi identity 
holds on any triple of generators: 

(1.31) {ui\{uj^u k }} - {u jiM {u ix u k }} = {{u iX Uj} x+IM U k } , for alli,j,k G /. 

Proof. First, notice that the sum in equation (1.29) is finite (cf. Definition 1.1), so that (1.29) 
gives a well defined C-linear map {• a '} : V(£> V ^ <C[A] ® V. Moreover, for / = U{, g = Uj, i,j G 
/, such map clearly reduces to the given polynomials {ui\Uj} G C[A] (g> V. More generally, for 

/ = u[ m \ g = Uj , equation (1.29) reduces to 

(1.32) K ( "M n) } = ("A) m (A + dY{u ixUj } . 

It is also useful to rewrite equation (1.29) in the following equivalent forms, which can be 
checked directly: 

(1.33) {f xg }= Yl 7%( A + 5 n/w= E KA + a5M-A-dr-^. 



or, using (1.32), in the following further form: 



(/w= E (^^KM"'}^. 

where the parentheses indicate that d in the exponent is acting only on ~t™p but ^ ac ts on A. 
Here and further we use the following identity for /(A) G V[A], coming from Taylor's formula: 

f(d + X)u={e d ^u)f(X). 
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We start by proving the sesquilinearity relations (1.24). For the first one we have, by the 
second identity in (1.33) and the commutativity relation (1.2): 

{dfxg} = £ {u iX+9 gU(-\-dr(-^df 

Y {u iX+d g}^(-X-d) m ( df ^ +d-^— 
^ i 'waj^i i i ( m _i) ( m ) 

9/ 



where, as in (1.2), we replace - by zero for m = 0. Similarly, for the first sesquilinearity 

condition, we can use the first identity in (1.33) and (1.2) to get: 

ihdg} = Yl (i^)(A+sn/w 

jeJ,nez+ du- 

= ( A + 5 ) E ^( X + d) n {fxu j } = (X + d){f x g}. 

In order to complete the proof of part (a) we are left to prove the left and right Leibniz rules 
(1.25) and (1.26). For both we use that the partial derivatives — f^r are derivations of the 

product in V. For the left Leibniz rule, we use the first identity in (1.33) to get: 

Ux9h}= £ (h-% + 9-^)(X + dr{f xUj } = h{f x g} + g{f x h}, 

and similarly, for the right Leibniz rule, we can use the second identity in (1.33) to get: 

{fg X h}=^{u lX+d hU(-\-dr( K ^g+-^f)={f x+d hUg + {g x+d hUf. 

iei,mez+ d u i ® u i 

We next prove part (b). With the notation introduced in equation (1.27) we have, by (1.29): 

(1.35) ^{g- X - d f} = e d &{g- X f} = e d & £ 4L)("A + a) m {«i_ A+e «i}^(A - d)»-^ . 

i,jei " u i ° u j 

It immediately follows from (1.30) that {uj x+9 Ui}^F = ±e _S dX ^[ui X Uj}F^j , for arbitrary 
F € V. Equation (1.35) then gives 

JiS-x-Bf} = ± E ^(X + dT{u ix+dUj U(-\-d) m -^ = ±{f x g}, 
i,jei du ) du \ 



21 



thus proving (skew-)commutativity. 

We are left to prove part (c). First, observe that, since, by (1.31), Jacobi identity holds on 
any triple of generators Ui,Uj,Uk, applying (— X) m (— fi) n (X + fi + d) p to both sides of (1.31) and 

using sesquilinearity, Jacobi identity holds on any triple of elements of type u\ m ^ for i E I and 
m € Z + : 

(1.36) x {uf „ ««}} - {«<»> a 4 P) }} = {{^ m) A ^ ] } x+ / k p) } • 

We have to prove that for every f,g,h&V Jacobi identity (1.28) holds as well. We will study 
separately each term in equation (1.28). For the first term in the LHS of (1.28) we have, by the 
first identity in (1.33), combined with sesquilinearity (1.24) and the left Leibniz rule (1.25): 



- k GUI 



k£l,p€Z+ ULL k 

<'- 37 > = E {/^K4 P, )+ E ^S){/^»r>}} 

The fist term in the RHS of (1.37) can be rewritten, using again the first identity in (1.33), as 

MeJ duf>dv%> 

p,q&+ 

Notice that this expression is unchanged if we switch / with g and A with fj,, therefore it does 
not give any contribution to the LHS of (1.28). For the second term in the RHS of (1.37), 
we apply twice the second identity in (1.33), combined with sesquilinearity (1.24) and the left 
Leibniz rule (1.25), to get: 



Oh 



ij- keI dv> k du) du) 



^ - e ^(^^)(^*.%)(". ,m M»r,4"}} 

i,j,kei du k du l du ) 

m,n,pd'Lj r 

dh ( d d df u ( m ) / a d dg 



Furthermore, if we use again both identities (1.33) and sesquilinearity (1.24), we get the fol- 
lowing identity for the last term in the RHS of (1.39), 



(1.40) Y^(e d £^^\(e d %-^)}{uV 4 p) }=V{4"\ h}Jf x -°°\. 

i,j,keP u k OU i OU j j€l,n<EZ+ ° u j 
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The second term in the LHS of (1.28) is the same as the first term, after exchanging / with g 
and A with \i. Therefore, combining the results from equations (1.37)-(1.40), we get that the 
LHS of (1.28) is 

{AW} - W/aM} 



i,j,kel ° u k ° u i ° u j 



m,n 



(1.41) + ^ n) ^+M+S /l }-{/A-^|y}- E 

We finally use (1.36), and then the first equation in (1.33), to rewrite the first term in the RHS 
of (1.41) to get: 

(1.42) + £ {«f) A+M+s ^{/ A -%}- £ K M A + , + a^{^-^y}. 

We next look at the RHS of equation (1.28) with v in place of A + fi. By the first identity in 
(1.33), combined with sesquilinearity (1.24), and using the right Leibniz rule (1.26), we get: 

{{h9hh } = E {^R{M n) }^}= E {jh^ h }jf^y 

jei,nez + au j jei,nez+ au j 

(1-43) + £ (<<**)«/»«<■ >}.*} . 

We then expand {/a^j- } using the second identity in (1.33) combined with sesquilinearity 
(1.24), and then apply the right Leibniz rule (1.26), to rewrite the last term in the RHS of 

(1.43) as 



E (<*:&)«/»«, w K*> 

E/ 84- d 9 \ ffd4r df \ f ( m ) ( n ), \ 



i,jei du ) du ^ 

m,nGZ_|_ 

dg \( d d 9 d df 



- m,nGZ_|_ 



( d d dg \( ( d d df \ "i (m) (n ) 



m,nGZ+ 
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Furthermore, the last term in the RHS of equation (1.44) can be rewritten, using equation 
(1.33) combined with sesquilinearity (1.24), as 

i,j€l uu j ua i 



iii.ii 



(1.45) = ]T {-^y, +a /l}^^K (m) A ^_ a5 } = - J] {-^^g/tj^^A^H}. 

For the last equality, we used that, by part (b), the A-bracket (1.29) satisfies the skew- 
commutativity (1.27). We can then put together equations (1.43)-(1.45) with v = A + /x, 
to get the following expression for the RHS of equation (1.28), 



(1.46) + {-r^j*+»+d h }^{f* u< r ] }- E {-M^^+d 11 }^^^}- 

To conclude, we notice that the first term in the RHS of (1-42) coincides with the first term in 
the RHS of (1.46). Moreover, it is not hard to check, using (1.33) and the fact that - ® m) and 

— ^-r commute, that the second and third terms in the RHS of (1-42) coincide, respectively, 

with the second and third terms in the RHS of (1.46). □ 

Theorem 1.15(a) says that, in order to define a A-bracket on an algebra of differential 
functions V extending R = C[u-™^ \ i £ I, n £ Z+], one only needs to define for any pair i,j £ / 
the A-bracket 

(1.47) {u ixUj } = HjiiX) £ V[A] . 

In particular A-brackets on V are in one-to-one correspondence with £ x ^-matrices H(X) = 
(Hij(\)) i j eI with entries Hij(X) = X^nLo ^ij;n^ n hi V[A], or, equivalently, with the corre- 
sponding matrix valued differential operators H(d) = (ffjj(<9)) i . j : V® e — ► V ■ We shall 
denote by {• \ -}h the A-bracket on V corresponding to the operator H{d) via equation (1.47). 

Proposition 1.16. Let H(d) = [Hij(d)) i . gJ be an £ x £ matrix valued differential operator. 

(a) The A-bracket {■ \-}h satisfies the (skew-)commutativity condition (1.27) if and only if 
the differential operator H{d) is self (skew-) adjoint, meaning that 

N 



(1.48) Htj{d) := £(-3) n o H ji;n = ±i^(<9) . 



n=0 
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(b) If H(d) is skew-adjoint, the following conditions are equivalent: 

(i) the X-bracket { -\- }h defines a PVA structure on V , 

(ii) the following identity holds for every i,j, k S /: 



E fip^ + «r**w - + sr*«oo) 

(1.49) heI,n£Z + \ au h 0u h / 

hei,n&z + 9u h 
(Hi) the following identity holds for every F,G€ y®^.- 

H(d)D G (d)H(d)F + H(d)D* H{9)F (d)G - H(d)D F (d)H(d)G 

(1.50) +H(d)D* F (d)H(d)G = D H{d)G (d)H(d)F - D H(d)F (d)H(d)G . 

Proof. Equation (1.30) is the same as equation (1.48), if we use the identification (1.47). Hence 
part (a) follows from Theorem 1.15(b). Next, equation (1.31) is easily translated, using the 
identification (1.47) and formula (1.29), into equation (1.49). Hence, the equivalence of (i) 
and (ii) follows immediately from Theorem 1.15(c). We are thus left to prove that, if H{d) is 
skew-adjoint, then conditions (ii) and (iii) are equivalent. Written out explicitly, using formulas 
for the Frechet derivatives in Section 1.2, equation (1.50) becomes (k £ /): 

E *«<»> (f% + C^fW^o,) 

i,j,hei \du)' x J V du\ J 

(1.5D -^^H ih{d) G h ) 4- { - BT (^{H jh{d )G h ))} 



J 



E ( - (a» g , j( a )gj )g^P 



i,j,h£l 



We then use Lemma 1.2 to check that the second term in the LHS of (1.51) is 

(i.52) £ jfc(fl)H»» ( (^#^) e, + |py (^0)0,; ' 



i,j,h£l 



the first term in the RHS of (1.51) is 



(1.53) £ (>^ W .) fe^l + £ ^(a)f|%(a^ ( a)F,)) , 

nGZ+ n€Z+ 

and the second term in the RHS of (1.51) is 

(1.54) - £ (Vi^)^) f^^W) - £ H fc ^)(^(^(^ 
i,j,h<=I \ ® u i J i,j,hel \^ u i 
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We then notice that the second term in (1.53) cancels with the first term in the LHS of (1.51), 
the second term in (1.54) cancels with the third term in the LHS of (1.51) and, finally, the second 
term in (1.52) cancels with the last term in the LHS of (1.51), since H{d) is skew-adjoint. In 
conclusion equation (1.51) becomes 



(1.55) 




Since the above equation holds for every F,G G V®^, we can replace d acting on Fi by A and d 
acting on Gj by \i and write it as an identity between polynomials in A and /i. Hence (1.55) is 
equivalent to (1.49). □ 

Definition 1.17. A matrix valued differential operator H(d) = (Hij(d)). which is skew- 
adjoint and satisfies one of the three equivalent conditions (i)-(iii) of Proposition 1.16(b), is 
called a Hamiltonian operator. 

Example 1.18. The Gardner-Faddeev-Zakharov (GFZ) PVA structure on R = C[u, u', u" , . . .} 
is defined by 

(1.56) {u x u} = Al 

(further on we shall usually drop 1). In fact, one can replace A in the RHS of (1.56) by any 
odd polynomial P(X) € C[A], and still get a PVA structure on R. Indeed, the bracket (1.56) 
is skew-commutative and it satisfies the Jacobi identity for the triple u, u, u, since each triple 
commutator in the Jacobi identity is zero. 

Example 1.19. The Virasoro-Magri PVA on R = C[u, u',u", ...], with central charge c G C, is 
defined by 

(1.57) {u x u} = (d + 2X)u + A 3 c . 

It is easily seen that the bracket (1.57) is skew-commutative and it satisfies the Jacobi identity 
for the triple u,u,u. 

Example 1.20. Let q be a Lie algebra with a symmetric invariant bilinear form ( . | . ), and let 
p be an element of g. Let {ui}i<zj be a basis for q. The affine PVA associated to the triple 
(fl? (" I ")>p) i s the algebra R = C[u^ \ i G I,n G Z+] together with the following A-bracket 

(1.58) {a x b} = [a, b] + (p\ [a, b}) + \(a\b) , a, b G . 

Note that, taking in the RHS of (1.58) any of the three summands, or, more generally, any 
linear combination of them, endows R with a PVA structure. 

Note that Example 1.18 is a special case of Example 1.20, when q = Cu is the 1-dimensional 
abelian Lie algebra and (u\u) = 1. 

The following theorem further generalizes the results from Theorem 1.15, as it allows us to 
consider not only extensions of R, but also quotients of such extensions by ideals. 
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Theorem 1.21. Let V be an algebra of differential functions, which is an extension of the 
algebra of differential polynomials R = C[u\ n ' \i G I,n G Z+]. For each pair i,j G /, let 
{uixUj} = Hji(X) G C[A] ® V, and consider the induced X-bracket on V defined by formula 
(1.29). Suppose that J C V is a subspace such that dJ C J, J ■ V C J, {J\V} C C[A] <8> 
J, {V\J} C C[A] (8) J, and consider the quotient space V/J with the induced action of d, the 
induced commutative associative product and the induced X-bracket. 

(a) The X-bracket on V/J satisfies the commutativity (resp. skew-commutativity) condition 
(1.27), provided that 

(1.59) {ui\Uj} =F <-{uj _x_gUi} G C[A] (8) J , for all i,j G I . 

(b) Furthermore, assuming that the skew-commutativity condition (1.59) holds, the X-bracket 
on V/J satisfies the Jacobi identity (1.28), thus making V/J a PVA, provided that 

(1.60) {uixiuj^Uk^-iuj^UixUk^-HuixUj}^^} G C[A,//]<g)J , for all i,j,k G I . 

Proof. By the assumptions on J, it is clear that the action of d, the commutative associative 
product and the A-bracket are well defined on the quotient space V/J. It is not hard to check, 
using equation (1.35) and the assumption (1.59), that {fxg} T*- {d-x-df} £ C[A]J for all 
f,g G V, thus proving part (a). Similarly, if we compare equations (1.41) and (1.46) and use 
the assumptions (1.59) and (1.60), we get that {fx{g^h}}-{g^{fxh}}-{{fxg}x+^h} G C[X, fi]J 
for all f,g,h G V, thus proving part (b). □ 

Example 1.22. Let A be a Lie conformal algebra with a central element K such that dK = 0. 
Then V k (A) = S(A)/(K — kl), k G C, carries the usual structure of a unital commutative 
associative differential algebra endowed with the A-bracket, extending that from A by the left 
and right Leibniz rules, making it a PVA. This generalization of Examples 1.18, 1.19, 1.20, may 
be viewed as a PVA analogue of the Lie-Kirillov-Kostant Poisson algebra S(q), associated to a 
Lie algebra q. 



1.4 The Beltrami A-bracket. All the examples of A-brackets in the previous section were 
skew-commutative. In this section we introduce a commutative A-bracket, called the Beltrami 
A-bracket, which is important in the calculus of variations. 

Definition 1.23. The Beltrami X-bracket is defined on the generators by, 
(1.61) {u ix Uj}B = b~ij , 

By Theorem 1.15, this extends to a well defined commutative A-bracket on any algebra of 
differential functions V, given by the following formula, 

a- 62 ) if^}B= £ ( -i r -^(A + a)^-^. 

iei 9u\ du\ 
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Both the variational derivative j- and the Frechet derivative Dp(d) have nice interpretation 
in terms of this A-bracket. Indeed we have 

(1.63) $L = {fxUi} B \ x=0 , D F (\) i - = {u jx F i } B . 

Moreover the whole A-bracket {fx u i} B c & n be interpreted in terms of the well-known higher 
order Euler operators E^ m) = £ ngZ (™)(-l) n d n ~ m 9 • 



du 



{f^} B = £ ^f- 

n£Z+ 

The Beltrami A-bracket does not satisfy the Jacobi identity (1.28). Indeed, when restricted 
to the subalgebra of functions, depending only on itj, i G J , the Beltrami A-bracket becomes the 
bracket: {f,g} B = E i£i ~8ul cHm' wrj i cri does not satisfy the Jacobi identity. The latter bracket 
is one of the well-known Beltrami operations in classical differential geometry, hence our name 
"Beltrami A-bracket" . 

On the other hand we have the following two identities for the triple A-brackets, which can 
be easily checked using (1.62) together with the sesquilinearity and the Leibniz rules: 

(1-64) {Ui\{f^g}B}B ~ {fti{Ui X g}B}B = {{uixf}B X+ ^9}B , 

(1-65) {fx{g^Ui} B }B + {^{/a^IbIb = {{hg}B X +n u i}B ■ 

Both Propositions 1.5 and 1.9 can be nicely reformulated in terms of the Beltrami A-bracket. 
Namely, equation (1.12) states that / E dV + C if and only if {fx u i}B\ x _ = 0, for all i £ /, 
where the only if part is immediate from sesquilinearity (1.24). Similarly, Proposition 1.9 states 
that F = {fxUi}s\ x=0 for some / G V if and only if {uj x Fi}s = {Fj x Ui}B, for all i, j G /. Again 
the only if part follows immediately by identity (1.65) with g = Uj and A = 0. 

Even though the Beltrami A-bracket does not make V into a PVA, it has the following nice 
interpretation in the theory of Lie conformal algebras. Consider the Lie conformal algebra 
A = (Bi£iC[d]iii with the zero A-bracket. Then V is a representation of the Lie conformal 
algebra A given by 



(1-66) u ix f = {u ix f}B ■ 

This is indeed a representation due to equation (1.64) with / = Uj, g G V. We will see in Section 
3 that, in fact, the whole variational complex can be nicely interpreted as the Lie conformal 
algebra cohomology complex for the Lie conformal algebra A and its module V. 

1.5 Hamiltonian operators and Hamiltonian equations. Lie conformal algebras, and, 
in particular, Poisson vertex algebras, provide a very convenient framework for systems of 
Hamiltonian equations associated to a Hamiltonian operator. This is based on the following 
observation: 

Proposition 1.24. Let V be a C[d]-module endowed with a X-bracket { \ } : V®V ^ C[A] <S> V, 
and consider the bracket on V obtained by setting A = 0: 

(1-67) {f,g} = {fxg}\x=o, f,ge V. 
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(a) The bracket (1.67) induces a well defined bracket on the quotient space V/dV. 

(b) If the X-bracket on V satisfies the Lie conformal algebra axioms (1-27) and (1.28), then 
the bracket (1-67) induces a structure of a Lie algebra on V/dV, and a structure of left 
V / dV -module on V. 

(c) If V is a PVA, the corresponding Lie algebra V/dV acts on V (via (1.67)) by derivations 
of the commutative associative product on V, commuting with the derivation d, and this 
defines a Lie algebra homomorphism from V/dV to the Lie algebra of derivations ofV. 

Proof. Part (a) is clear from the sesquilinearity conditions (1-24) for the A-bracket. Clearly, 
the bracket (1.67) on V satisfies the Jacobi identity for Lie algebras (due to the Jacobi identity 
(1.28) for Lie conformal algebras with A = fi = 0). Moreover, the skew-commutativity of the 
bracket on V/dV induced by (1-67) follows from the skew-commutativity (1.27) of the A-bracket. 
Furthermore, it is immediate to check that, letting {fh, u} = {h\u}\\ = o for fh = h+dV G V/dV 
and u G V, we get a well-defined left V/<9V-module structure on V, proving (b). Claim (c) follows 
from the left Leibniz rule and the Jacobi identity with A = fi = 0. □ 

Proposition 1.24 motivates the following definition and justifies the subsequent remarks. 

Definition 1.25. (a) Elements of V/dV are called local functionals. Given / 6 V, its image 
in V/dV is denoted by f f. 

(b) Given a local functional fh G V/dV, the corresponding Hamiltonian equation is 

(1.68) ^- = {h x u}\ x=0 = {fh,u} 

and {fh,-} is the corresponding Hamiltonian vector field. 

(c) A local functional J f G V/dV is called an integral of motion of equation (1.68) if -± = 
mod dV, or, equivalently, if 

{fhjf} = 0. 

(d) The local functionals fh n , n G Z + are in involution if {jh m , fh n } = for all m,n£ Z+. 
The corresponding hierarchy of Hamiltonian equations is 

du 

— = {h nX u}\x=o =: {Jh n ,u} , ti£Z + . 
ctt n 

(Then all fh n 's are integrals of motion of each of the equations of the hierarchy.) 

From now on we restrict ourselves to the case in which the PVA V is an algebra of dif- 
ferential functions in the variables {nj}j e / (cf. Definition 1.1). In this case, by Theorem 
1.15, the A-bracket {-\-}h is uniquely defined by the corresponding Hamiltonian operator 
H(d) = (Hij(d)) . - eJ (cf. equation (1.47)), and the A-bracket of any two elements f,g G V can 

be computed using formula (1.29). In this case the Hamiltonian vector field {fh, •} is equal to 
the evolutionary vector field X H(d -,sh (defined by (0.6)) and the Hamiltonian equation (1.68) 
has the form 
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where |^ = G V®^ is the variational derivative of h (defined by (0.7)). Moreover, the 

corresponding Lie algebra structure of V/dV is given by 

Remark 1.26. Since Jh i— ► X H ^sh is a Lie algebra homomorphism, local functionals in invo- 
lution correspond to commuting evolutionary vector fields. If a sequence Jh n £ V/dV is such 
that G V®^ span infinite-dimensional subspace and dimKer H(d) < oo, then the vector 
fields X H ^j sh n span an infinite-dimensional space as well. 

Definition 1.27. The Hamiltonian equation (1.68) is called integrable if there exists an infinite 
sequence of local functionals Jh n , including Jh, which span an infinite-dimensional abelian 
subspace in the Lie algebra V/dV (with Lie bracket (1.67)), and such that the evolutionary 
vector fields X H , Q s sh n span an infinite-dimensional space (they commute by Remark 1.26). 



1.6 Center of the Lie algebra V/dV. Consider a Hamiltonian operator H (d) = (Hij(d)) . . g 

V® £ — ► V^i where V is an algebra of differential functions extending R = C[uj n ' ) | i G /, n G 
and consider the corresponding PVA structure on V, defined by equations (1-29) and (1.47). 
Recall that two local functionals J f and f g are said to be in involution if {f f, fg} =0. We 
have, in fact, three different "compatibility conditions" for the local functionals ff,fg€L V/dV 
or the corresponding Hamiltonian vector fields X H , 9 ,s£, X H(d ,sg_, namely: 

(iii) / [X Hm si,X H(d) sg] (h) = for every h G V. 
By definition the action of the Hamiltonian vector field X H( ^si on V is given by 




hence condition (ii) above is equivalent to saying that {ff,fg} is in the kernel of the rep- 
resentation of the Lie algebra V/dV on the space V (see Proposition 1.24(b)), namely to the 
equation H(d)-^{f f,Jg} = 0, while condition (iii) is equivalent to saying that {//, fg] lies 
in the center of the Lie algebra V/dV. In particular, we always have the implications (i) ==>■ 
(ii) ==^ (iii). The opposite implications, in general, are not true. To understand when they 
hold, we need to study the kernel of the representation of the Lie algebra V/dV on the space V 
and the center of the Lie algebra V/dV. Clearly the former is always contained in the latter, 
and we want to understand when they are both zero. 

Lemma 1.28. Let i £ I and let f G V be a non-zero element such that 
(1.71) uff G dV , for every n G Z + . 
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Then ord (/) = 1, and in fact f = da(ui), for some a(ui) G C(uj) C V. In particular, if £ > 2 
and condition (1.71) holds for every i £ I and n G Z + , i/ien / = 0. 

Proof. Let / G V\{0} satisfy (1.71) for some i £ I. If we let n = in (1.71) we get that f = dg 
for some g G V. Let ./V = ord (g) (so that ord (/) = N + 1). 

Let us consider first the case iV = 0, namely g G Vo- For f = 1 we have i = 1 and / = dg 
with g G Vo = C(ui), which is what we wanted. For i > 2, we have from condition (1.71) that, 
for every n G Z+, there exists g n G Vo such that 



u^dg = dg n , 

9w'. 



and, if we apply j G I, to both sides of the above equation, we get, using (1.2), that 



■ j 



Usin s that sgiSj = e J ' we get from (L72) that 

9 = ±U dg 



duj V 1 <9tij / dui V * 3uj 

which gives J^- = for every j G I different from i. Hence g G C(nj), as we wanted. 

Next, we will prove that, for N > 1, no element / G V of differential order N + 1, satisfying 
condition (1.71) for every n G Z+, can exist. By assumption (1.71), for every n > 1 there exists 
an element g n G V of differential order iV such that 



(1.73) ufdg = dg n . 

iju 



If we apply to both sides of the above equaiton we get, using (1.2), that 
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do 

n > 1 



which is equivalent to saying, after performing integration in 1*^^, j G /, that the elements 
g n G R are of the form 

9n = u"g + r n , 

for some v n G Vjv— l* Substituting back in equation (1.73) we thus get 

ufdg = a«g + r„) , 

which immediately gives 

(1.74) u n i9 u\ = --L^dr n+1 G dV , 

for every n G Z + . Suppose first that N = 1. Then ord (g) = 1 and equation (1.74) with n = 
implies that gu^ G dV. But this is is not possible since gu' { , as function of u^, . . . ,u'«, is not 
linear, namely gu[ does not have the form (1.8) with N = 0. Next, suppose, by induction, that 
no element in Vat exists which satisfies condition (1.71) for every n G Z+, and let / G Vat+i\Vjv, 
be as above. Then equation (1.74) exactly says that gu'^ which has differential order TV, also 
satisfies assumption (1.71), thus contradicting the inductive assumption. □ 
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Remark 1.29. Consider the pairing VxV-> V/dV given by (1.4) with r = 1. Lemma 1.28 can 
then be stated in terms of the following orthogonality condition: 

(1.75) (CM) X C dC{ Ui ) 

and the above inclusion becomes an equality under the additional assumption that : C(ui) — > 
C{ui) is surjective. 

Lemma 1.28 will be useful in Section 2.2 to prove that the integrals of motion of the KdV 
hierarchy span a maximal abelian subalgebra of the Lie algebra R/dR, where R = C[u^ \ n G 
Z + ], with Lie bracket (1.70) with H(d) = d (cf. Theorem 2.15). An analogous result will be 
needed to prove that the integrals of motion of the "linear" KdV hierarchy span a maximal 
abelian subalgebra. This is given by the following 

Lemma 1.30. Let i £ I and let df G dV be an element such that 

(1.76) u? df G dV , for every n G Z + . 

Then df is a C-linear combination of the monomials u^ 2n+1 \ n G Z + . In particular, if £ > 2 
and condition (1.76) holds for every i G I, then df = 0. 

Proof. We start by observing that, under the assumption (1.76), / must be a function only of 

df 

statement to the case £ = 1. Indeed suppose, by contradiction, that, for some j ^ i and some 
du(N) + and Qu[ 

£(^a/) = -B^($«r i »)=o. 

J p=0 j 

Here we are using the fact that, in any algebra of differential functions V, dV C Ker 4-. If we 
then take n large enough (such that 2n > ord (/) > N), by looking at the terms of highest 
differential order, the above equation has the form 

(-l)" +1 ^«^ +1) +r = 0, 
du)j ' 

where r G V2n+N- If we then apply — (2n ^ JV+1) to both sides, we get d f N) = 0, contradicting 

1 3 

our original assumption. 

Next, we consider the case £ = 1. We need to prove that if / G V is such that 



Ui and its derivatives, namely — = 0, for all j / i,n 6 Z+. In other words, we reduce the 

du) 

lee 

N > 0, we have 9 / N) ^ and d f n) = Ofor all n> N. As a consequence of (1.76) we get that 

1£ ■ 



(1.77) fu {2n+1) G 9V , for all n G 



then 9/ G Span c {u( 2m+1 ), m G Z+}. 

It is clear from (1.77), that / has even differential order. Indeed, if ord (/) = 2n + 1 is odd, 
u (2n+i) j j g no j. Q £ f orm (1.8) and hence it cannot be in dV, contradicting (1.77). We thus 
let ord (/) = 2iV and we will prove the statement by induction on N > 0. 
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Let us consider first the case N = 0, namely / = f(u) G Vo = C(u). If we let n = 1 in 
(1.77), we have u"'f(u) G dV, so that 



'/) = l^-5 3 (/(n)) =0. 



Expanding 9 3 (/(n)) and looking at the coefficient of u", it follows that = 0, namely 
9/ = qm' for some a G C, as we wanted. 

Next, let iV > 1. Condition (1.77) says that u'f G <9V, hence, by (1.8), / must have the 
form 

(1.78) f = au^ N '> + b, 

where a,b G V27V-1. We want to prove that, necessarily, a G C and b G V2N-2 (so that, using 
the inductive assumption, we can conclude that df G Span c {?/, u^ 3 \ . . . , u^ 2N+1 ^}). Condition 
(1.77) with n = N gives u^ 2N+l ^ f G dV. Using (1.78) and integrating by parts we have 



(1.79) \(da)(4 2N ^f + (db)u^ G 8V . 



In particular the LHS of (1-79) must have the form (1.8). It follows that the coefficients of 
(^( 2A 0) an( j ^(2^)^ j n (1.79) must vanish, namely 

(L80) cd^) =0 < l da+ dJ^ = °- 

The first of the above equations says that a G V 2 n-2- It follows that da, viewed as a function 
of uC 2 ^- 1 ), is linear, and therefore, by the second equation in (1.80), it follows that b, viewed 
as a function of u( 2Ar_1 \ has degree at most 2, namely 

(1.81) b = a + Pu( 2N -V + I^u^" 1 )) 2 ,-lda = + 7 u^ 2N ^ , 

where a, /3, 7 G V 2 n-2- Next, let us apply condition (1.77) with n = N — 1, namely u^ 2N ~^ f G 
dV. By equation (1.78) and integration by parts, this condition gives 

A {da) ^N-i ) ^ + bu{2N . 1) G 5V) 

which, by (1.81), is equivalent to 

(1.82) cm^" 1 ) + 2(3( U ( 2N -^) 2 + ^(u^ 1 )) 3 G dV . 

To conclude we observe that the LHS of (1.82), being in dV, must have the form (1.8), and this 
is possible only if (3 = 7 = 0, which is exactly what we wanted. □ 

Remark 1.31. Lemma 1.30 can be stated, equivalently, in terms of the pairing V x V — > V/dV 
given by /, g ^ J fg, as an orthogonality condition similar to (1.75). Let E = Span c {l ; u[ 2n \ n G 
Z_|_}. Then: 

(1.83) E L = dE 
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We next use Lemma 1.28 and Lemma 1.30 to study the center of the Lie algebra V/dV. 

Corollary 1.32. (a) If ' I > 2, the center of the Lie algebra V/dV with Lie bracket (1.70) coin- 
cides with the kernel of the Lie algebra representation ofV/dV on V defined in Proposition 
1.24(b), namely it consists of the local functionals J f G V/dV such that 

(1.84) H{d) 5 4- = 0. 

ou 

In particular, the "compatibility conditions" (ii) and (Hi) above are equivalent. 

(b) If £ = 1, the center of the Lie algebra V/dV with Lie bracket (1.70) consists of the local 
functionals Jf G V/dV such that 

(1.85) H(d)^- = au , where aeC. 

ou 

Proof. Suppose that the local functional J f is in the center of the Lie algebra V/dV. Then 
necessarily, for every i £ I and n G Z+, 

jei 3 

We can thus apply Lemma 1.28 and Lemma 1.30 to deduce that, if t > 2, then (1.84) holds, 
while, if £ = 1, then 

(1.86) H(d) 5 4- G dC(u) n Span c {u( 2n+1 ) | n G Z+| = {au' \ a G C} . 

ou 

To conclude we just notice that, integrating by parts, fir-u' = Jdg = for every g G V, so 
that, if condition (1.85) holds, then J f lies in the center of the Lie algebra V/dV. □ 

In conclusion of this section, we discuss two examples: in Proposition 1.33 we consider the 
GFZ PVA from Example 1.18, for which the center of the Lie algebra R/dR is strictly bigger 
than the kernel of the representation of R/dR on R, and in Example 1.35, we discuss a simple 
case in which the "compatibility conditions" (ii) and (iii) hold, but (i) fails. 

Proposition 1.33. Consider the Hamiltonian operator H(d) = d on R = C[u^ \ n G Z+], so 
that the corresponding Lie algebra structure on R/dR is given by the bracket 



I Mr 

and the representation of the Lie algebra R/dR on the space R is given by 



{//,,} = *„„(,)= £(ff~Vv* 



5uJ <9u( n ) 
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(a) The kernel of the Lie algebra representation of R/dR on R is two-dimensional, spanned 
by 

fl, fu G R/dR. 

(b) The center of the Lie algebra R/dR is three-dimensional, spanned by 

fl, fu, fu 2 G R/dR. 

(c) The "compatibility" conditions (i), (ii) and (Hi) for the local functionals ff,fg€ R/dR 
are equivalent. 

Proof. A local functional ff is in the kernel of the representation of R/dR on R if and only if 
dj£ = 0, which is equivalent to saying that ff G Span c { J 1, fu}. For the last assertion we use 
the fact that Ker 4- = C + dR (see Proposition 1.5). Similarly, by Corollary 1.32(b), a local 
functional ff G R/dR is in the center of the Lie bracket if and only if d^ t G Cu', which, by 
the same argument as before, is equivalent to saying that ff G Span c { f 1, fu, f u 2 }. 

We are left to prove part (c). We already remarked that (i) implies (ii), which implies (hi), 
and we only need to prove that (hi) implies (i) . Recall that condition (hi) is equivalent to saying 
that {f f, fg} lies in the center of the Lie algebra R/dR. Hence, by part (b) we have 

(1.87) J ^d^ G Span c {/1 ,fu,fu 2 } C R/dR. 

On the other hand, R admits the decomposition 

(1.88) R = C[u]®R 9 , 

where R 9 C R is the ideal (of the commutative associative product) generated by the elements 
u( n \ n > 1, and clearly dR C R 9 . Hence, by (1.87), we get 

6 -^-d 5 4- G R 9 n (C[d + dR) = dR , 
ou ou v ' 

which implies fi^d^- = 0, as we wanted. □ 

Remark 1.34. By the same argument as in the proof of Proposition 1.33, it is not hard to show 
that, if V = R and H(d) = (Hij(d)) : R® e — > R^ is an injective Hamiltonian operator (or, 

more generally, if Ker H (d) film 4- = 0), then the commutativity condition of Hamiltonian vec- 
tor helds [X H(d sSf_ , ^ H /g) s g ] = 0, is equivalent to the following condition for the corresponding 

local functionals J f, fg G R/dR: 

(1.89) {ffjg} G C/l. 

Moreover equation (1.89) is clearly equivalent to {f f, fg} = if the operator H(d) does not 
have a constant term (by taking the image of both sides of equation (1.89) via the quotient 
map R/dR -» fC). 
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Example 1.35. Take R = C[p( n \ q( n \ \ n G Z+], and define the A-bracket on R by letting z 
to be central and 

{p x q} = z = -{q\p} , {p x p} = ~{q\q} = . 

Let H (<9) be the corresponding Hamiltonian operator. In this example we have I = 3, therefore 
by Corollary 1.32(a) we know that the kernel of the representation of R/dR on R and the center 
of the Lie algebra R/dR coincide. In particular conditions (ii) and (hi) above are equivalent. 
Next, take the local functionals Jp,Jq G R/dR. By the above definition of the A-bracket on 
R we have {Jp, Jq} = / 0, namely condition (i) fails, but Jz is in the center of the Lie 
algebra R/dR, namely condition (iii) holds. 

1.7 Poisson vertex algebras and Lie algebras. In this section we shall prove a converse 
statement to Proposition 1.24(b). In order to do so, we will need the following simple lemma: 

Lemma 1.36. Let V be an algebra of differential functions in the variables {«i}j g /. 

(a) IfPeV e is such that J = for every ff e V/dV, then P = 0if£>2, and P G Cu^ 
if 1 = 1. 

(b) IfH{d) = {Hij(d)). . eI and n G 1+ are such that H(d)^ G V n ® 1 for every Jfe V/dV, 
then H = 0. 

Proof. Note that -^j--^jpi = dij 11 ?, and ^ { u t^) 2 ) = &i,j u ^? n ' '■ Hence, part (a) is an 

immediate corollary of Lemmas 1.28 and 1.30. For part (b), suppose that Hij(d) = X^=o hijnd 71 
is a differential operator such that hijN ^ 0. Let then h = — j — (nj A/ ^) 2 , so that -j^ = 
S kjuf M) . If M is sufficiently large, we have ( £ +N) {H(d)j^) i = h ijN / 0. □ 

Proposition 1.37. Let V be an algebra of differential functions in the variables {ui}i^i, en- 
dowed with a X-bracket {-a - }- Suppose that the bracket on V/dV induced by (1.67) is a Lie 
algebra bracket. Then the X-bracket satisfies skew-commutativity (1.27) and Jacobi identity 
(1.28), thus making V into a PVA. 

Proof. Recall that the A-bracket on V is given by (1.29), hence, the induced bracket on V/dV 
is given by (1.70) where Hji(d) = {uigUj}^. Integrating by parts, the skew-symmetry of this 
bracket reads 

fl (H{d) + H . {d)) s± = , 

for every Jfjg G V/dV. By Lemma 1.36(a), it follows that (H (d) + H*(d))^ G Vf* for 
every ff G V/dV. Hence, by Lemma 1.36(b), we have H*(d) = —H{d). Recalling (1.48), this 
identity can be equivalently written in terms of the A-bracket as {ut\Uj} = —{uj_ x _QUi}, which 
implies the skew-symmetry of the A-bracket by Theorem 1.15(b). 

Next, let us prove the Jacobi identity for the A-bracket, which, by Theorem 1.15, reduces 
to checking (1.31). We will use the Jacobi identity for the Lie bracket: 

(1-90) {//, {Jg, Jh}} - {jg, {ff, jh}} = {{//, Jg}, jh} . 



36 



Recalling that the bracket (1-70) is the same as (1-29) with A = 0, we have 



i,j,P,qeI J 9u j du Q ° Up ° Ul 



Note that the first term in the RHS is symmetric in / and g, hence it cancels out in the LHS 
of (1.90). By Lemma 1.2, the last term in the RHS of (1.91) is 



-»K(4^)H'<)) 

J 



Integrating by parts, we can replace X]nez + Q w ^" by the variational derivative Hence 
the LHS of (1.90) reads, 

Similarly, the RHS of (1.90) reads 

Ej^X-ar((^)(^a)|9 

5g (dHji{d) 8f \ t ( d Sf\( TTmta , 5g 



(L94) + 5 Uj { du f) 6uJ + \ du (n) 6uJ { H ^ 9) 5 Uj 

where in the last term we used Lemma 1.2. The third term in (1.93) is equal to the first term 
in (1.94), due to (1.17). Similarly, the last term in (1.93) is equal to the last term in (1.94), due 
to skew-adjointness of H(d). Hence, the Jacobi identity (1.90) is equivalent to the following 
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relation 



nGZ+ 

V Sg fdH pi (d) 6f 



Since this equation holds for every Jf,Jg,Jh G V/dV, we conclude, using Lemma 1.36, that 

£ fer<* + *>***<*) - ^f#(" + »>***oo 



-fl-^A + Zi + aK-A-M-Sr^^] =0, for allege J, 



where we replaced d acting on -J- by A, and d acting on ^ by j-L. The latter equation is the 
same as (1.31), due to equation (1.29). □ 

A more difficult proof of this proposition, in terms of Hamiltonian operators, is given in [O]. 



2 Compatible PVA structures and the Lenard scheme of integrability. 

2.1 Compatible PVA structures. In this section we discuss a method, known as the Lenard 
scheme of integrability, based on a hi- Hamiltonian pair (H, K), to construct an infinite sequence 
of local functionals Jh n G V/dV, n G Z+, which are pairwise in involution with respect to both 
Hamiltonian operators H(d) and K(d) (cf. Definition 1.25(d)). This scheme will produce an 
infinite hierarchy of integrable Hamiltonian equations 

(2-1) ^ = {K\u}h\ x=0 ( = {h n +i x u}K\ x=0 ^ , n£Z +1 

such that the local functionals Jh n , n G Z+, are integrals of motion of every evolution equation 
in the hierarchy. 

Definition 2.1. Several A-brackets { . \ . } n , n = 1, 2, . . . , N, on a differential algebra V are 
called compatible if any C-linear combination { • a • } = ^2n=o Cn i • A ■ }n of them makes it a PVA. 
If V is an algebra of differential functions, and H n (d), n = 1, 2, . . . , N, are the corresponding 
to the compatible A-brackets Hamiltonian operators, defined by (1.47), we say that they are 
compatible as well. A bi- Hamiltonian pair (iJ, K) is a pair of compatible Hamiltonian operators 
H(d),K(d). 

Example 2.2. (cf. Examples 1.18 and 1.19) Let R = C[u, u', u", . . .]. The A-brackets 

{u x u}i = (<9 + 2A)u, {u A u} 2 = A, {u x u} 3 = A 3 , 
are compatible. The corresponding Hamiltonian operators, defined by (1.47), are: 

H 1 (d)=u' + 2ud, H 2 (d) = d, H 3 (d) = d 3 . 
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Example 2.3. (cf. Example 1.20) Let g be a Lie algebra with a non-degenerate symmetric 
bilinear form (.|.), let {uj}j e / be an orthonormal basis of g and [uj,Uj] = J2k c ij u k- Let 

R = C[u[ n ^\i G I , n € Z}. The following A-brackets on R are compatible: 

{uixUj}' = ^2 c ij u k , {ui x Uj}" = XSij , {u ix Uj} k = cfj , kel. 
k 

The corresponding Hamiltonian operators, H',H" and H k , k £ I, are given by: 

H> J (d) = -^2c^u k , H? j (d) = 5 ij d, fl*(0) = -c*. 
kei 

Definition 2.4. Let V be an algebra of differential functions and let H(d) = [Hij^S)^ . gJ , 
K(d) = (Kijip)^ be any two differential operators on V®^. An (H, K)- sequence is a 
collection {F n } <„<Ar C V® £ such that 

(2.2) if (d)F n+1 = H (d)F n , < n < N - 1 . 
I£ N = oo, we say that {F n } n ^z + is an infinite (H,K)- sequence. 

Remark 2.5. Equation (2.2) for an infinite (H, if)-sequence can be rewritten using the gener- 
ating series F(z) = Ylne%+ F n z~ n ^ as follows: 

(2.3) (H(d) - zK{d))F{z) = -zK(d)F° . 
Assuming that 

(2.4) K(d)F° = 0, 

and taking the pairing of both sides of (2.3) with F(z), we get the equation 

(2.5) F(z) ■ (H{d) - zK(d))F(z) = 0. 

Note that for a skew-adjoint differential operator M(d), we have jF ■ M{d)F = for every 
F E V®^. Therefore, if both operators H(d) and K(d) are skew-adjoint, the LHS of the equation 

(2.5) is a total derivative, hence we can reduce its order. If, in addition, I = 1 and K(d) is a 
differential operator of order 1, this produces a recursive formula for the sequence F n , n > 1. 
Explicitly, for K(d) = d(k) + 2kd, where k £ V, equation (2.5) becomes 

1 n i n 

(2.6) d(kF°F n+1 ) = F n ~ m H{d)F m - - ^ d(fcF ra+1 - m F m ) . 

m=0 m=l 

Notice that the RHS is a total derivative (since H(d) is skew-adjoint), so that, given F° sat- 
isfying (2.4), the above equation defines recursively F n+1 up to adding a constant multiple of 
(fcF )- 1 . 

Note that in the special case in which F n = for some fh n £ V/dV, equation (2.2) can 
be written in terms of the A-brackets associated to the operators H{d) and K(d) (cf. (1.47)): 

{h n+ \ \ u}k\ x=0 = {h n \u}H\ x=0 - 
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Lemma 2.6. Suppose that the operators H(d), K(d), acting on V® e , are skew- adjoint. Then 
any (H, K) -sequence {-F n }o<n<Ar satisfies the orthogonality relations: 

(2.7) fF m ■ H(d)F n = jF m ■ K(d)F n = , < m, ra < iV . 

Proof. If m = n, (2.7) clearly holds, since both H (d) and K(d) are skew-adjoint operators. We 
may assume, without loss of generality, that m < n (hence n > 1), and we prove equation (2.7) 
by induction on n — m. By (2.2) and the inductive assumption we have 

jF m ■ K(d)F n = jF m ■ H(d)F n ~ 1 = , 

and similarly, since H(d) is skew-adjoint, 

jF m ■ H(d)F n = -jF n ■ H{d)F m = -fF n ■ K(d)F m+1 = . 

□ 

The key result for the Lenard scheme of integrability is contained in the following theorem, 
which will be proved in Section 4. 

Theorem 2.7. LetV be an algebra of differential functions in the variables {uj}j g /, and suppose 
that H(d) = (i?jj(9)). . j and K(d) = (ify(3)) i . j form a bi-Hamiltonian pair. Assume 
moreover that K{d) is non- degenerate, in the sense that K{d)M{d)K(d) = implies M(d) = 
for any skew- adjoint differential operator M{d) of the form Df(B) — D* F {d), F G V . Let 
{i ?n } n= o,--- ,v C V®^ (N > 1, and it can be infinite) be an (H, K)-sequence, and assume that 
F ° = ' Fl = 1&> f or two local Junctionals jh , Jhi G V/dV. Then all the elements of the 
sequence F n , n = 0, • • • ,N, are closed, i.e. (1.16) holds for F = F n . 

Proof. This theorem is a corollary of the results in Section 4, see Remark 4.26. □ 

Remark 2.8. If I = 1, the non-degeneracy condition for the differential operator K{d) in Theo- 
rem 2.7, i.e. that K{d)M{d)K(d) = implies M{d) = 0, holds automatically unless K(d) = 0. 
It is because the ring of differential operators has no zero divisors when I = 1. For general 
£ a sufficient condition of non-degenaracy of K(d) is non-vanishing of the determinant of its 
principal symbol. 

Theorem 2.7 provides a way to construct an infinite hierarchy of Hamiltonian equations, 
4~ = {Jh n , u}h i and the associated infinite sequence of integrals of motion Jh n , n G Z+. In 
order to do this, we need to solve two problems. First, given a bi-Hamiltonian pair H{d),K(d) 
acting on V®^, we need to find an infinite (H, if )-sequence {F n } n£ z + , such that F°, F 1 G Im 4-;. 
Second, given an element F G V®^ which is closed, we want to prove, that it is also "exact" , i.e. 
that F = || for some local functional Jh G V/dV, and we want to find an explicit formula for 
the local functional J h G V/dV. Then, if we have F n = ^ for all n G Z+, by Lemma 2.6, the 
corresponding local functionals Jh n are pairwise in involution with respect to the Lie brackets 
associated to both Hamiltonian operators H(d) and K{d): 

(2.8) {Jh m , Jh n } H = {Jh m ,Jh n } K = , for all m, n G Z + . 

The solution to the second problem is given by Propositions 1.9 and 1.11. We shall discuss, 
in Propositions 2.9 and 2.10, a solution to the first problem: we find sufficient conditions for 
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the existence of an infinite (H, FC)-sequence {F n }nez + > and we describe a practical method to 
construct inductively the whole sequence. 

Let V be an algebra of differential functions. Given a subset U C V® , let 

U x := {P G V e | JF • P = , for all F 6 [/} C V e , 

and similarly, given a subset C7 C V , let 

{/^ := {F G V® 1 | /F • P = , for all P € U} C V® £ . 

Proposition 2.9. Let H(d) = (Hij(d)) i - eV K(d) = {K%j {d)) i eI , be skew-adjoint operators 
on V® 1 , and let {P n }o<n<A r be a finite (H, K)- sequence satisfying the following condition: 

(2.9) (SpanciF^oKnKN) 1 C Im K{d) . 

Then it can be extended to an infinite (H, K) -sequence {F n } ne z + ■ 

In particular, letting U = Span c {F n | n G Z + } C V®^ and f7' = Span c {F n I n > 1} C V®^, 
we /iaue: 

(2.10) U'cU , H{d)U = K{d)U' , U x C Im K{d) . 

Proof. We prove, by induction on n > N + 1, that F" G V®^ exists such that K{d)F n = 
H(d)F n ~ l . Indeed, by the inductive assumption {F fc }o<fc<n-i is an (H, F)-sequence, therefore 
by Lemma 2.6 H(d)F n ~ 1 _L F m for every m G {0, ...,N}, and hence by assumption (2.9) 
H(d)F n ~ 1 G Im K(d), as we wanted. The second part of the proposition is obvious. □ 

Proposition 2.10. Let V be an algebra of differential functions in the variables {uj}j e / and 
let H(d), K(d), be skew-adjoint operators on V® e . Suppose that U,U' C V®^ are subspaces 
satisfying conditions (2.10). Then: 

(a) Ker K{d) C (H(d)U) ± , 

(b) for every F° G [H(d)U) ± there exists an infinite (H, K) -sequence {F n } n& z + , and this 
(H, K) -sequence is such that Span c {F n } n£ z + C (H(d)U) ± , 

(c) if F° G Ker K{&) and {F n } n& z + , {G n } n& z + are two infinite (H,K)- sequences, then they 
are "compatible" , in the sense that 

(2.11) jF m ■ H{d)G n = jF m ■ K{3)G n = , for all m, n G Z+ . 

Proof. Since K(d) is skew-adjoint, we clearly have Im K(d) C ( Ker K ((9)) . Hence, from 

H(d)U = K{d)U' C Im K{d) c ( Ker K(8))^, it follows that Ker K{d) C (H(d)U) ± . For part 

(b), we first observe that, if {F n }o<n<7V is any (H, FC)-sequence starting at F° G (H(d)U) ± , 

then F n G (H(d)U\ for every n G {0, . . . , N}. Indeed, let 1 < n < N and suppose by induction 

that F n ~ l G {H{d)U) L . If G G U, by assumption (2.10) we have that H(d)G = K(d)G x for 
some G\ G U' C U. Hence 

fF n ■ H(d)G = fF n ■ K(d)G\ = -fd ■ K{d)F n 

= -fGx ■ H{d)F n ^ = jF n ~ l ■ H(d)G\ = , 
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hence F n G (H{d)U~) , as claimed. Next, we need to prove the existence of an infinite (H,K)- 
sequence starting at F°. Suppose, by induction on N > 0, that there exists an (H, ET)-sequence 
{F n }o<„,<Ar. By the above observation, F N G [H(d)U) ± , or equivalently H(d)F N G U 1 , so 
that, by the third condition in (2.10), there exists F N+l G V® 1 such that H(d)F N = K(d)F N+1 , 
as we wanted. We are left to prove part (c). We will prove equation (2.11) by induction on 
m G Z+. For m = we have fF° ■ K(d)G n = -jG n ■ K(d)F° = 0, and JF° ■ H(d)G n = 
JF° ■ K(d)G n+1 = 0, for every n G Z + . We then let m > 1 and we use inductive assumption 
to get, for every n G Z+, 

fF m • K(d)G n = -fG n ■ K{d)F m = -jG n • H(d)F m " 1 = fF m ^ • H(d)G n = , 

and similarly 

jF m ■ H{d)G n = fF m ■ K{d)G n+1 = . 

□ 

Remark 2.11. If H(d) and K(d) are pseudodifferential operators on V , i.e. the corresponding 
A-brackets on V are "non-local", then the statements of Propositions 2.9 and 2.10 should be 
modified to make sure that the elements F n , n > 0, belong to the domains of the operators 
H(d) and K{9). In Proposition 2.9 we then need the further assumption 

(2.12) (Span c {H(d)F n } < n < N ) ± C BomH(d) . 
Similarly, in Proposition 2.10 we need the further assumption 

(2.13) (H{d)U) L C DomH(d) . 

By Proposition 2.9, a sufficient condition for the existence of an infinite (H, fC)-sequence 
{F n } ne z + , is the existence of a finite (H, K)-sequence {-F n }o<n<Af satisfying condition (2.9). By 
Proposition 2.10, if a finite (H, if)-sequence satisfying condition (2.9) exists, then any element 
F° G Ker K(&) can be extended to an infinite (H, if)-sequence {F n } n( zz + . Therefore, in some 
sense, the best place to look for the starting point of an infinite (H, _fT)-sequences is Ker K(d). 

Moreover, by part (c) of Proposition 2.10, the infinite (H, X)-sequences starting in Ker K{d) 
are especially nice since they are compatible with any other infinite (H, if)-sequence. Notice 
that the space Ker K(d) is related to the center of the Lie algebra V/dV with Lie bracket 
{•, -} K - Indeed, if F° = ^ belongs to Ker K(d), then the local functional fh G V/dV is 
central with respect to the Lie bracket {• , -}k, and the converse is almost true (cf. Corollary 
1.32). Proposition 2.10(c) says that, if F° G Ker K(d), then the corresponding infinite (H,K)- 
sequence {i ?n }nez + is "central" among all the other infinite (H, if)-sequences, in the sense of 
the compatibility condition (2.11). 

We conclude this section with the following result, which incorporates the results of Theorem 
2.7 and Propositions 1.9, 2.9 and 2.10. It will be used in the following Sections 2.2 and 2.3 to 
implement the Lenard scheme of integrability in some explicit examples. 

Corollary 2.12. Let V be a normal algebra of differential functions. Let H(d), K(d) be a 
bi-Hamiltonian pair acting on V® e , with K(d) non- degenerate. Let [ho, [hi G V/dV satisfy 
the following conditions: 
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(i) H(d) s -£=K(d) s -£, 

(ii) (Span c { s -t,*t}y ClmK(d). 

Then there exists an infinite sequence of local Junctionals Jh n , n G Z + , which are pairwise in 
involution with respect to the Lie brackets associated to both H{d) and K(d), i.e. (2.8) holds. 
Such hierarchy of compatible Hamiltonian functionals can be constructed, recursively, using 
Proposition 1.9 and the following equations: 

(2.14) H (d)^ = K(d) 6 -^ , neZ + . 

du ou 

Proof. By the conditions (i) and (ii), the elements {^fj^j} form a finite (H , i^)-sequence 
satisfying condition (2.9). Hence, by Proposition 2.9, they can be extended to an infinite 
(H, /^-sequence {F n }n&+ such that F° = ^ and F 1 = Since, by assumption, (H, K) 
form a bi-Hamiltonian pair and K is non-degenerate, we can apply Theorem 2.7 to deduce that 
every element F n , n G Z + , is closed. Hence we can apply Proposition 1.9 to conclude that F n 
is exact, namely there is h n G V such that F n = □ 

Proposition 2.13. LetV be a normal algebra of differential polynomials in one variable u. Let 
K(d) = d and let H(d) be a Hamiltonian operator forming a bi-Hamiltonian pair with d. Then 
there exists Jh± £ V/dV such that 

(2.15) H(d)l = d^-. 

ou 

Furthermore, we can extend fho = fu and Jh\ to an infinite sequence fh n ,n£ Z + , satisfying 
(2.14). 

Proof. First, notice that, since H(d) is a skew-adjoint operator, jH(d)l = 0, namely there 
exists F G V such that H(d)l = dF. Next, the condition that the Hamiltonian operators H{d) 
and K = d form a bi-Hamiltonian pair reads 

{u x H(jj)} k - { Ufl H(X)} K = {H(X) x+ ,u} K . 

Using equation (1.29) and letting A = we then get 

Z-/ Q u (n) VP ) g uin) 

Since H(0) = dF, using the commutation rule (1.2), we obtain, after simple algebraic manip- 
ulations, that fi(d + ii)(Df([i) — Dp(—n — <9)) = 0. This implies that Dp{d) is a self-adjoint 
differential operator, i.e. F G V is closed. By Proposition 1.9, F is a variational derivative. 
Thus, all the assumptions of Corollary 2.12 hold, which proves the proposition. □ 
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2.2 The KdV hierarchy. Let R = C[u, u', u", . . . ] with the following two compatible A- 
brackets (see Example 2.2): 

{u\u}h = {u\u}i + c{n A n} 3 = (d + 2X)u + cA 3 , 

(2.16) {u\u} K = {u x u} 2 = \, 

where c £ C, so that the corresponding Hamiltonian operators are: 

(2.17) H{d) =u' + 2nd + cd 3 , K{d) = d . 

We shall use the Lenard scheme discussed in the previous Section, to construct an infinite 
hierarchy of integrable Hamiltonian equations, which includes the classical KdV-equation 

(2.18) ^=3uu' + cu"'. 

By Proposition 2.13, we conclude that there exists a sequence of local functionals Jh n £ 
R/dR, n £ Z+, where ho = u, which are in involution with respect to both Lie brackets {• , -}h 
and {• , and such that the corresponding variational derivatives F n = , n £ Z + form 
an infinite (H, ET)-sequence (cf. Definition 2.4). 

In fact, we can compute the whole hierarchy of Hamiltonian functionals jh n , n £ Z + , 
iteratively as follows: assuming we computed Jh n , n £ Z + , we can find Jh n+ i by solving 
equation (2.14). Clearly F n = is defined by equation (2.14) uniquely up to an additive 
constant. To compute it explicitly we have to find a preimage of H(d)F n ~ 1 under the operator 
d, which is easy to do, inductively on the differential order of H{d^F n . Similarly, Jh n is 
defined uniquely up to adding arbitrary linear combinations of jl and J u, and to compute it 
explicitly we can use Proposition 1.11, since Ker (A + 1) = 0. For example, for the first few 
terms of the hierarchy we have: 

13 11 

F° = 1 , h = u; F 1 =u , hi = -u 2 ; F 2 = -v 2 + cu" , h 2 = -u 3 + -cuu" ; 

(2.19) F 3 = -u 3 + bcuu" + -c(u') 2 + c 2 u^ , /* 3 = V + -cu 2 u" + -cuu' 2 + -cW 4 ) . 
v; 2 2 v; ' 83 6 2 

By Corollary 2.12 the local functionals Jh n , n £ Z + , are in involution with respect to both 
Lie brackets associated to H(d) and K(d), i.e. (2.8) holds. Hence, if we let J h-\ = Jl, we have 
an increasing filtration of abelian subalgebras of R/dR with Lie bracket {• , -}k (or {• , •}#), 
f)" 1 C £)c C S)l C 9) 2 C C . . . , where $j™ = Span c {/ hk}^ = _i- Notice that the space ft™ does not 
depend on the choice of the local functionals Jh n , n £ Z + , solving equation (2.14). We also 
denote Sj c = [} n&+ ft™ = Span c {//i n }^ = _ 1 . 

The local functionals Jh n , n > are integrals of motion of the hierarchy of evolution 
equations (2.1). The th and the 1 st equations are gjp = and jfj- = §j, and they say that 
u(x, to, ti, . . . ) = n(x+ti, t2 ; • • • )> the 2 nd equation is the classical KdV equation (2.18) provided 
that c / (then c can be made equal 1 by rescaling), and the 3 rd equation is the simplest higher 
KdV equation: 

— = — u 2 u' + Wcu'u" + 5cuu"' + cV 5) , 
dt 3 2 

provided again that c ^ (it can be made equal to 1 by rescaling). 
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Remark 2.14. Equation (2.6) for k = 1/2 gives an explicit recursion formula for the sequence 
F° = 1, F , F 2 , . . . , which in this case reads 

n ^ n 

(2.20) F n+l = (cF n - m {d 2 F m ) - ^(dF n ~ m )(dF m ) + uF n " m F m ) - - ^ F n+1 " m F m , 

m=0 m=l 

up to adding an arbitrary constant. This gives an alternative way to compute the F n, s. 

If we put c = in the KdV hierarchy (2.19) we get the so called dispersionless KdV hierarchy. 
We denote in this case F and Jh n in place of F n and jh n . The recursion relation (2.14) for 
the elements ~F n = n E Z+, is 

so if we put, as before, P° = 1, it is immediate to find an explicit solution by induction: 

— n (2ra-l)H „ r - (2n-l)!!„„, 1 

(2.21) F = { - ' u n , Jh n = { {n + i ! )[ Ju n+1 ,nGZ + , 

where (2n — 1)!! := 1 • 3 • • • (2n — 1). The corresponding dispersionless KdV hierarchy is 

du (2n-l)!! 



dt n n\ 



-d{u n ) , n G Z_ 



As before, we denote Jh-i = J 1, and we have a uniquely defined increasing filtration of abelian 
subalgebras, with respect to the Lie bracket {• , -}k (or {• , -}h with c = 0): Sj^ 1 C C -fjg C 
• • • C £o, where f$ = Span c {/7i fc }^ = _ 1 and S) = IXL-i^o- 

Finally, as c — > oo the KdV hierarchy turns into the linear KdV hierarchy, corresponding 
to the pair H(d) = d 3 , K(d) = d. We denote in this case FJ 1 and Jh n in place of F n and fh n . 
The Lenard recursion relation (2.2) becomes clF n+1 = d 3 F^, n > 1, and a solution of such 
recursion relation, such that all F_ n are in the image of 4r, is obtained starting with F} 



u. 



In this case, F n+1 is defined by F_ n uniquely up to adding an arbitrary constant, while fh n+ i 
is defined uniquely up to adding an arbitrary linear combination of J 1 and J u. An explicit 
solution is given by 

(2.22) = , fh n+1 = l -Juu^ = LiL/( u (»)) 2 , n G Z + , 

giving the linear KdV hierarchy 

Ctt n 

By Theorem 2.7, the local functionals //i n , n > 1, are in involution with respect to both the 
Lie brackets {• , -}h and {• , If we let Jh_i = jl and jh = ju, we thus get, as before, a 
uniquely defined increasing filtration of abelian subalgebras, with respect to both {• , -}k and 
{• , -}h- fi^ C Sj ^ C Sj 1 ^ C • • • C £oo, where = Span c {//i fc }^ = _ 1 and ijoo = U£L-i#£>- 
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Theorem 2.15. The spaces $j c = Span c {fh n } ( ^ = _ 1 , c / 0, % = ^P an c{I^ l n} C ^ = _ l , and 
fioc = Span c {J "h n )°^__ l , defined by (2.19), (2.21) and (2.22) respectively, are maximal abelian 
subalgebras of the Lie algebra R/dR with the bracket 

(2-23) { JfJ g} = J*L d V (of. (1.60)). 

Moreover, the local Junctionals jh n , n > — 1 (resp. jh n , n > —1, and fh n , n > —1), form a 
maximal linearly independent set of integrals of motion of the KdV (resp. dispersionless KdV, 
or linear KdV) hierarchy. 

Proof. We consider first the case of the dispersionless KdV hierarchy Jh n , n > 1. Notice that, 
by (2.21), we have S)o = C[u] C R/dR. Let J" / G R/dR be a local functional commuting with 
all elements Jh n , n > 1. We have 

f u n d 6 4- = -^—{ ff, h n+1 } = , for every n G Z+ . 
J ou n + 1 

Hence by Lemma 1.28 we get that <9g G 0C[u], namely (since Ker (d) = C) | S C[n]. On the 
other hand C[u] = j^(C[u]), so that we can find a(u) G C[u] such that |^ = . Equivalently, 
/ - a(u) G Ker (^) = C + d# (see Proposition 1.5), so that f edR + C[u], as we wanted. 

Next, we consider the linear KdV hierarchy jh n ,n > —1. Let // G R/dR be a local 
functional commuting with all elements JT^, re > — 1, namely, since —5^- = tt^ 2n \ such that 

every n G Z + . 

on 

By Lemma 1.30 this is possible if and only if G d Span c {u^ 2n ^ , n G Z+}, which is equivalent 
to saying that ff G Span c {/l,/u, jW 2n )}, as we wanted. 

We finally prove the general case, for arbitrary value of the central charge c G C, by reducing 
it to the case c = 0. According to the decomposition (1.88), the elements F n , n G Z+, defined 
by the recurrence formula (2.20), can be written, for arbitrary c G C, in the form 

(2.24) F n = F n + G n , 

where F G C[u] is given by (2.21) and G n G R 9 . Likewise, the local functionals fh n £ R/dR 
in (2.19) admit a decomposition 

(2.25) h n = h n + r n , 

where h n G C[u], is given by (2.21), and r n G R d . Indeed, let tto : R -» C[u] be the projection 
defined by the decomposition (1.88). We need to prove that 7ro(F n ) = F and no(h n ) = h n 
for every n G Z+. If we take the image under no of both sides of equation (2.20), we get the 
following recurrence relation for the elements iro(F n ) G C[u]: 

n 1 n 

MF n+1 ) = ^nvro(F"- m )7ro(F m )--^7ro(F" +1 - m )7ro(F m ). 

m=0 m=l 
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Since in the above equation the central charge c does not appear anywhere, its solution 7TQ(F n ) £ 
C[u] is the same for every value of c, in particular it is the same as the solution for c = 0, i.e. 
7To(-F n ) = F . Furthermore, by Proposition 1.11, a representative for the local functional 
fh n £ R/dR can be obtained by taking h n = A~ 1 (uF n ), where by A" 1 we mean the inverse of 
A, restricted to the space of differential polynomials with zero constant term. Therefore such 
h n clearly has a decomposition as in (2.25) since A-^ttF 1 ) = h n G C[u], while A~ l (uG) G R a 
for every G G R 9 . 

The basic idea to reduce to the case c = is to consider the change of variable x ^ x/e 
in u = u(x) and take the limit e — > 0. In order to formalize this idea, we define the map 
ir £ : C[A] ® R -> C[e, A] ® R given by 

ir £ (F(\;u,u',u",...)) = F{e\;u,eu ,e 2 u" , ...) . 

Clearly, we have induced maps ir £ both on the space of differential polynomials, ir £ : R — > 
C[s] ® R, defined by the natural embedding R C C[A] <8> R, and on the space of local functionals 
vr e : R/dR -► C[e] ® ~ (C[e] ® i?)/d(C[e] ® R), defined by vr e (//) := firj. The 

decomposition (2.25) can then be restated by saying that 

(2.26) n £ (h n ) = h n + er £tU , 

where h n is given by (2.21) and r £tTl G C[e] <8> R- 

It is clear by the definition of ir £ that 7r £ = g n 7r £ , so that Jjvr e = vr e ^-. In particular, 
recalling the definition (2.23) of the Lie bracket on R/dR, we have 

(2-27) {/*./, jTTrf} = e-^UfJg}, 

for every f,g £ R. Let then J/ G R/dR, for some / G i?, be a local functional commuting with 
all /ft „, ra G Z_|_, and assume, by contradiction, that ff$ ^P an c{/^ i ra} n= _ 1 - Notice that, by 
the decomposition (2.25) and by equation (2.21), we can subtract from / an appropriate linear 
combination of the elements h n , so that, without loss of generality, / G R 9 and J f 7^ 0. In 
other words, we have 

(2.28) n £ f = e k (J + er £ ). 

where k > 1 and / G is such that f / 7^ 0. By equation (2.27) and by the assumption on 
J* /, we have {J 7r E /, J* vr e /i n } = for all n G Z+, namely 

(2.29) {vr £ / A vr £ /i n }| A=0 G C[e] for all n G Z+ . 

On the other hand, by the decompositions (2.26) and (2.28), the minimal power of e appearing 
in the LHS of (2.29) is e k , and looking at the coefficient of e k we get that 

{/A^-n}| A=0 G dR, for all ra G Z + , 

namely {//i/^n} = for all n G Z+. By the maximality of i^o ; it follows that J"/ G f)0) 
namely / G C[u] + OR, contradicting the assumption that / G R d \dR. 

Finally, all local functionals Jh n are linearly independent. Indeed, since = F n , it suffices 
to show that the elements F n , n G Z + are linearly independent in R = C[u^ ra ^ | n G Z+]. For 
c = (respectively c = 00) this is clear from the explicit expression of F in (2.21) (resp. 
F n in (2.22)). For c G C\{0}, it follows by the decomposition (1.88) and by the fact that 
ttq(F u ) = F" G C[u], ra G Z+. □ 
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In view of Remark 1.26, the proof of Theorem 2.15 shows that the KdV hierarchy is inte- 
grable. 

The following result shows that the only hierarchy that can be obtained using the bi- 
Hamiltonian pair (H = d 3 ,K = d) is fh n defined in (2.22), for any choice of Jho, fh\. We 
believe that the same should be true for all values of c. 

Proposition 2.16. All the solutions Jf,Jg€ R/dR of the following equation 

(2.30) 5 -f = d 25 f. 

ou ou 

are obtained by taking linear combinations of the pairs fh n +i, fh n , n > 0, defined in (2.22). 

Proof. Let N = maxjord (/), ord (g)+l} : so that / = f(u, u',. . . , u'^'), g = g(u, u', . . . , u^ N ~^). 
By induction, we only need to show that 

(2.31) / = 7 ( u W) 2 + /l) 9 = _ 7 ( U (^D) 2 +51 moddi?, 

for some 7 G C and f\,gi G R such that ord(/i) < N — 1, ord (gi) < N — 2. We start by 
observing that, if (ft S R has ord {4>) < N — 1, then there exists ip £ R with ord {ip) < N — 1 such 
that (fiu( N ^ = ip mod dR. Indeed, in general we have 4> = Sfc>o 4'k( u ^ N ~ 1 " > ) k with ord^^) < 
N — 2, and we can then take ip = — ^2 k>0 -Ar[(d4>k) (u^ N ~^) +X . Therefore, since / and g are 
defined modulo dR, we can assume, without loss of generality, that /, as a polynomial in u^ N \ 
and g, as a polynomial in u^ N ~^, do not have any linear term. Under this assumption, claim 

(2.31) reduces to saying that 

d 2 f d 2 q 

(2.32) — = G C. 

By assumption ord(/) < N, hence we can use the commutation rule (1.2) and the definition 
(0.7) of the variational derivative to get 

d o~f VV -nn „n 0/ / iv!V &*f 



Q U (2N) 5u I Q U (2N) Q u (n) *■ > g {N) 2 ' 
n=0 

g ^ = fv-p* 9 d n — = (-l)^iV^L 

n=0 



and similarly, since ord (g) < N — 1, we have 



7V-1 



2 



n=0 

a 2 ^ = Vf-ir - d n+2 -^- = (-D N - 1 (N + i)d 929 



d u (2N-l) Su > Q U (2N-1) Q u (n) V ' V ■ , fl (JV-1) 

n=0 

The first identity in (2.33) and the first identity in (2.34), combined with assumption (2.30), 
give 

, s d 2 f 8 2 g 
2.35 J -^ = , 
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while the second identities in (2.33) and (2.34) give 

d 2 f d 2 q 
(2.36) N8— A, = -(N + l)8 y 9 ■ 

Clearly, equations (2.35) and (2.36) imply (2.32), thus completing the proof. □ 



2.3 The Harry Dym (HD) hierarchy. We can generalize the arguments of Section 2.2 to 
the case in which we take {\}h and { a }k to be two arbitrary linearly independent linear com- 
binations of the A-brackets { \ } n , n = 1,2, 3, in Example 2.2. Namely, in terms of Hamiltonian 
operators, we take 

(2.37) H{d) = a x (u 1 + 2nd) + a 2 d + a 3 d 3 , K(d) = fa (u + 2nd) + f3 2 d + f3 3 d 3 , 

acting on V, a certain algebra of differential functions extending R = C[u^ \ n G which 
will be specified later. 

Clearly, if we rescale H(d) or K(d) by a constant factor 7 G C\{0}, we still get a bi- 
Hamiltonian pair and the theory of integrability does not change: any (H, K) sequence {F n }o< n <N 
(cf. Definition 2.4) can be rescaled to the ('yH, if)-sequence {7 n i ?n }o<n<Af ; or to the {H,^K)~ 
sequence {7 _n i ?n }o<n<Af- The corresponding local functionals Jh n , n G Z+, satisfying F n = 
can be rescaled accordingly. Furthermore, if (H,K) is a bi-Hamiltonian pair, so is 
(H + aK, K) for any a G C. It is immediate to check that, if {^"jo^n^Af is an (H, iT)-sequence, 
then 

Fa =Y^( r f)a k F n - k , 0<n<N, 
k=o ^ ' 

defines an (H + aK, i^)-sequence. In particular, if {fh n } n ^z + is an infinite hierarchy of local 
functionals satisfying equation (2.14), hence in involution with respect to both Lie brackets 
{• , -}h and {• , -}k (or any linear combination of them), then J/i Q , n = Ylk=o (fc) a ^ fhn-k, n> G 
Z+, defines an infinite hierarchy of local functionals satisfying (2.14) with H replaced by H + 
aK, hence still in involution with respect to both the Lie brackets {• , -}h and {• , -}k- The 
corresponding increasing filtration of abelian subalgebras, f)° C f) 1 C f) 2 C ...V/dV, given 
by $) N = Span c {//i Qjn } 0<n<Ar , is independent of a. In conclusion, from the point of view 
of integrability and of the Lenard scheme as discussed in Section 2.1, nothing changes if we 
replace the bi-Hamiltonian pair (H,K) by any pair of type (aK + j3H^K), with a, (3, 7 G C 
and /3,7 7^ 0. Namely the integrability depends only on the flag of Hamiltonian operators 
CK c CH + CK. 

We next observe that, in order to apply successfully the Lenard scheme of integrability, we 
must have f3 3 = 0. Indeed, suppose by contradiction that (3% 7^ 0. By the above observations 
we can assume, after replacing H with H — ^K, that 03 = 0. In order to apply the Lenard 
scheme, we need to find a sequence F n , n G Z+, solving equations (2.2). By looking at the 
differential order of the various terms in equation (2.2), unless F n G C, we get that ord(F ra ) + 
3 = ord (F n ~ l ) + 1, namely the differential order of F n decreases by 2 at each step, which is 
impossible since differential orders are non-negative integers, unless F n G C for n » 0, i.e. the 
Lenard scheme fails. 
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Let then K(d) = fa(v! + 2nd) + fad. If fa = 0, we may assume that K(d) = d and 
H(d) = ax(u' + 2nd) + a 3 d 3 , which gives the KdV, the dispersionless KdV or the linear KdV 
hierarchies, if 0103 7^ 0, ot\ 7^ and 03 = 0, or a\ = and 03/0 respectively. These cases 
were studied in detail Section 2.2. In this section we consider the remaining case, when /3i 7^ 0. 
Then we may assume that H (d) = ad + (3d 3 , K(d) = v! + 2ud + jd, for a, j3, 7 6 C, and a, /? 
not both zero. 

As a last observation, we note that the case 7 ^ can be reduced to the case 7 = with 
the change of variable 

(2.38) u ^ u - I , vF> ^ , n > 1 , 

in the polynomial algebra = C[u^ \n G or in an algebra of differential functions V 

extending R. Indeed, it is immediate to check that (2.38) defines an algebra automorphism 
(p of R commuting with the action of d, hence the operator H{d) commutes with 99, while 
pKid)^- 1 = u' + 2nd. 

We thus arrive at the following one-parameter family (parametrized by P 1 ) of bi-Hamiltonian 
pairs: 

(2.39) H{d) = ad + (3d 3 , K(d) = v! + 2nd , 

where a, (3 S C are not both zero. We want to apply the Lenard scheme of integrability, namely 
Corollary 2.12, to construct an infinite hierarchy of local functionals Jh n , n £ Z + , which are 
pairwise in involution with respect to the Lie brackets {• , -}h and {• , -}k (thus defining an 
infinite hierarchy of integrable Hamiltonian equations (2.1)). 

As we pointed out in Section 2.1, the best place to look for the starting point F° of an 
infinite (H, iT)-sequence {i ?n } n ez+ is the kernel of the operator K(d). Namely we want to find 
solutions of the equation 

(2.40) K(d)F° = (u 1 + 2ud)F° = . 

If we write K(8) = 2u 1 / 2 d o -u 1 / 2 , we immediately see that, up to a constant factor, the only 
solution of equation (2.40) is F° = 1/ y/u. Therefore, if we want to have a non-zero kernel of the 
operator K(d), we are forced to enlarge the algebra R = C[u^ \ n £ Z + ] to a bigger algebra of 
differential functions V, which includes 1/^/u. The smallest such algebra (which is invariant by 
the action of d) is 

(2.41) V = C[u ± 5,M>",...]Di?. 

The degree evolutionary vector field A is diagonalizable on the algebra V with eigenvalues 
in |Z. If we let U = ua Cfit^ 1 , u' , u", . . . } C V, the operator A is diagonalizable on U with 
eigenvalues in 1>+\, so that both A and A+l are invertible. We also have K(d)~ l (H(d)U) C U. 
We then let 

(2.42) F° = , 

and we look for F 1 G V such that K(d)F 1 = H(d)F° . By the definition (2.39) of the operators 
H(d) and K(d), the equation for F 1 is 

2u^d{u^F l ) = (ad + (3d 3 )u~^ , 
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which, by the identities u zdu 2 = ^du 1 and u 2<9 3 u 2 = 2d(u id 2 u 4), gives, up to 
adding elements of Cu~2 5 

(2.43) F 1 = jau'^ + /?u~t<9 2 u~3 . 

Clearly F°, F 1 G U, and it is not hard to show that both F° and F 1 are in the image of the 
variational derivative j^. In fact we can use Proposition 1.11 to compute preimages ho, hi EU 
explicitly. Since F° (respectively F 1 ), is homogeneous of degree — \ (resp. — |), we have 
F n = ^f, n = 0, 1, where 

(2.44) h = 2u^ , hi = -^au~^ - 2f3u~^d 2 u~^ . 

We finally observe that the orthogonality condition (ii) of Corollary 2.12 holds. Indeed we have 
(Span c {F°, F 1 }) ± C (Cu - *)- 1 = u^dV = Im K{d) . 

We can therefore apply Corollary 2.12 to conclude that there exists an infinite sequence of local 
functionals Jh n , n G Z + , with h n in some normal extension of V, which are in involution with 
respect to both Lie brackets {• , and {• , -}k, and such that the corresponding variational 
derivatives F n = n G Z + form an (H, if)-sequence. In fact, we can compute the whole 
hierarchy of Hamiltonian functionals jh n , n G Z + , iteratively by solving equation (2.14). Since 
K \d)~ l (H (d)U) C U, we can choose all elements h n in U C V. 

In the special case a = 1, (3 = 0, we have a closed formula for the whole sequence F n , n G Z + , 
and the corresponding local functionals Jh n , n G Z + , which can be easily proved by induction: 

(2.45) F n = %^Er»-* , K = - (2n " 3) -,> -^ , n G Z + . 
v ; 2 n (2n)!! ' 2 n ~ 1 (2n)!! + 

For arbitrary a, (3 G C, the first two terms of the sequence are given by (2.42), (2.43) and (2.44), 
but already the explicit computation of the next term, namely F 2 , and the corresponding local 
functional J*/i2i is quite challenging if we try to use the equation K(d)F 2 = H(d)F 1 . 

Remark 2.17. Equation (2.6) for k = u gives an explicit recursion formula for the sequence 
F° = u"2, F , F 2 , . . . , which in this case reads 

1 n 1 1 

F n+l = -u~^ {-aF n ~ m F m + (3F n - m (d 2 F m ) - -[3{dF n ~ m ) (dF n 

m=0 



1 H 

(2.46) vh Y F n+l-m F m ^ 



m=l 

up to adding a constant multiple of u~?. It is clear from this recursive formula that F n has 
degree — 2w ^ 1 for every n G Z+, so that, by Proposition 1.11, we have F n = where 

(2.47) h n = -—^—uF n , n G Z+ . 

An — 1 
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We can use equations (2.46) and (2.47) to compute F 2 and h%. Recalling that F° and F 1 are 
given by (2.42) and (2.43) respectively, we have, after some lengthy algebraic manipulations, 

(2.48) F 2 = ^-a 2 u-i + ^-a(iu-id 2 u-^ + \l3 2 u-ld i u-^ , 

32 12 6 

7 1 2 — ^ " n — o2 — — 1 n 2 _ 1 4 _ 3 

no = a u 2 ap« 4# z u 4 tr?/ 4cTy, 4 . 

16 18 9 

The local functionals Jh n , n > 0, are integrals of motion of the integrable hierarchy of evolution 
equations, called the Harry Dym hierarchy: 

(2.49) ^ = ( a d + [3d 3 )F n , n G Z+ . 

The first three equations of the hierarchy are then obtained by (2.42), (2.43) and (2.48) respec- 
tively. The first equation for a = 0, (5 = 1 is the classical Harry Dym equation: 

— = a (u 2 . 
dt y ' 

Remark 2.18. As in the case of the KdV hierarchy, we point out that the abelian subalgebras 
£>»,/3 = {fh-n, n G Z + } C V/dV are infinite-dimensional for every a,/3 G C. Indeed, since 
= F n , it suffices to show that F n , n G Z+, are linearly independent in V. Since the 
operators H(d) and if(<9) in (2.39) have degree and 1 respectively, then, by the recursive 
formula (2.2), F n has degree — n — ^. Therefore, to prove that they are linearly independent it 
suffices to show that they are non-zero. It is clear from the definition (2.39) and some obvious 
differential order consideration, that Ker H{d) = C. In particular F n is in Ker H(d) if and 
only if F n = 0. But then we use equation (2.2) and induction to deduce that, if F n ~ l ^ 0, then 
F n_1 ^ Ker H(d), so that F n ^ 0. Hence, in view of Remark 1.26, the Harry Dim hierarchy is 
integrable. 

2.4 The coupled non-linear wave (CNW) system. Let V = C[4 n \ uW | n G Z+] be the 

algebra of differential polynomials in the indeterminates (u,v). Then the formulas 

(2.50) {u x u} H = (8 + 2A)n + cA 3 , {v x v} H = , 

-juA""}// = (0 + A)u , {v x u} H = Xv , 

and 

(2.51) {u\u} K = {v\v} K = A , {u x v} K = {v\u} K = , 

define a compatible pair of A-brackets on V for any c G C. Indeed, formulas (2.50) define 
a structure of a Lie conformal algebra on C[(9]u © C[<9]t> © C, corresponding to the formal 
distribution Lie algebra, which is the semi-direct sum of the Virasoro algebra and its abelian 
ideal, defined by the action of the former on Laurent polynomials C[t, t -1 ]. Moreover, adding 
a constant multiple of A in the first line of (2.50) corresponds to adding a trivial 2-cocycle [K]. 
The claim now follows from Theorem 1.15(c) or Example 1.22. 
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The Hamiltonian operators (on V® 2 ) corresponding to the A-brackets {• \ and {• \ -}k 



arc 



,„ KO s rj,^ ( cd 3 + 2nd + v! vd \ / 9 

(2 - 52) ^ ) = ( ttfW o J ' K(a) = Vo 



5feo \ / n \ / 



\ „ 1 ^ / 1 



Let ho = v, h\ = u £ V. We have ^SJ = ^ij'^' = 

and H(d)F° = K(d)F 1 = ( jj J . Moreover, obviously (CF° © CF 1 )^ = Im(d). Hence 

conditions (i) and (ii) of Corollary 2.12 hold, and the sequence J ho, J hi can be extended to an 
infinite sequence jh n G V/9V, n G Z + , satisfying (2.14). Thus, we get an infinite hierarchy of 
Hamiltonian evolution equations (n G Z+): 

d f n \ / §hu \ 

' = H(d) t n = K{d) 



dt \ v 




for which all Jh n are integrals of motion for both brackets {• , -}h and {• , -}k on V/dV, defined 
as in (1.67). One easily computes the whole hieararchy of Hamiltonian functionals Jh n and 
their variational derivatives F n G V® 2 by induction on n, starting with the given J ho and J hi, 
as we did for KdV and HD in Sections 2.2 and 2.3. For example, we have: 



F°=(°) , fto = »; F' = ( J) , >>!=«; 

f! = (") ' "H<» 2 +» 2 ); 

F d = I 2 2 I , /i 3 = -(cuu + u d + uv ) , 

V uv J 2 

and the first four equations of the hierarchy are: 

±( u )=0 ±( M = —( u ) = ( u ' 

dt \ v J ' dti \ v J ' dt 2 \ v J \ v' 

d ( u \ _ ( cu "' + 3uu' + vv' 
dt 3 \ v J \ d(uv) 

The last equation is called the CNW system [I], [D]. It is easy to see by induction that the 
first coordinate of F n G V® 2 has total degree n — 1 if n > 1. It follows that the F n , and hence 
the Jh n , are linearly independent for all n G Z_|_. Thus, in view of Remark 1.26, the CNW 
hierarchy is integrable. 



2.5 The CNW of HD type hierarchy. Recall that the three compatible Hamiltonian 
operators H n (d), n = 1,2, 3, from Example 2.2 were used to construct, in Section 2.2, the KdV 
hierarchy, and, in Section 2.3, the HD hierarchy. Here we describe an analogous picture for an 
algebra of differential functions V in two variables u and v. 
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We have three compatible A-brackets in two variables, {• \ •}„, n = 1, 2, 3, defined as follows: 

{uxujx = (d + 2X)u , {u A u}i = 0, {n A i;}i = (5 + X)v , {u A ^}i=Af, 
{ti A M} 2 = {v\v} 2 = A , {u A v} 2 = {v\u} 2 = , 

{""A^h = A 3 » { U A«}3 = {U\V}3 = {V\U} 3 = . 

The corresponding Hamiltonian operators H n (d) : V® 2 — > V 2 , n = 1,2,3, are given by the 
following differential operators: 

">< 8 > - ( u '/:: a a i ) ■ *«> = ( o s ) ■ *w = ( o 3 o 

It is clear from the discussion at the beginning of Section 2.4 that the above operators are 
all compatible, namely any linear combination of them is again a Hamiltonian operator. The 
CNW system is constructed using the bi-Hamilonian pair H = Hi + cH^, K = H2, and it can 
be viewed as the analogue, in two variables, of the KdV hierarchy. It is therefore natural to look 
for a two- variable analogue of the HD hierarchy, which should be associated to the 1-parameter 
family (parametrized by P 1 ) of bi-Hamiltonian pairs H = aH 2 + /3i?3, K = H\, namely 

, , / ad + /3d 3 \ fa . ( u' + 2ud vd 

(2 - 53) H ( d ) = { ad ) ' K{d)= { v' + vd 

with a, (3 £ C. Note that, using arguments similar to those at the beginning Section 2.3, any 
other bi-Hamiltonian pair (H, K) obtained with linear combinations of the operators Hi,H 2 , H3 
can be reduced, from the point of view of finding integrable systems, either to (2.52) or to (2.53). 

We next apply the Lenard scheme of integrability to construct an infinite hierarchy of 
integrable Hamiltonian equations associated to the Hamiltonian operators H(d) and K{d) in 

(2.53) , which we can call the "CNW of HD type" hierarchy. 

First, by Remark 2.8, the operator K(d) is non-degenerate. According to Proposition 2.10 
and the subsequent discussion, the best place to look for the starting point of an integrable 
hierarchy is the kernel of K{d). All the solutions of K(d)F = are linear combinations of the 
following two elements: 

(2.54) F ° = (?) ' Fl= (X) GVS2 - 

To make sense of them, we choose an algebra of differential functions V which contains ^, for 
example we can take 

(2.55) V = C[«,u ±1 ,«',i/,« /, y,...]- 
Note that F° and F 1 are the variational derivatives, respectively, of 

(2.56) Jho = Jv , fhi = J^€V/8V. 

Moreover, we clearly have H(d)F° = 0, so that condition (i) of Corollary 2.12 holds trivially. 
Next we want to check condition (ii) of Corollary 2.12, namely 

(2.57) (Span c {F° , F 1 }) 1 - C Im K(d) . 
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Let P = ( M G (Span^P^P 1 }^. Since jP • P° = 0, we have that q = d(vf) for some 
/ 6 V. Let r = ± (p - it'/ - 2u0/) , so that 

/ u'f + 2ua/ + w \ 
" I a(«/) y • 

In order to prove that P G Im we only have to check that r G <9V. Since P _L P , we have 

= fP • P 1 = f-(u'a + 2uda + w) - ^diva) = [ (r + d(-i 

V V 2 J \ V 



so that r G <9V, as we wanted. This proves (2.57). 

According to the above observations all the assumption of Corollary 2.12 hold. Hence there 
exists an infinite sequence of local functionals fh n E V/dV, n G Z+, with h n in some normal 
extension of V, which are in involution with respect to both Lie brackets {• , and {• , 
and such that the corresponding variational derivatives F n G V® 2 , n G Z+, form an (H,K)- 
sequence. In fact, we can compute the whole sequence of Hamiltonian functionals fh n , n G Z+, 
iteratively by solving equation (2.14). The corresponding hierarchy of evolution equations is 
given by (2.1), namely 



(2.58) 



& = {ad + ^) & -t 



Is the above hierarchy integrable? According to Definition 1.27, in order to establish inte- 
grability, it suffices to prove that the elements K(d)F n , n G Z + span an infinite-dimensional 
subspace of V 2 , and since dimKer K{d) = 2 < oo, it suffices to prove that the P n 's are linearly 
independent. 

As for the HD hierarchy, we observe that H(d) and K(d) in (2.53) are homogenous of degree 
and 1 respectively. Hence, by the recursive formula (2.14), we can choose F n G V® 2 , n G Z+, 
homogeneous of total degree — n. Two claims follow from this observation. First, we can choose 
the Jh n in V, due to Proposition 1.11. Second, to prove linear independence, we only have to 
check that the P n 's are not zero. 

We consider separately the two cases a = and a ^ 0. For a = we can compute explicitly 

the whole sequence: P° and P 1 are as in (2.54), P 2 = [ 3 , v /\2 v „ J , and P n = for all 

n > 3. In particular the P n 's span a finite-dimensional space and the Lenard scheme fails. 

If a we may rescale so that a is replaced by 1 and we let c = — . We want to prove 
that, in this case, F n ^ for every n > 0, namely integrability holds. Notice first that 
Ker (d) = Ker (d + c<9 3 ) = C. In particular, since each component of F n has degree —n, if, for 
n > 0, P™ 0, then necessarily (H(d)F n )i ^ 0, i = 1, 2. Moreover, we have by the expression 
(2.53) of K(d), that, for n G Z + , 



(2.59) (K(d)F n ) 2 = d(vF? 
while, if Pf = 0, then 

(2.60) (K(d)F n ) 1 = vdF? 
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Suppose that F£ 7^ for some n > 1. It follows that {H{d)F n )2 7^ and hence, by the recursive 
formula (2.2) and equation (2.59), that F"' +1 / 0. On the other hand, if Fg = and F? / 
for some n > 1, namely {H{d)F n )2 = and (H(d)F n )i ^ 0, we have, by the recursive formula 
(2.2) and equations (2.59) and (2.60), that = and F 2 n+1 + 0. In conclusion, F n + for 

every n € Z+, thus, in view of Remark 1.26, proving integrability of the hierarchy (2.58). 

We can compute explicitly the first few terms of the hierarchy. The first two Hamiltonian 
functionals [ho, Jh± are given by (2.56), the corresponding variational derivatives F°, F 1 are 
in (2.54), the 0-th evolution equation is trivial and the 1-st is 

_d_( u \ _ ( (d + cd s )(l) 

dti \v J \ -e(£) 

We call this equation the CNW system of HD type. We next compute F 2 by solving the 
equation H(d)F l = K{d)F 2 . The result is 



F 



/ - u 

3« 2 _H I IAv') 2 _ r v" 



which is the variational derivative of the Hamiltonian functional 

In 2 11 1 {v') 2 



f h * = -o^3+o- + 



-c- 



2v 3 2v 2 v 3 
The corresponding evolution equation is 



3 The variational complex. 

3.1 g-complex and its reduction. Given a Lie algebra q, we can extend it to a Lie super- 
algebra q containing g as even subalgebra as follows: 

(3.1) = 0[£] xO 5 = e £©Cc\, 

where £ is an odd indeterminate such that £ 2 = 0. A Q-complex is a Z + -graded vector space 
= ® ng 2 with a representation of on tt : q — > End 17, such that, for a E g, the 
endomorphisms 7r(a), vr(a£), 7r(<9g) have degrees 0, —1 and +1, respectively. 

Given an element d E 0, we can consider its centralizer g 9 = {a € | [5, a] = 0} C 0, and 
the corresponding superalgebra 9 = 9 © a £ © C<9g C 0. Notice that q 9 = q 9 . Let O be a 
0-complex. Clearly Tr(d)fl = © ragZ+ 7r(<9)f2 n C (l is a graded submodule with respect to the 
action of g 9 . We can thus take the corresponding quotient 9 -module Q = ©„ eZ+ where 
$7" = f2 n /ir(d)VL n , which is called the reduced (by 5) 5 -complex. 



56 



3.2 The basic and the reduced de Rham complexes f2 and over V. As in Section 
1, we let R = C[uf^ \ i G /, n G Z+] be the C-algebra of differential polynomials in the variables 
Ui, i G {1, . . . ,£} = I, and we let V be an algebra of differential functions extending R. Again, 
as in Section 1, we denote by q the Lie algebra of all vector fields of the form (1.10). 

The basic de Rham complex 0, = f2(V) is defined as the free commutative superalgebra over 

in) 

V with odd generators 5u\ , i G I,n G Z + . In other words £1 consists of finite sums of the form 

(3.2) £ ^ a • • • a *t4r } . cr sv . 

and it has a (super)commutative product given by the wedge product A. We have a natural 
Z + -grading O = © fcgZ+ defined by letting elements in V have degree 0, while the generators 

5uf^ have degree 1. The space Q k is a free module over V with basis given by the elements 

Su^™ 1 ^ A • • • A Su\ mk \ with (mi,ii) > • • • > (rnk,ik) (with respect to the lexicographic order). 

In particular £7° = V and ft 1 = n< =% + VSu^ . Notice that there is a V-linear pairing 

fi 1 x g — > V defined on generators by (fa -"^ , ^ ) = SijS mjn , and extended to fl 1 x g by 

du ] 

V-bilinearity. 

Furthermore, we want to define a representation tt : g — > Der f2 of the Lie algebra g on the 
space f2 given by derivations of the wedge product in Q, which makes into a 0-complex. We 
let 5 be an odd derivation of degree 1 of Q, such that Sf = J2iei nez+ & u i f° r f ^ V> an d 



<5(<5u- ) = 0. It is immediate to check that 5 2 = and that, for Cj G £l k as in (3.2), we have 

(3.3) 5(0)= £ £ A^fef A ^A...A^. 

mi,- ,mfcGZ-|_ 

For X G fl we define the contraction operator lx ■ f2 — > fi, as an odd derivation of Q of degree 
-1, such that i x (f) = for / G V, and t x (Su\ n) ) = X(u\ n) ). If X G is as in (1.10) and Q G U k 
is as in (3.2), we have 

(3.4) = x; D-rv^-rv^feir'A^A^'. 

ii,— r=l 
mi,-" ,m.fe£Z_|- 

In particular, for / G V we have 

(3-5) i x (Sf) = X(f). 

It is easy to check that the operators Lx, X G 0, form an abelian (purely odd) subalgebra of 
the Lie superalgebra Der 0, namely 

(3.6) [ix,iy] = ix°iy + iy°ix = 0. 

The Xze derivative along X G is the even derivation of f2 of degree 0, defined by Cartan's 
formula: 

(3.7) L x = [5, l x ] = 5 o l x + l x ° S ■ 
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In particular, for / G V = f2° we have, by (3.5), Lx(f) = i-x(3f) = X(f). Moreover, it 
immediately follows by the fact that 5 2 = that 

(3.8) [5, L x ] = 5 o L x - L x o 5 = . 
We next want to prove the following: 

(3.9) [l x ,Ly] = lx ° Ly - L Y o L X = l>[X,Y] • 

It is clear by degree considerations that both sides of (3.9) act as zero on Cl° = V. Moreover, it 
follows by (3.5) that [t x ,Ly](Sf) = LxStySf - L Y 5L X Sf = X(Y(f)) - Y{X(f)) = [X,Y](f) = 
t[x,y](5/) for every / £ V. Equation (3.9) then follows by the fact that both sides are even 
derivations of the wedge product in Q. Finally, as an immediate consequence of equations (3.8) 
and (3.9), we get that 

(3.10) [L X ,L Y ] = L x o Ly - L Y o L x = L[ X ,y] ■ 

We then let vr(5g) = S, ir(X) = Lx, n(X£) = tx, for X e g, and equations (3.6), (3.7), (3.8), 
(3.9) and (3.10) exactly mean that ir : q — > DerO is a Lie superalgebra homomorphism. Hence 
the basic de Rham complex fHs a g-complex. 

Remark 3.1. To each differential /c-form Q E Cl k we can associate a V- multilinear skew-symmetric 
map w:gx...xg^Vby letting 

k 

k 



= ^sign(.)n^^ 



aU h ° U h ^S k r=l 



J m r ,n crr . ? 



and extending it by V-linearity in each argument. We can then identify the complex Q k with 
the space of such maps Homy(AvS' V). The action of the differential 5 on Q G £l k is then given 
by the usual formula for the Lie algebra cohomology differential: 

(3.11) Soj(X 1 ,...,X k+1 )= Yl (~l) i+1 X i {Q(X 1 ,-*;X k+1 )) 

l<i<k+l 

+ J2(-l) i+j cj([X i ,X j ],X 1 ,- l --- J --,X k+1 ), 

i<j 

the corresponding formula for the contraction operator Lx is 

(3.12) (L X Cj)(X u ...,X k _ 1 )=Co(X,X 1 ,...,X k _ 1 ) if fc > 1, 
and that for the Lie derivative Lx is 

k 

(3.13) (L x Cj)(X u ...,X k )=X(Q(X u ...,X k ))-^ t u(X 1 ,--- ,[X,X r },--- ,X k ). 

r=l 
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As in Section 1, we denote by g" the subalgebra of consisting of evolutionary vector fields, 
namely the centralizer of d G g, given by (1.1). By an abuse of notation we also denote by d 
the Lie derivative Lg : Cl k — ► fl k . Since d commutes with 5, we may consider the corresponding 
reduced g^-complex Q = Q,/d£l = @kez+ Q k - This is known as the reduced de Rham complex 
£1 = O(V), or the complex of variational calculus, or the variational complex. 

Elements u G = tl k /dCl k are called local k-forms. In particular Q° = V/dV is the space 
of local Junctionals (cf. Definition 1.25). By an abuse of notation, we denote by 5 and, for 
X G g d , by lx, Lx, the maps induced on the quotient space £l k by the corresponding maps on 
£l k . In particular Lq acts as zero on the whole complex Q,. The contraction operators induce 
an f2°-valued pairing between Q 1 and g d : 

(3.14) (X,lo) = (u,X) = t x {uj) , for X G g d , ueSl 1 . 

We have the following identity, valid in any g-complex, which can be easily checked using 
equations (3.7) and (3.9): 

L [X,Y] = Lx^Y — LyiX + St<Yt>X ~ lYl>X& ■ 

In particular, for u G ft 1 and 1,7 6 g 9 , we get the following identity, which will be used later: 

(3.15) (u,[X,Y])=L x (u,Y)-Ly(u,X)-iYLxSu. 
Another useful form of the above identity is, for uefi 1 , 

(3.16) (w, [X, Y]) = L x (u, Y) - (L x uj, Y) . 



3.3 Cohomology of Cl and Q. In this section we compute the cohomology of the complexes 
Q and f2 under two different assumptions on the algebra of differential functions V. In Theorem 
3.2 we assume that V is normal, while in Theorem 3.5 we assume that the degree evolutionary 
vector field A is diagonalizable on V. 

We first need to extend filtration (1.7) of V to the whole complex Cl. For n G Z + and % G I, 
we let 

(3.17) n k i = \u>eti k i a Cj = L_§_uj = 0, for all (m,j) > (n,i) \ , 

3 3 

where the inequality is understood in the lexicographic order. We let Q k = ^ for n > 1, 
and S7q,o = ^fc,oC Clearly 0° 4 = V n ,i- Note that, consists of finite sums of the form 

(3.18) s= e cr^A-A^r'. cr^.- 

il,— mi, ■■■ ,rn k eZ+ 

(n,i)>(mi,ii)>—>(m k ,i k ) 

Given i G /, n G Z+, we have t a (ft*,) C ft*" 1 , L a (0*,) C and «5(^) C f^J 1 . 
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Explicitly, if <£> is as in (3.18), we have 



i a \oo) 

«.,(») 



il,*" ,i k (zl, mi,"' ,mj.GZ+ 
(n,i)=(mi,ii)>— >{m k ,i k ) 



p. rm\---m k 

(s.m^) = y: J f&r 5 < x) a • • • a ^r fc) e ^ ■ 

s "i i 1 ,-,i k eI,mi,-,m k eZ+ au i 

(n,i)>(mi,ii)>— >{m k ,i k ) 

flrmi—m k 

s(u) y y Jii -;\ 5uf ] a K (mi) a • • • a <y m&) G • 

v y / / / / p. (m) J *i *fc »V 

(n,i)>(m 1 ,i 1 )>--->(m k ,i k ) (n,i)>(m,j) 

Suppose next that the algebra of differential functions V is normal. Recall that this means that, 
for / G V n i, there exists g G V n j, denoted by f duf 1 ^ f, such that — %y = /; the antiderivative 

fdvf- f is defined up to adding elements from Ker ^ = V n ,i-i- Now we introduce "local" 
homotopy operators h n ^ : f2^ i — ► ^j" 1 by the following formula 

(3.20) km = £ (/*4 B) C^D ^f 2) a • • • a G 



(n) #mi -m^ ^,(m 2 ) A A ^,(m fe ) c fifc-l 

i ' 

il,-" ,i k £l ,mi,--- ,m k £Z+ 
(n,i)=(mi,ii)>—>(m k ,i k ) 



where u is as in (3.18). The element (3.20) is defined up to adding elements from It is 

easy to check that 

(3.21) i a {li„ ,{Co)) = , L_a_(h n! i(u>)) = t_e (&) . 

Theorem 3.2. Lei V be a normal algebra of differential functions. 

(a) For Co G £l k L { , we have h n ^{8G)) + 5h nt i(Co) — Co G 0, k li _ 1 . 

(b) H k (n,5) = S kfi C, 

(c) H k (n,5) = 5 kfi C. 

Proof. Clearly all three elements h n ^(8u)), Sh ni i(Co) and Co lie in f^j- Hence, to prove (a) we 
just have to check the following two identities: 

i a (h n i (5Co) + 8h n i (uj) - Co G ,_] ) = 0, 

on: 

L a (hn^SQ) + Shn^u) - U) G fl* ,-_*,) = 0. 

This is a straightforward computation using (3.19), (3.20) and (3.21). Part (b) immediately 
follows from (a) by induction on the lexicographic order. We next claim that 

(3.22) H k (dtt) = , for all k > . 
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Indeed, if 5{duj) = 0, then 6u = since 5 and d commute and Ker <9|^fc+i = for k > 0. Hence, 
by (b), uj = 5fj, where fj G Cl k ~ 1 , if k > 1, and u G C if fe = 0. In both cases we have <9o) G <5(<90), 
thus proving (3.22). Next, consider the following short exact sequence of complexes: 

This leads to the long exact sequence in cohomology 

• •• -^H k (dn) -> F fe (o) -» # fc (o) H k+1 (dn) ■•• 

By part (b) and equation (3.22), we get H k (Q,S) ~ H k (£l,5), for all fc G Z + , proving (c). □ 

Remark 3.3. It follows from the proof of Theorem 3.2 that, for any algebra of differential 
functions V, given a closed £;-cocycle, we can write it as Slj, where uj is a k — 1-cocyle with 
coefficients in an extension V of V obtained by adding finitely many antiderivatives. 
Remark 3.4. Theorem 3.2 holds in the more general setup of [DSK] where an algebra of dif- 
ferential functions is defined as a differential algebra V with commuting partial derivatives 
- , i G I, n G Z+, such that all but finitely many of them annihilate any given element of V, 

and the commutation relations (1.2) hold. In this case, in parts (b) and (c) of this theorem, 
one should replace C by C/(C n dV). 

We next construct a "global" homotopy operator, different from (3.20), based on the as- 
sumption that the degree evolutionary vector field A, defined by (1.21), is diagonalizable on 
V. 

Let V = ® Q V[a] be the eigenspace decomposition of V with respect to the diagonalizable 
operator A. Then the Lie derivative La is diagonalizable on £l k for every k G Z+. In fact, 
we have the eigenspace decomposition £l k = [a + k], where £l k [a + k] is the linear 

span of elements of the form f5u™ A • • • A Suf" 1 ^ with / G V[a]. We define the operator 
La-i : fi fc — > fi fc , letting it act as zero on the subspace O fc [0] = Ker LaU^ and as the inverse 
of La on = ©Q^.fc" [a + k]. Equivalently, La-i is uniquely defined by the following 
equations: 

(3.23) La o La-i = La- 1 ° La = 1 — 7ro , 

where ttq : £l k — > £l k denotes the projection onto £7 fc [0]. We also notice that 6 preserves the 
degree of a differential form u), hence it commutes with both La and La-i- Moreover, since 
A is an evolutionary vector field, the total derivative d commutes with ^a. La and L^-i. 
In particular, we have an induced eigenspace decomposition for A on the reduced complex, 

n k = ® a n k [ a + k}. 

Theorem 3.5. Let V be an algebra of differential functions, and assume that the degree evo- 
lutionary vector field A is diagonalizable on V. Consider the map h = La-i ° la '■ & k ~ ► 
£i k ~ l , k > 0. It satisfies the following equation: 

(3.24) h(5u>) + 5h(u>) = u> - 7r u> , u €&. 
Hence 

H k {n) ~ Ker (5 : n k [0] -> fl k+1 {0}) /^^[O] , for all k G Z+ . 

In particular, if V[0] = C and A does not have negative integer eigenvalues on V, we have 
H k (n) ~ 5 kj0 C. 
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Proof. Equation (3.24) is immediate from the definition of h, the fact that 5 commutes with 
L A -i and equations (3.7) and (3.23). For the second statement of the theorem, we notice 
that, since 5 and La commute, we have the corresponding A-eigenspace decomposition in the 
cohomology: 

H k (Q) = ^H k (Q[a + k}), 

a 

and we only need to show that H k (Q[a + k]) = when a + k ^ 0. For this, let u> S Q k [a + k], 
with a + k / 0, be such that 5uj = 0, and let u £ £l k [a + k] be a representative of u>, so that 
5uj = dfj for some fj G fi fc+1 [a + &]. Since 9 commutes with h, by (3.24) we have 

uj = h{5Cb) + 5h{uj) = dh(fj) + Sh(u) = 5h(uj) mod 317 . 

Hence cj € Im <5, as we wanted. □ 



3.4 The variational complex as a Lie conformal algebra cohomology complex. 

Let us review, following [BKV], the definition of the basic and reduced cohomology complexes 
associated to a Lie conformal algebra A and an ^-module M. A k-cochain of A with coefficients 
in M is a C-linear map 

7: A® fc ^M[Ai,...,A fc ] , oi®---®o fc i-^7A lj „.,A fc (oi > ...,a fc ) J 

satisfying the following two conditions: 

1- 7Ai,...,A fc (oi, • • • ,da i: ...,a k ) = -\a\ u ... ) \ k (a 1 , . . .,a k ) for all i, 

2. 7 is skew-symmetric with respect to simultaneous permutations of the a^s and the A,;'s. 

We let t k (A, M) be the space of all fc-cochains, and t(A, M) = fc>o f k (A, M). The differential 
5 of a fc-cochain 7 is defined by the following formula (cf. (3.11)): 

(3.25) (5^) Xl ,...,x k+1 (a 1 ,. . . ,a k+1 ) = V (-l) l+1 a iAi 7 * (o 1} • r -,a fc+ i) 

l<i<fc+l Ai,...,A fc+1 

x — > . . ■ « i 

+ 2^(- 1 ) 7 u ([aiA l «j],«i,- r -- : -,Ofc+i)- 

j<j Ai+Aj,Ai, ,A fc+ i 

One checks that 5 maps f k (A,M) to f fc+1 (^,M), and that 5 2 = 0. The space f(A,M) with 
the differential S is called the tesic cohomology complex associated to A and M. 
Define the structure of a C[<9]-module on T(A,M) by letting 

(3.26) (c^)A 1 ,...,A fc (ai,---,afc) = (9 + Ai H h Afc)7A 1 ,...,A fc (ai, • • • ■ 

One checks that 5 and 9 commute, and therefore dT(A, M) C T(^4, M) is a subcomplex, and 
we can consider the reduced cohomology complex T(A,M) = T (A, M) / dT (A, M) . 

For us, the interesting example is A = ©j 6 / C[<9]u;, with the zero A-bracket, and M = V, 
where V is an algebra of differential functions in the variables iti, i £ I, and the representation 
of A on M is given by the Beltrami A-bracket, as in (1.66). In this case we consider the subspace 
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ffin = ©fcez -^fin ^- V), where fg n consists of /c-cochains 7 with finite support, namely 
such that — ,\k( u hi " ' 1 u ik) = ® f° r an Du ^ finitely many choices of the indices ij, . . . , if. € /. 
Clearly, Tfi n is a subcomplex of r(A, V), and it is a C[<9]-submodule. Hence, we can consider 
the associated reduced complex rg n = T{[ n /dr{i n . The following theorem was proved in [DSK] 
(the first part of Theorem 4.6), in the case when / is finite. However the proof is the same for 
arbitrary I. 

Theorem 3.6. The map 4> : Tfi n — > f2, given by 

(3.27) m = I £ ft^ A -.. A ft^> , 

where f™*" i 'J nk £ V is the coefficient of A™ 1 • • • A™* m 7Ai,...,A fc ("iij ■ ••)*%)> * s bijective, and it 
commutes with the actions of 5 and d. Hence it induces isomorphisms of complexes £1 — ► Tfi n 
and Q — > r^ n . 

3.5 The variational complex as a complex of poly-A-brackets or poly-differential 
operators. Consider, as in Section 3.4, the Lie conformal algebra A = Q ieI C[d]v,i, with trivial 
A-bracket, and recall the Beltrami A-bracket (1.62) on an algebra of differential functions V. 
Recall from [DSK] that, for k > 1, a k-X-bracket on A with coefficients in V is, by definition, a 
C-linear map c : A® k -> C[Ai, . . . , A fc _i] ® V, denoted by 

a\ <g> • • • (g> a fc 1 ^ {a iAl • • • ak-i Xki a k } c , 

satisfying the following conditions: 

Bl. {a 1Xl • • • {da,i) Xi ■ ■ ■ a fc _ Ufc _ i a fc } c = -Ai{a Ul • • • a fc _ Ufc _ i a fc } c , for 1 < i < k - 1; 

B2. {a 1Xl • • • a fe _i Afc _ 1 (<9a fc )} c = (Ai H h A fe _i + <9){ai Al • • • afc-i^a^d 

B3. c is skew-symmetric with respect to simultaneous permutations of the aj's and the Aj's in 
the sense that, for every permutation a of the indices {1, . . . , k}, we have: 

{a iAl • • • a k . lXk i a k } c = sign(<r)K (1) ^ ■ ■ • ^-D^M*)}' AfcH , A t ■ 

The notation in the RHS means that A^ is replaced by X k = — X)j=i ~~ ^ ^ occurs, 
and d is moved to the left. 

We let C° = V/dV and, for k > 1, we denote by C fc the space of all fc-A-brackets on A with 
coefficients in V. For example, C 1 is the space of all C[<9]-module homomorphisms c : A — > V 
(which in turn is isomorphic to V ). We let C = © fcgZ+ C fc , the space of all poly X-brackets. 

We also consider the subspace Cfi n = © fcgZ+ Cg n C C, where Cg n = C° and, for A; > 1, Cg n 
consists of fc-A-brackets c with finite support, namely such that {^£ lAl ■ ■ ■ Ui k _ ix ^ Ui k }c = for 

all but finitely many choices of the indices ij, . . . , i& £ I. For example, Cg n ~ V®^. 

We next define, following [DSK], a differential d on the space C of poly A-brackets such that 
d{C k ) C (7 fc+1 and d 2 = 0, thus making C a cohomology complex, and such that Cfi n C C is a 
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subcomplex. The difference with [DSK] is that here we use the Beltrami A-bracket in place of 
the equivalent language of left modules defined by the Beltrami A-bracket. 

For J f G C° = C5 = V/dV, we let df f £ C 1 be the following C[<9]-module homomorphism: 

(3.28) (dff) (a) ( = {a} d j f ) := {a„ d f} B . 

For c G C k , with k > 1, we let dc G C fc+1 be the following A; + 1-A-bracket: 

k . 

{fllAi ' ' ' a kX k a k+l}dc '■= ^2(-^Y +1 ^ a iX,{ a lXi a kX k a k+l} ^ B 

i=l 

(3-29) +(-l) fe {{a 1Al ---a fc „ lAfc i a fc } c Al+ ... +Afc a fc+ i} B . 

For example, for a 1-cocycle, namely a C[<9]-module homomorphism c : A — > V, we have 
(3.30) {a A 6} rfc = {a x c(b)} B - {c{a)\b} B . 

We define, for fe > 1, a C-linear map V fc : F fc — > C k , as follows. Given 7 G T fc , we define 
^) . A®*->C[Ai,...,A fc _i]®V,ty: 

(3-31) {a Ul ■■■a fe -lA fe _ 1 a fe}^ fc (7) = 7 Ali ... |Afe _ liA t(air-- ' afc )' 

where, as before, \' k = — Sj=i Aj ~~ ^ an d <9 is moved to the left. The following theorem is a 
special case of Theorem 1.5 from [DSK]. 

Theorem 3.7. The identity map on V/dV and the maps tj) , k > 1, induce isomorphisms of 
cohomology complexes T — > C and Tfi n — > Cfi n . 

Due to Theorems 3.6 and 3.7, we can identify the variational complex f2 with the complex 
of poly-A-brackets Cfj n . Explicitly, we have: 

Corollary 3.8. Let k > 1 and let u G (l k be as in (3.2), where we assume that the coefficients 
f™}" i ' mk are skew-symmetric with respect to simultaneous permutations of the lower and upper 
indices. Let p k (oj) be the k-X-bracket defined by 

(3.32) {u ilXl ---u ik _ lXk _u ik }^ ) =k\ Y, Ar---A™V(-Ai-----A fc -i-9rV™. 1 ,r' 

mi,— 

and extended to A® k by sesquilinearity. Then, the identity map on V/dV and the maps p k : 
£l k — > C k , k > 1 factor through an isomorphism of complexes: p : O — > Cfi n . 

Proof. It follows immediately from the definitions (3.27) of <f> and (3.31) of ip k . □ 

Next, for Xp G Q 9 , P G V e , we define, the contraction operator lx p '■ Cg n — > Cg" 1 , > 1, 
making p an isomorphism of g^-complexes. For c G Cg n , we let 

(3.33) l Xp {c) = fY,Pic( Ui ) G V/dV = C° . 
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For c G Cg n with k > 2, we let 

(3.34) {a 2 A 2 • • • a k -i Xk i a k } bXp{fi) = ^ {n i9 a 2 A 2 • • • a k -i Xh _ 1 a k } c P , 

where, as usual, the arrow in the RHS means that d is moved to the right. The following 
proposition is a more precise form of the second part of Theorem 4.6 from [DSK]. 

Proposition 3.9. The contraction operators lx p defined by (3.33) and (3.34) endow the coho- 
mology complex Cfi n with a g d -structure, and the isomorphism of complexes p : $7 — > Cfi n (de- 
fined in Corollary 3.8) is an isomorphism of Q d -complexes, namely p k ~ 1 (ix P {oj)) = ix P {p k {^)) 
for all wetf and Xp G g d . 

Proof. If uj G fT is as in (3.2), with f™]l " i ' mk skew-symmetric with respect to simultaneous 
permutations of lower and upper indices, we have, 

t XP (Q)=k y: ( e c':;rM 5mi ^))^r 2) A---A^. 

i2,...,i k £l i 1 eI,mi€Z+ 

Hence, applying equation (3.32), we get 

{u i2M ■ ■ ■ Ui k _ 1Xk _ i U ik } p k-i( liXp (Cj)) 

= m 53 xr ■ ■ ■ a:_v (-a 2 x k -i - d) m * (./;:::;:"" (^j) . 

mi,— ,m k ez + 
hei 

On the other hand, by (3.32) and (3.34), we have 

{ u t2\ 2 ' ' ' Ui k-l\ k _ 1 U ik}t, Xp (p k (u)) = E^ Uil 9 ni2 A2 ' ' ' Ui k-l\ k _ 1 U ik}p k (u) ->Pil 

hei 

= k\ 53 x ™ 2 ■ ■ ■ C fc r 1 (-A 2 \ k -i - d) mk (/::':;;;"" (a mi P n )) . 

mi,— ,rrifeGZ+ 

hei 

□ 

The space C k of fc-A-brackets, for k > 1, can be equivalently described in terms of maps 
S : (V^)^" 1 — ► V®^ of differential type, namely of the form: 

(3.35) ^P 1 ,..-^- 1 )^ = ^ C'X^^ 1 ^)---^ 1 ^)' 

mi,— ,m fe _ 1 eZ + 

n,- i*fc-ie/ 

where / 4i j'"^ 1 £ V is non-zero only for finitely many choices of the upper and lower indices. 
Indeed, given a /c-A-bracket c G C fc , the corresponding map of differential type is (3.35), where 
fh - tfc *7fc * s ^ ne coefficient of A™ 1 • • • Aj.^ 1 in {ui lXi ■ ■ ■ Ui k _ ix ^ iH k }c- It is easy to see that 
this gives a bijective correspondence between the space C k of fc-A-brackets and the space of 
maps S : (V^) fc_1 — > V®^ of differential type which are skew-symmetric, namely such that 

(3.36) fS{P a{1 \ . . . 5 i^(fc-i))p<fc) = S ign(a)/5(P 1 , . . . , P k ~ 1 )P k , 
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for every P 1 , • • • , P k £ V e and every permutation a. 

An equivalent language is that of skew-symmetric poly-differential operators. Recall from 
[DSK] that a k- differential operator is a C-linear map 5 : (V i ) k — ► V/dV, of the form 

(3.37) sop 1 , • • • ,p fc ) = | 53 /;::::;;" f^-/^ • • • (^j*) . 

The operator S is called skew-symmetric if 

J S(P\--- ,P fe ) = sign(a) | S(P CT « ••• ,P CT W), 

for every P 1 , • • ■ ,P fc G V and every permutation a. Given a A:-A-bracket c £ C fc we associate 
to it the following ^-differential operator: 

(3.38) s(p\- ,p k ) = J yi cx^ 1 ^)---^ 1 ^:*)^ 

ni,- ,n fe _ 1 eZ+ 

where /™!.'. jlTi it ^ s as m (3-35). It is easy to see that the skew-symmetry property of the 
/c-A-bracket is translated to the skew-symmetry of the /c-differential operator. Conversely, inte- 
grating by parts, any ^-differential operator can be written in the form (3.38). Thus we have 
a surjective map from the space of fc-symbols to the space S fc of skew-symmetric ^-differential 
operators. By Proposition 1.3, the ^-differential operator S can be written uniquely in the form 

(3.38) , hence this map is an isomorphism. 

The differential, the contraction operators, and the Lie derivatives for the variational com- 
plex f2 have an especially nice form when written in terms of poly-differential operators. The 
following result was obtained in [DSK, Section 4.5]. 

Theorem 3.10. The variational complex (17 = £l k ,S) is isomorphic to the complex of 

poly-differential operators (S = @fc g x + ^ k ->d), where 

k+i 

(3.39) (dS)(P\--- ,P fc+1 ) = 53(- 1 ) S+1 ( X ^* S )( pl '---' pfc+1 )- 

s=l 

In this formula, if S is as in (3.37), XpS denotes the k- differential operator obtained from S 
by replacing the coefficients /"*.'.'.' by X P (f^,['^ k ). 

Furthermore, the -structure on £1 is translated to the following q® -structure on E: 

(3.40) (i Xpl S)(P 2 ,--- ,P k ) = S(P\P 2 ,--- ,P k ), 

k 

(3.41) (L Xq S)(P\--- ,P k ) = (X Q S)(P\--- ,P k ) + Y,S(P 1 ,--- ,A>(Q),-.. ,P k ). 

s=l 

Equation (3.39) is similar to formula (4.9) in [D]. 
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We can write explicitly the bijection between the spaces and E, which gives the isomor- 
phism of g^-complexes. If u G the corresponding ^-differential operator in G E fc is given 
by 

(3.42) sup 1 , ■■■ ,P k ) = ix pk ■ ■ ■ lx p1 u g n° . 

This is clear by applying (3.40) iteratively. Due to the identification (3.42) of O with E, we can 
think of an element oj & O as a map cu : given by 

(3.43) w(X!, ■ • • = i Xk ■ ■ ■ i Xl u G ^° • 

Clearly, E° = V/dV. Moreover, for k > 1, using the bijections between the space E fc of 
^-differential operators and the space of maps (V^) fc_1 — ► V®^ of differential type, we get that 
E 1 ~ V®^, and E 2 is naturally identified with the space of skew-adjoint differential operators 
S : V — > V® . With these identifications, d, ix P and L Xp , for P G V , are as follows. For 
// G S° = V/aV we have 

(3.44) 4/7) = g G ye*, ^ p (//)=0, L Xp {Jf) = fX P (f) eV/dV. 

For F G E 1 = V® £ , we have 

(3.45) d(F) = £> F (d) " ^f(^) G s2 > ^p(P) = / E ie / ^ G V/9V , 

(3.46) Lx P (dF) = D F (d)P - D* F (d)P £V ee , d{i Xp F) = D* F (d)P + D* P {S)F G V®^ , 

(3.47) L Xp (F) = D F (a)P + D* P (d)F G V® £ , 

where Dp(d) is the Frechet derivative (1.13). Finally, a skew-adjoint operator S(d) = (Sij(d))ij^i 
G E 2 is closed if and only if the following equation holds (cf. (0.28)): 



(3.48) {ui^S kj (X)} B - {uj \S ki (n)} B = {Sij(X) \ +/J ,u k } 
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Definition 3.11. A matrix valued differential operator S(d) : V — > V®^ is called symplectic 
if it is skew-adjoint and satisfies equation (3.48). 

Propositions 1.5 and 1.9 provide a description of the kernel and the image of the variational 
derivative for normal V, based on Theorem 3.2, namely the exactness of the variational complex 
at k = and 1. The next corollary of this theorem provides a description of symplectic 
operators, and it is based on the exactness of the variational complex at k = 2. 

Proposition 3.12. If V is a normal algebra of differential functions, then any symplectic 
differential operator is of the form: Sf(0) = Dp(d) — D* F {d), for some F G V®^. Moreover, 
Sf = Sg if and only if F — G = |^ for some f G V. 

Example 3.13. The simplest symplectic differential operators in the case t = 1 are: 
Si u(n) {d) = d n for n odd , S uu ,(d) =u' + 2nd, Si u #(d) = u" + 2u'd. 
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Example 3.14. Let V be the localization of the algebra of differential polynomials C[u, u', • ■ ■ ] 
at u' . Then, for F = — ^7 6 V, the corresponding symplectic operator 

S F (d) = ± 7 do± 
u u 

is the Sokolov symplectic operator. For G = — <9 2 ^^7^ G V, the corresponding symplectic 
operator 

s G (d) = do^-do^-d 

u u 

is the Dorfman symplectic operator. As we will see in Section 4, these operators play an 
important role in the integrability of the Krichever-Novikov equation. 



3.6 Projection operators. Another way to understand the variational complex uses projec- 
tion operators Vk, k > 1, of £l k onto its subspace complementary to d£l k . We can then identify 
the space of local A;- forms Q k with the image of Vk in 

For i £ / we define the operators U : Q k — > fi* - k € Z+, as 

(3.49) I, = (~9) n ° * fl/ft4 n) , 

and, for k > 1, we let "P^ : fc — > f2 fc be given by 

(3.50) P fc (cD) = ^ A ^ • 

iei 

This formula expresses the fact that, integrating by parts, one can bring any differential form 
Co £ £l k uniquely in the image of V ■ The projection operators V k on the space of differential 
forms f2 are analogous to the projection operators on the space of polyvectors, introduced in 
[B]. 

The following theorem shows that Vk is a projection operator onto a subspace Im (Vk) C Q k 
complementary to dQ k , for k > 1. 

Theorem 3.15. Let k > 1. 

(%) 7j o 9 = /or aZZ i £ I. In particular Vk d = 0. 
(t»; 7/w€ O fc , t/iera V h (Cb) -Co £ d& k . 
(Hi) Ker V k = dVt k and V k 2 = V k - 

Proof. By (1.2) and (3.9) we have that [ L g/g u (n),d\ = l q/q u (-"-i), where the RHS is stipulated 
to be zero for n = 0. It follows that 

u o d = (-or ° i eM * ° d = E(~ d r ° - £ (-^) n+1 = > 

n>l nGZ + 
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thus proving (i). Notice that, if Co G £l k , we have 

(3.51) ]T 6u^AL d/du M(u) = kCj. 

i£l,n£Z + 

On the other hand, since d commutes with 5 and it is a derivation of the wedge product in f2, 
we have that 

(3.52) ]T 5ut ] M d/du ^{u)= S Ui A(-drt d/du ^(u;) mod dtl k . 

Part (ii) is then an immediate consequence of equations (3.51) and (3.52). Finally, (hi) follows 
from (i) and (ii). □ 

Let cD = E h,.,i k a f""^^ A ' ' ' G where the coefficients J'^X'" G V 

mi....,mfeGZ_|_ 

are skew-symmetric with respect to exchanging simultaneously the lower and the upper indices. 
We have 

(3.53) vm = £ (-ar* (fz^ 5 ^ a • • • a Ktr i} ) a <k • 

We associate to the RHS of (3.53) the following map S : (V^)^'" 1 -> V® £ of differential type 
(defined in Section 3.5): 

(3.54) SiP 1 ,...,!*- 1 )^ = (-dr^f^ir^P^-'-id^Ptl)) ■ 

U...jfe_l£/ 
mi...mii€Z-|. 

It is immediate to check that S is skew-symmetric in the sense of (3.36). Hence, using the 
projection operators, we get a more direct identification of the spaces Q k and T, k . 

For example, if Co = ^2 ieI n£ z + fFfi u i e the corresponding element in is V\(Cj) = 
Yliei (Snez + ( — ff)^ u i- Therefore the projection operator gives an explicit identification 
of V® £ with fi 1 , which, to F G V® £ , associates 

(3.55) w F = J^FjfJui G P 1 ^ 1 ) ~ ft 1 . 

As pointed out in the previous section, the space of skew-symmetric maps of differential type 
S : (V^)^" 1 — ► V®^ is canonically identified with the space S fc of skew-symmetric /c-differential 
operators S : (V e ) k V/dV. In view of Remark 3.1, the differential 5 on S fc is given by 
formula (3.11), where Co G & k is replaced by S G S fc and Xj G fl 9 . Using formula (1.15) for 
the commutator of evolutionary vector fields, it is straightforward to rewrite 5 on S fc in the 
form (3.39) (cf. [D], p. 61). Likewise, it is immediate to derive (3.40) and (3.41) from (3.12) and 
(3.13) respectively. 
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Remark 3.16. We can use the isomorphism g" ~ V to conclude that the space of local fe-forms 
ft fc can be identified with the subspace of Homc(A fc_1 g 9 , ft 1 ) consisting of skew-symmetric (in 
the sense of (3.36)) linear maps S : A k ~ 1 g® — > ft 1 acting by differential operators (in the sense of 
(3.35)). In order to find the expression of the map S : A k ~ 1 g d — > ft 1 corresponding tow£ ft fc , 
we notice that equation (3.54) can be equivalently written as ^2 ieI S(P 1 , . . . , P k ~ 1 )i5u{ = 
kl L X p i " " " L x pk _ 1 ^- We thus conclude that the map 5 : A kl g 9 — > ft 1 corresponding to u G ft fc 
is given by 

(3.56) S : X x A • • ■ A X k ^ ^ —l Xi ■ ■ ■ tjf^w • 

Remark 3.17. Via the identifications g 9 ~ and ft 1 ~ V®^, the pairing (3.14) between g 9 
and ft 1 becomes (1.4). Since, by Proposition 1.3, this pairing is non-degenerate, we get the 
embedding ft 1 ■— > Homc(g 9 , V/dV), discussed in Section 3.5, which associates to w G 1 the 
corresponding map of differential type S : g® — > V/dV. However, this embedding is not an 
isomorphism. For example the linear map <j) E Homc(g s , V/dV) given by <j>(Xp) = Yliei I §IT' 
is not in the image of ft 1 . This is the same as saying that there is no element F G V® e such 
that fF-P = JYsi dPi/dm, for all P eV l . Indeed, by taking Pj = S itj , we get that Fj G 9V 

for all i, and by taking Pj = S id u[ N) , N > 1, we get that / i^uj^ = for all k G J, N > 1. By 
Lemma 1.30 this is possible only for F = 0, which is a contradiction. 

4 Dirac structures and the corresponding Hamiltonian equations. 
4.1 Derived brackets and the Courant-Dorfman product. 

Definition 4.1. Let L be a Lie superalgebra with the bracket [• , •]. Given an odd derivation 
d : L — ► L such that d 2 = 0, the corresponding derived bracket on L is defined by the following 
formula: 

(4.1) [a,b] d = (-l) 1 -^)^)^, 
where p(a) G Z/2Z denotes the parity of a G -L. 

This construction goes back to Koszul [KS]. It is not hard to check that the derived bracket 
satisfies the Jacobi identity with respect to the opposite parity on L, namely 

(4.2) [a, [b, c] d ] d - (-l)(i-ri«))a-p(6)) [ &) [ a? c y d = [[ 0j b]dj c]d . 

However, the derived bracket on L with the opposite parity does not define a Lie superalgebra 
structure since, in general, skew-commutativity fails. 

The relevant special case for us of the above construction is the following. Consider the 
space ft = ft(V) with the parity induced by its Z + -grading and let L = Endc(ft) be the Lie 
superalgebra of endomorphisms of the superspace ft. Then d = ad 5 is an odd derivation of L 
satisfying d 2 = 0, and we can consider the corresponding derived bracket [• , -J^ on L. 

Note that we have a natural embedding of ft 1 © g — > L obtained by Q i— > cDA , J k tj. 
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Proposition 4.2. (a) The image of ti 1 © g in L is closed under the derived bracket [• , -]d, 
and we have 

(4.3) (u © X) o (?) ®Y) = (L x (fj) - i y {5oj)) © [X, Y] , 

for every X,Y £ g and £j,fj G Cl 1 . 

(b) The bracket (4.3) induces a well-defined bracket on fi 1 ® g 9 . 

Proof. First notice that, since 5 is an odd derivation of the wedge product, we have d{£uA) = 
[<5,u)A] = <5(u))A for every Cj G fi , and, by equation (3.7), we have d(t x ) = [8,ix] = Lx for 
every J £ g. Hence we have, by the definition (4.1) of the derived bracket, that 

[uA,fjA] d = [S(u)A,fjA] = , [uA,L X ] d = [5(uj)A,i x ] = -ix(Su)A, 
[Lx,u)A] d = [L x ,tiA] = L x (u>)A , [ix,hr] d = [Lx^y] = L[x,y]- 

In the last identity we used equation (3.9). Equation (4.3) follows from the above relations, 
proving (a). For part (b), we just notice that Q 1 © g 9 C © q is a subalgebra with respect to 
the bracket (4.3), and dQ l © q 9 C n 1 © g 9 is an ideal in this subalgebra. □ 

4.2 Definition of a Dirac structure. Dirac structures were introduced independently by 
Courant and Dorfman (see [D] ) , as a generalization of the notions of Hamiltonian and symplectic 
operators. We will discuss their relation to these operators in Sections 4.7 and 4.8. 

Let V be an algebra of differential functions. Recall from Proposition 1.3 and equation 
(3.14) that we have a non-degenerate pairing g 9 x fi 1 — ► V/dV, given by 

(X P , up) = (lo f , X P ) = i Xp (u F ) = y^J PiFi ■ 

iei 

We extend it to a symmetric fi°-valued bilinear form ( , ) : (O 1 (Bg 9 ) x (fi 1 (&g 9 ) — > 0° = V/dV 
by letting J7 1 and g 9 be isotropic subspaces, i.e. (u>, rj) = (X, Y) = 0, for all u, n G Q 1 , X, Y G g 9 . 
In view of Proposition 4.2(b), we define the Courant- Dorfman product on © g 9 (with even 
parity) by 

(4.4) (to © X) o (7/ © Y) = (L x (v) ~ iy(8u)) © [X, Y] . 

By (4.2) and Proposition 4.2, this product satisfies the left Jacobi identity [A, [B, C]] — [B, [A, C]] 
= [[A, B],C], for all A, B, C G Q 1 © g 9 , but it is not skew-commutative. 

Definition 4.3. A Dirac structure is a subspace £ C fi 1 ® / which is maximal isotropic with 
respect to the bilinear form ( , ), 

(4.5) C = L L = { A G n l © g 9 \ (A, B) = , for all B G C] , 
and which is closed with respect to the Courant-Dorfman product, 

(4.6) A,B£C => AoBeC. 
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Remark 4.4. It is clear that condition (4.6) can be equivalently expressed, using the first con- 
dition (4.5), by the relation 

(4.7) (4oB,C)=0, for all A, B, C £ C . 

Proposition 4.5. Assuming (4.5), condition (4.7) is equivalent to the following "integrability 
condition" for a Dirac structure: 

(4.8) (L Xl W2,Xa) + (Lx a W3,X 1 ) + (Lx a wi,X 2 ) = , for all o>i ®X lt u 2 ®X 2 ,uj 3 ®X 3 £ C. 

Proof. In order to derive (4.8) from (4.7) we notice that, for A = u\ © X\, B = uj 2 © X 2 , C = 
oj 3 © X 3 £ Q, 1 © g 9 , we have 

(4.9) (A o B, C) = (L Xl (u 2 ),X 3 )- (ix^Su^Xs) + (w 3 , [X U X 2 ]) . 
The second term in the right hand side of (4.9) is 

-(lx 2 (8uji),X 3 ) = (lx 3 {5ux),X 2 ) = (L X3 ui,X 2 ) - L X2 (ix 3 (^\)) ■ 
Similarly, the last term in the right hand side of (4.9) can be rewritten as 

(w 3 , [Xi,X 2 ]) = -t[x 2 ,Xi]^3 = 1X1(^X3(^3)) - ix 2 (iXi(^3)) = (Lx 2 UJ 3 ,Xi) - Lx 2 (<Xi(^3)) • 

We can thus put the above equations together to get 

(AoB,C) = (L Xi uj 2 ,X 3 ) + (Lx^X^ + (L X3 0Ji,X 2 ) 
-L Xa {(ui,X 3 ) + (u 3 ,X 1 )). 

On the other hand, by the assumption A,B,C £ £ = C , we have (A, C) = 0, which precisely 
means that (u±,X 3 ) + (u 3 ,X%) = 0. This proves the equivalence of (4.7) and (4.8), assuming 
(4.5). □ 

Remark 4.6. Define the following important subalgebras of a Dirac structure C: 

Keri£ = {u © £ C\u £ n 1 } , Ker 2 £ = {0 © X £ C \ X £ g 9 } , 
5C = {5(J f) ©X££|//£l2°,X£ fl 9 }. 

Note that the Courant-Dorfman product on Keri£ is zero, and on 5C D Ker 2 C is as follows: 

(4.10) {5(J f) © X) o (5(J g) © Y) = 5fX(g) © [X, Y] . 
Consequently, (5(f f)@X)o (5(J g) © Y) £ Ker 2 £ if / X(#) = 0. 
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4.3 Hamiltonian functionals and Hamiltonian vector fields. Let C C fl 1 © Q d be a 

Dirac structure. We say that fh G Q° is a Hamiltonian functional, associated to a Hamiltonian 
vector field X, if 

6(fh)®X eC. 

We denote by .F(£) C fl° the subspace of all Hamiltonian functionals, and by TL(£) C g d the 
subspace of all Hamiltonian vector fields, i.e. 

;F(£) = j/fr G S(Jh) © X G C for some Xeg 8 }, 

H(£) = {x G g a © X G £ for some fh G V/dv} . 

Lemma 4.7. 71(C) C g 9 is a Lie subalgebra, and J-(C) C 0° is an 71(C) -submodule of the 
q 9 -module QP with action given by Lie derivatives. 

Proof. It is immediate from (4.10) and the fact that, by definition, a Dirac structure C is closed 
under the Courant-Dorfman product. □ 

Lemma 4.8. Let X,Y G 7i(C) be two Hamiltonian vector fields associated to the same Hamil- 
tonian functional Jf G F(C), i.e. S(ff) © X, 5(ff) ®Y G C. Then 

JX(g) = fY(g) , 

for every Hamiltonian functional j g G J~(C) . 

Proof. By assumption 5(f f) © X, 5(ff) © Y G C, and 5(fg) © Z G C for some Z G g 9 . In 
particular, they are pairwise orthogonal with respect to the bilinear form ( , ): 

(*(//) © X, 8(Jg) ®Z) = (6(Jf) © Y, S(fg) © Z) = , 

which means that j'X(g) = -f Z(f) = j'Y(g). □ 

Lemma 4.8 guarantees that the action of the Lie algebra 7i(C) on J~(C), given by Lemma 
4.7, induces a well-defined bracket 

{ , } c : F(C) x F{C) - , 

called the Lie bracket of Hamiltonian functionals, given by 

(4.11) {//Js^pfG?), H5(ff)®XeC. 

Proposition 4.9. T/ie bracket { , }x defines a structure of a Lie algebra on T(C). 

Proof. Skew-commutativity of { , }c follows by the fact that C C f^ 1 © is an isotropic 
subspace. Indeed, if S(f f) © X, 5(f g) © Y G C, we have 

{!f>fs} £ + {Sg,Sf} £ 

= JX(g) + fY(f) = (S(ff) © X, 5(J g) © Y) = . 
Suppose now <5(J* /j) © Xi G £, i = 1, 2, 3. We have 

{.//.. {/./2.j7: ! }J i . = /*i(*2(*0) = (L x MIh)lXi). 

Hence Jacobi identity for { , }c immediately follows from the integrability condition (4.8). □ 
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Let C C fi 1 © g 9 be a Dirac structure, and let Jh G be a Hamiltonian functional 
associated to the Hamiltonian vector field X G g 9 , i.e. S(Jh) © X £ C. The corresponding 
Hamiltonian evolution equation is, by definition, 

(4.12) !-*(„). 

An integral of motion for the evolution equation (4.12) is, by definition, a Hamiltonian functional 
/ / € .F(£) such that [jh, Jf} c ( = JX(f)) = 0. 

Remark 4.10. One can think about the RHS of (4.12) as <yjh, fj , but for a general Dirac 

structure such a notation can be misleading since, for / E V, the "bracket" { f h, f}c depends 
on the particular choice of the Hamiltonian vector field X associated to Jh. In other words, the 
analogue of Lemma 4.8 fails (it only holds under the sign of integration). On the other hand, as 
we will see in Section 4.7, if £ is the graph of a Hamiltonian map H : Vt 1 — > q 9 , then to every 
Hamiltonian functional jli e ^(C) there is a unique Hamiltonian vector field X = H( y 5{jK)) 
associated to it, hence the notation {Jh,f}c makes sense. 

Definition 4.11. Given an element S(Jh) © X G £, the evolution equation (4.12) is called 
integrable (with respect to the Dirac structure C) if there exists an infinite sequence 5(jh n ) © 
X n G £, including the element S(Jh) © X, such that all Jh n are integrals of motion for all 
evolution equations 

(4.13) ^ = Xn (u), 

all evolutionary vector fields commute, i.e. [X m ,X n ] = for all m,n G Z+, and the integrals 
of motion Jh n , respectively the vector fields X n , span an infinite-dimensional subspace in 
0°, respectively in q 9 . In this case the hierarchy (4.13) of evolution equations is also called 
integrable. It follows from the definition of a Dirac structure that the above sequence 5(J/i n ) © 
X n , n > 0, spans an isotropic subspace of C with zero Courant-Dorfman product. 



4.4 Pairs of compatible Dirac structures. The notion of compatibility of Dirac struc- 
tures was introduced by Gelfand and Dorfman [GD1],[D]. 

Given two Dirac structures C and we define the relation Nc,c C Q d © Q 9 by 

(4.14) M c ,c> = {X © X' G Q 9 © Q 9 | r] © X G C, r) © X' G £ for some n G O 1 } , 
and the adjoint relation ^, C ft 1 © fi 1 by 

(4.15) N£ )C , = {uj © J G n 1 © n 1 | (w, X) = (lu', X') for all X © X' G A/" £j£ /} . 

Definition 4.12. Two Dirac structures £, £' C VL l © g s are said to be compatible if for all 
X,X',Y,Y' G 9 , u,u/,u" G ft 1 such that X © X' , Y © F' G A/" £j£ /, w © w', J © w" G AA^ £ „ 
we have 

(4.16) (a;,[X,y])-( W ',[X,y'])-( W ',[X',y]) + (a;",[X',y']) = 0. 

Note, of course, that X © X' G Mc,c if an d only if X' © X G Afc,c, hence (£, £') is a pair 
of compatible Dirac structures if and only if (£',£) are compatible. 
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4.5 Lenard scheme of integrability for a pair of Dirac structures. Recall that we 
think of oj G Q k as a &;-linear skew-symmetric function uj(X\, ■ ■ ■ ,Xk), defined by (3.43). The 
following theorem is due to Dorfman [D]. 

Theorem 4.13. Let (£,£') be a pair of compatible Dirac structures. Suppose that Xq is 
a Hamiltonian vector field with respect to both C and C , that is, there exists functionals 
J ho, j hi G 0° such that 5(Jho) © Xq G C and S(Jhi) © Xq G £ . Assume moreover that: 

(i) ifuJGQ 1 is such that 5lo(X', Y') = for all X', Y' G tt 2 (M c ,c') (where vr 2 : g 9 © g 9 -> g 9 
denotes projection on the second factor), then 5uj = 0; 

(ii) there exists a sequence of 1-forms to n G fi , n = 0, 1, . . . , iV + 1 (N > 0, and it can be 
infinite), and of evolutionary vector fields X n G g 9 , n = 0,1,..., N (starting with the 
given Xq), such that ujq = S(fho), uj\ = S(Jhi) and 

0J n © X n G C , u Tl +i © X n G C! , < n < N . 

Then: 

(a) uj n © u n+ i G A/^ a for all n = 0, . . . , N + 1; 
^ w n is closed for all n = 0, . . . , N + 1, i.e. 5uj n = 0. 
Proof. By definition of Nc,c an d A/^ £/, (a) is equivalent to saying that, for < n < N, 

(uj n ,X) = (uj n+ i,X') , 

for all X, X' G g 9 such that 77 © X G £, 77 © X' G for some r\ G fi 1 . By assumption (ii) we 
have cj n © X n G £, w n+ i © X n G hence, since C and £' are isotropic, we get 

(uj n ,X) = -(j],X n ) = (u n+ i,X') , 

as we wanted. We next prove (b) by induction on n > 0. By assumption u n = 6(J f n ) for 
n = 0, 1. Furthermore £,£' is a pair of compatible Dirac structures, namely (4.16) holds. In 
particular, since by (a) us n © oJ n +i, ^>n+i © <^n+2 G A/^ £/, we have, by Definition 4.12, that 

(4.17) {oj n , [X, Y]) - {u n+1 , [X, Y']) - (u n+1 , [X', Y}) + (u; n+2 , [X 1 , Y'}) = , 

for all X(B X', Y(BY' G Mc,c- Using equation (3.15) and recalling the notation (3.43), equation 

(4.17) can be written as follows: 

- {5uj n ){X,Y) + (5io n+l )(X,Y') + (5io n+1 )(X',Y) - (5u n+2 )(X' ,Y') 

(4.18) + L x (u n ,Y) - L Y (u n ,X) - L x (u n+1 ,Y') + Ly,(u n+1 ,X) 

- L X '(uj n+ i,Y) + Ly(w n+ i, X') + L X '(u n +2,Y') - L Y >(uJn+2, X') = 0. 

Since uj n © uj n +i, u n +i © w n +2 G N£ c , and X © X', Y © Y' G Nc,c, we have 

(U! n , X) = (iU n+ i,X') , (LO n+ l,X) = (u n +2,X f ) , 

(u n , Y) = (un+i, Y') , (u n+ i,Y) = (uj n+2 , Y') , 

hence the last eight terms in the LHS of (4.18) cancel. Moreover, by the inductive assumption 
uj n and uJn+i are closed, so that equation (4.18) reduces to (Su n+ 2)(X' , Y') = 0, which, by 
assumption (i), gives 8u n+ 2 = 0. □ 
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Remark 4.14. Let (£,£') be a pair of compatible Dirac structures. Suppose that Tti(C) = fi 1 . 
In this case ^{■N'cc) = 7T 2(£')- Hence the non-degeneracy condition (i) of Theorem 4.13 reads: 
if 8u(X', Y') = for every X 1 , Y< G tt 2 (C'), then Su = 0. 

The following definition and proposition are a generalization, in the context of Dirac struc- 
tures, of Definition 2.4, Lemma 2.6 and Proposition 2.9. 

Definition 4.15. Let C and £ be two subspaces of ft 1 (Bg 9 - Two sequences {w n } n= o,... ,j\r+i C 
Q, 1 and {^n}n=o,— ,N C form an {C, C 1 )- sequence if 

w n © X n G C and w n +i © X„ G £' for n = 0, ■ • • , N . 

Proposition 4.16. Let L and £ be isotropic subspaces of&Qg 8 ', and suppose that {oj n }n=o,--- ,N+i 
C SI 1 and LXn} n= o,...,jv C 9 /orm an (£,£') -sequence. 

^aj T/ien we have: (uj m ,X n ) = /or all m = 0, • • • , iV + 1 and n = 0, • • • ,N. 

(b) Suppose in addition that 

(4.19) (Span c {uj n }o<n<N+i) 1 ~ C 7r 2 (£') , (5pan c {X n }o<n<iv) C 7ri(£) , 

where 7Ti and 7T2 denote projections on J7 1 and g a respectively, and the orthogonal com- 
plement is with respect to the pairing of and g®. Then we can extend this sequence to 
an infinite (C, C')- sequence {u n }n& + , {X n } n( zi + . 

Proof. Since C and C are isotropic, we have 

K, X n ) = for all n = 0, • • • , N , (w m , I n ) = -(w„,, X m ) for all m, n = 0, • • • , N , 
(<jj n+ i,X n ) = for all n = 0, • • • , JV , (w m +i, X re ) = -(w n+ i,I ra ) for all m, n = 0, ■ ■ ■ , iV . 

From the above relations we get 

{u} m ,X n ) = (u m +i,X n -i) , 

hence part (a) follows by an easy induction argument, as in the proof of Lemma 2.6. We then 
have, by (a), 

G (Span c {X n } <n<7v) , 

so that, using the second inclusion in (4.19), we find Xjy + i G Q 9 such that u>n + i © Xjy + i G C. 
On the other hand, we also have 

Xn+i G (Spancl^nlcKn^TV+i)" 1 , 

and using the first inclusion in (4.19) we find an element lun+2 G Q 1 such that u>n+2®Xn + i G C, 
and (b) follows by induction. □ 

Corollary 4.17. Let V be a normal algebra of differential functions, and let (£,£') be a pair 
of compatible Dirac structures, which is non- degenerate in the sense of condition (i) of The- 
orem 4-13. Let, for some N > 0, {u n }n=o,--- ,n+i C and {X re }n=o,- ,N C g d be such that 
loq = 5(Jho), = S(Jhi) G 5^1°, uj n © X n G C, © X n G £' /or all n = 0, • • • , N , and 
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conditions (4.19) hold. Then there exist a sequence of Hamiltonian functionals {Jh n } n ^z + 
such that 5(Jh n ) = uo n for n = 0, • • • ,N + 1, and a sequence of Hamiltonian vector fields 
{X n } ne z + extending the given one for n < N, such that 5(fh n )® X n £ C, 5(fh n+ i)®X n £ £ . 
Furthermore, 

(4.20) {j h m,Jh n } c = \I h m,Jh n J cl = , for all m,n £ Z + . 

Namely the Hamiltonian functionals J h n ,n £ Z+, are integrals of motion for all the vector 
fields X m £ q 9 , m £ Z+, and these vector fields are Hamiltonian with respect to both C and C . 
Finally, all vector fields X m , m £ Z+, commute, provided that Ker^C D Ker%C' = 0. 

Proof. Proposition 4.16 guarantees the existence of an infinite (£, £')-sequence {uj n , X n } ne % + , 
extending the given one. By Theorem 4.13, all elements uj n £ Q 1 are closed, hence, by Theorem 
3.2, they are exact, namely u n = 5(f f n ). Finally, equation (4.20) follows by Proposition 4.16. 
The last claim of the Corollary follows from Remark 4.6. □ 

4.6 The non-linear Schrodinger (NLS) hierarchy. Let V = C[u^ n \v^ \ n £ Z+]. Con- 
sider the following pair of subspaces of © g 9 ~ V® 2 © V 2 : 

-{(4>Uf A)} ' tf -{(0-(7)}- 

where /, g are arbitrary elements of V for while for C they are such that |, £ V. 
Proposition 4.18. (a) (£,£') is a compatible pair of Dirac structures. 

(b) Condition (i) of Theorem 4-13 holds. 

(c) Let ho = 0, h\ = \{u 2 + v 2 ) £ V, let F ^ 1 £ V® 2 be their variational derivatives, and 
let P° = ( ~ U J £ V 2 ~ Q 9 . Then 

F° © P° £C , F 1 © P° £ £ , 

namely ^F^^F 1 ^ and {P } form an (C, C) -sequence, i.e. condition (ii) of Theorem 4-13 
holds (forN = 0). 

(d) The orthogonality conditions (4.19) hold (for N = 0), namely 

(Ci^ffiCF 1 )^ C7T 2 (£') , (CP )^ C 7Tl(£). 

Proof. Clearly, C is a Dirac structure, since it is the graph of the operator H = 

V® 2 — > V 2 , which is a Hamiltonian operator. We next prove that £ is a Dirac structure. In 
particular, we start by showing that C C V® 2 © V 2 is a maximal isotropic subspace. Namely, 
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given B G V® 2 © V 2 , we want to prove that (A, B) = for all A G C if and only if B G C. An 
arbitrary element of C has the form 

and we can always write B G V® 2 © V 2 in the form 

B 

Clearly, B G C if and only if the following conditions hold: 




(4.21) 



6 — a = dr G <9V , c = 3 



4ot , d = d 



+ dr 



Aur . 



On the other hand, the equation (A, B) = gives 



— c H a + ( a 

u 1i 



4up 



+ d 



f + dg 



+ 4ti5 



0, 



which, after integration by parts, becomes 



c d 1 „ fa 
- + a - 

V U V \v 



u \u 



+ g[A{b-a)-d 



+ d -d 



u \ u 



Since the above equation holds for every /, g, we get the two identities: 



(4.22) 
(4.23) 



d 1 



V U V 



- + d -)--d - 



l 



A(b-a)-d - )+d l-d - 



1 



0. 



The second equation, in particular, implies that b — a = dr G dV, which is the first condition 
in (4.21). Equations (4.21) and (4.22) then become 

c 1 



(4.24) 
(4.25) 



— -dl- 

V V \v 

d !* 
4r h -d 

u u 



+ 4r = 0, 
0. 



which are the same as the second and third conditions in (4.21). This proves that £ is a maximal 
isotropic subspace of V® 2 © V 2 . 

In order to show that £ is a Dirac structure, we are left to prove that it is closed under the 
Courant-Dorfman product (4.4). Let A and B be two elements in C\ they can be written in 
the form 



(4.26) A 

(4.27) B 



FffiP 



G®Q, G 



a+db i 



P 



d [2 



Avb 



qI a±db\ +4ub 



-dd 



, Q 



d £ 







c+dd 



- Avd 
+ 4ud, 



dF + A 



dG + A 



u 



d. 
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Let A o B = H R, with H=[ \e V wz and R G V 2 . In order to prove that Ao B e £, 
we need to show that 

-v 



(4.28) uw - vz = dr G dV , R = dH + 4 



u 



for some r G V. By formula (4.4) for the Courant-Dorfman product and using equations (1.15), 
(3.46) and (3.47), we have: 

(4.29) H = D G {d)P + D* P (d)G - D F (8)Q + D* F (d)Q , 

(4.30) R = D Q (d)P - D P (d)Q . 

Using the formula for P in (4.26) we get, using Lemma 1.8, 

( be— ad+bdd— ddb 
ad U -bc 
v 

Hence, the vector H G V® 2 in (4.29) becomes 

( \ / be— ad+bdd— ddb 

w J =D G {d)P-D F {d)Q + A(^ 

With similar computations involving Lemma 1.8 we also get 

/ D c (d)P _ c /Qa+db 4 b \ \ 



and 

/ Dg(d)Q _ a (xc+dd , 4 a \ 

Putting together equations (4.32), (4.33) and (4.34), we get, after a straightforward computa- 
tion, that 

uw-vz = d (D d (d)P - D b (8)Q + ° 



uv 

namely the first condition in (4.28) holds with 

cdb — add 



(4.35) r = D d {d)P - D b {d)Q + 



uv 



Next, let us compute the vector R G V 2 defined in (4.30). Using the last formulas in (4.26) and 
(4.27), and Lemma 1.8, we get, after a straightforward computation, 



(4.36) R = d{D G (d)P-D F (d)Q)+4^ J \ {D d {9)P-D h {d)Q) +4 ( 



dd± - bd% 
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Puttin together (4.32), (4.35) and (4.36), we get 



c+dd jpi a+db \ / bc—ad+bdd—ddb \ / add—cdb 

u , I 1 I „. u 



R - dH ~ 4 [ U ) r = A [ bd dd±-t% U )~ A9 y J-*ycdb=add 

It is immediate to check that the RHS above is zero, thus proving the second condition in (4.28). 

In order to complete the proof of part (a), we are left to show that £ and £' are compatible, 
namely that they satisfy the integrability condition (4.16). By the definition (4.14), the relation 
Ncc C V 2 © V 2 consists of elements of type P®P' , where P is as in (4.26) and P' = JF, with 

J = f ^ ^ j . Moreover, it is not hard to check that the adjoint relation J\f£ ^ C V® 2 © V® 2 

consists of elements of type F ®F' , where again F is as in (4.26) and F' = —JdF + A f U ~\ b. 

Therefore, the integrability condition (4.16) reads as follows: 

(4.37) (H, [P, Q]) - (H', [P, Q']) - (H 1 , [P', Q]) + (H", [P\ Q']) = , 

where 

a+db 
u 
c 

V 

c+dd 



(4.38) P = dF + 4J J b , P' = JF , F = 

(4.39) Q = dG + 4J ( M d , Q' = JG , G = 

(4.40) H = ^ f ^ dg ^ , H' = -JdH + 4 ^ " ^ g , H" = -JdH' + 4 



and where f,g,h£ V are related by the following identity 

(4.41) ± Oh+ l a L + l a l±°£. = . 

UV V V u u 

In (4.37) [ •, • ] denotes the bracket of V 2 ~ Q 9 , namely (cf. (1.15)), [P, Q] = D Q (d)P - D P {d)Q, 
and ( •, • ) denotes the usual pairing V® 2 x V 2 — ► V/dV. Notice that, by construction, we have 

(4.42) F © P , G © Q , H®JH', H' © JiT" G £ . 
Since >C is closed under the Courant Dorfman product, we have that 

(F © P) o (G © Q) = {D G (d)P - D F (d)Q + D* P (d)G + D* F (d)Q) © [P, Q] G £ . 
Hence, since -ff © Jiif' G £ and £ is a maximal isotropic subspace of V® 2 © V 2 , it follows that 

{H,[P,Q]) = -(D G (d)P-D F (d)Q + D* P (d)G + D* F (d)Q,JH') 

(4.43) = (H',JD G (d)P-JD F {d)Q + JD P {d)G + JD F {d)Q). 

In the last identity we used the fact that J : V® 2 — ► V 2 is a skew-adjoint operator. Recalling 
the identities (4.38) and (4.39), and using Lemma 1.8, we also have 

(4.44) (H',[P',Q]) = (H',D Q (d)JF-JD F (d)Q), 

(4.45) (H',[P,Q'}) = (H',JD G (d)P-D P (d)JG). 
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Furthermore, using the fact that J is skew-adjoint, we can easily get 

(4.46) (H", [P\ Q']) = (Dp(d)JG - D G (d)JF, JH") . 

We can then put together equations (4.43), (4.44), (4.45) and (4.46), to rewrite the integrability 
condition (4.37) as follows 

(4.47) (H',Dp{d)JG-D Q (d)JF+JD*p{d)G+JD* F {d)Q)+{D F (d)JG-D G {d)JF,JH") = 0. 

Since H' © JH" G C and C is maximal isotropic, in order to prove equation (4.47) it suffices to 
show that 

(4.48) (D F {d)JG - D G {d)JF) ® (D P (d)JG - D Q {d)JF + JD* P (d)G + JD* F (d)Q) G C , 

whenever F © P, G@Q G C. We start by proving that D F {d)JG - D G (d)JF G tq(£). Indeed, 
it follows by a straightforward computation involving Lemma 1.8, that 

/ Dg{d)JG-D c {d)JF \ 

D F (d)JG - D G (d)JF = I Da ^Q)j G _ Dc ^Q)j F+ Q^ Db ^Q-)j G _ Dd {Q) JF ^ l> 
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namely the RHS has the form ^ z+^m J G ~K\{C), with 

(4.49) z = D a (d)JG- D c {d)JF , w = D b {d)JG - D d {d)JF . 
In order to prove (4.48) we are thus left to prove the following identity: 

D P {d)JG - D Q {d)JF + JD* p {d)G + JD* F (d)Q = d(D F (d)JG - D G {d)JF) 

(4.50) +4J ( U v J (D b (d)JG - D d {d)JF) . 

Using equations (4.38)-(4.39) and Lemma 1.8, it is not hard to check that 

(4.51) D P (d)JG = d(D F (d)JG) + 4 J ^"J D b {d)JG-4Gb, 

(4.52) Dq{8)JF = d(D G (d)JF) + 4 J f M £> d (0) - 4Fd . 

Hence, using (4.51) and (4.52), equation (4.50) becomes 

(4.53) JD* P {d)G + JD* F (d)Q + 4Fd - 4G6 = . 

Equation (4.53) can be now checked directly using equation (4.31). This completes the proof 
of part (a). 

Let us next prove part (b). From the expression of C and C, and from the definition (4.14) 
of Nc,c-> it follows that ^(A/^/y) consists of elements of the form 

/ f+df) \ 

(4.54) ^ z « J GV 2 , 
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where f,g are such that 9 g V. Suppose now that M(d) is a 2 x 2 matrix valued 

differential operator such that 

(4.55) fP-M(d)Q = 0, 

for all P,Q G t^i(Mc £>)■ Recall that P G KiiNcc) ls an arbitrary element of the form (4.54). 
Therefore, if (4.55) holds for every such P, we deduce, by Proposition 1.3, that M(d)Q = 

const, f U j . Since this has to be true for every Q of the form (4.54), we conclude, by a simple 
differential order consideration, that M(d) = 0, thus proving (b). 

For part (c), we just notice that -F°©P°=^||^©^ ^ is the element of £ correspond- 
ing / = and g = \, while F 1 © P° = [ U | © ( I if the element of £ corresponding to 

\ v J V u J 

f = u and g = v. 

We are left to prove part (d). The first orthogonality condition is trivial, since ^(C) = V 2 . 

For the second one we have that, if F = ( Z , | G (CP ) ± , then 

y w/u J 

jF.P° = J(w-z) = 0, 
namely w — z G dV, which exactly means that F G iri(£). □ 

Corollary 4.17 allows us to extend the (£, £')-sequence {F°, F 1 } , {P } to an infinite (£, £')- 
sequence {F n } ne z + , {P n }n& + - Moreover, by Proposition 4.18, all elements F n G V® 1 ~ 
n 1 , n G Z+, are closed, hence, since V is normal and by Theorem 3.2, they are exact: F n = 

Sh n \ 

$ff n I for some h n G V, for every n G Z + . This gives us an infinite hierarchy of Hamiltonian 
~st J 
equations (n G Z + ): 

/ du \ 

(4-56) f = P" , 

\ dt / 

for which all J/i m > m G Z+, are integrals of motion. It is easy to compute the first few integrals 
of motion and evolutionary vector fields of the hierarchy: 

ft2 =i ( w-„'„), f 2 =( ""„ +2 : ( f 2 +,,2 u 

2 2 V ; \ -u' ) ' \ -u" - 2u(u 2 + v 2 ) J 

U _ 1 C/2 ,2 ( 2 , „2%2\ p3 _ / - 2«(U 2 + V 2 ) 



pa 



V" 



2d(u(u 2 + v 2 )) - 4u(W - to') 
- 2d{v(u 2 + v 2 )) + 4u(W - m/) 

b -k V - v'u") + 2(v 3 u' - u 3 v>) F 4 — H" ~ 2d ( v{u2 + ^9 + 4U{VU> ~ UV>) ^ 
hi- 2 (uv vu) + 2{vu u v),* -^ w + 2d(u(u 2 + v 2 )) +4v(vu>-uv>)) 
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Therefore the first three equations of the hierarchy (4.56) are: 

J_( u \ _ ( ~ v \ _L( u \ - ( u ' \ A ( u \ - ( v " + 2v ( u2 + 
dt v v )~ v « / ' * v u / " v ^ J ' rf*2 v w / v _u " ~ 2n (^ 2 + y2 ) 

The third equation is called the NLS system. The first four integrals of motion are h n , n = 
1,2,3,4, written above. It is not hard to see by induction on n that the 1-forms F n , n > 1, 
have differential order n — 1, hence they, and consequently the h n , are linearly independent. 
Finally, Ke^/Z' = 0, hence all evolutionary vector fields Xpn commute, and they are linearly 
independent for n > 1 since the F n, s are. Thus the NLS system (and the whole NLS hierarchy) 
is integrable. 

In fact, using the conditions F n © P n £ C and F n+1 © P n G it is not hard to find a 
recursive formula for F n and P n , n G Z+. We have 

/n. \ / fn + l+dgn + l \ 

A ) » P "= ^ , n>0, 

u / \ u / 

and the elements / n , g n G V are given by the recursive equations: /o = 0, go = 4, and 



( ^ n + ® 9n \ + 4uvg n 
ogn+i = -ud — - vd 



4.7 The Hamiltonian case. In Section 1.3 we gave a definition of a Hamiltonian operator 
in relation to Poisson vertex algebras (see Definition 1.17). As we shall see in Theorem 4.21, 
the following definition is equivalent to it in the case of differential operators. 

Definition 4.19. A Hamiltonian operator is a linear map H : il 1 — > such that its graph 

g(H) = {u® h{oj) | uj e n 1 } c n 1 © g 9 , 

is a Dirac structure. 

In the case of the Dirac structure given by the graph 0(H) of the Hamiltonian operator H, 
the corresponding space of Hamiltonian functionals J-(Q(H)) is the whole space f2°, while the 
space of Hamiltonian vector fields 7i(Q(H)) is the Lie subalgebra H(SQP) C Q 9 . Hence, by the 
discussion in Section 4.3, the space Q, = V/dV acquires a structure of a Lie algebra with the 
bracket 

(4-57) {Ifj9} H = fH(6(ff))(g) = (h (6 (J f)) , 5 (J g)) £ 0°, 

and the evolution equation associated to a Hamiltonian functional Jh £ $7° takes the form 

(4-58) 1/ = H(5(Jh))(f). 

Notice that, in the special case of a Hamiltonian operator H : = V®^ — > g 19 = given by a 
matrix valued differential operator -ff(<9) = (Hij(d))i t j^i, the above formula (4.57) for the Lie 
bracket reduces to equation (1.70), and the evolution equation (4.58) reduces to (1.69). This 
coincidence is not surprising, as will be clear from Theorem 4.21. 
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Proposition 4.20. Let H : Q 1 — > q 9 be a linear map, and let G{H) = |cj © H{u) \ u E SI 1 } C 
© g 9 6e iis graph. 

(a) The space Q(H) C ft 1 (B Q 9 is isotropic with respect to the bilinear form ( , ), i.e. Q(H) C 
Q(H)- 1 -, if and only if it is maximal isotropic, i.e. Q(H) 1 - = G(H), which in turn is 
equivalent to saying that H is skew- adjoint: 

(4.59) (w, H(r])) = -(n, H(u)) , for all u, i/etf. 

(b) A skew-adjoint operator H : fi 1 — > g 9 is a Hamiltonian operator if and only if one of the 
following equivalent conditions holds: 

(i) for every 0^2,^3 E fi 1 , we have 

(4.60) (L h ^u 2 ,H(uj3)) + (L h ^uj 3 ,H(ux)) + (L H ^ui, H(u) 2 )) = 0, 
(mJ for every u),rj E SI 1 we Zioue 

(4.61) H(L H[u)) (r,) - l h[v) (8lo)) = [H(u), H(r,)] . 

Proof. Equation (4.59) is equivalent to (u> © H(uj),t) © H(rj)) = for every u>,n E , which 
in turn is equivalent to G{H) C 0(11)-*-. On the other hand, if if is skew-adjoint, we have 
(u; © X, 7] © H(r])) = (n,X — H(u))), and this is zero for every n E if and only if X = H[uj). 
Hence Q(H) 1 - C G(H), thus proving (a). For part (b), we notice that equation (4.60) is the same 
as the integrability condition (4.8) for the Dirac structure Q{H). Hence, by Remark 4.4, H is a 
Hamiltonian operator if and only if (i) holds. Moreover, by definition of the Courant-Dorfman 
product, we have 

(w9%)) o (v®H(t])) = (l H (u,)(v) - t>H( v )(Suj) © [H(uj),H(7])] , 
and this element lies in G(H) if and only if (4.61) holds. □ 

Theorem 4.21. Let H : Q 1 = V®^ — ► Q d = V e be a matrix valued differential operator, i.e. 
H = (Hij(d)) . Then Definitions 1.17 and 4-19 of a Hamiltonian operator are equivalent. 

Proof. Due to Propositions 1.16 and 4.20, the only thing to check is that equation (4.61) reduces, 
if H is a differential operator, to equation (1.50). Let uj = u>f, V = wg for F,G6 V , be as in 
(3.55). Equation (4.61) then reads 

(4.62) h(l h(u , f) (lo g ) - l h{uig) {8uj f )) = [H{u F ),H(u G )] . 

Using equation (3.47) for Lx p {&f), the element in V® e corresponding to L h ^ f -^(ujg) E SI 1 is 
DG(d)H(d)F-\-D* H ^ F (d)G. Similarly, using equation (3.46), the element in V® 1 corresponding 
to Ih(u g ){o~ojf) E ft 1 is Dp(d)H(d)G — D* F (d)H(d)G. Furthermore, recalling the correspon- 
dence (1.11) and formula (1.15), the element in V e corresponding to [H(uf),H(ug)] E Q d is 
D H (Q^ G (d)H(d)F — Djjig\F(d)H(d)G. Putting together the above observations, we can rewrite 
equation (4.62) as (1.50), and vice versa. □ 
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Next, we discuss how the notion of compatibility of Dirac structures is translated in the 
case of Hamiltonian operators. First we recall the definition of the Schouten bracket [D]. 

Definition 4.22. If H, K : Q l — > are skew-adjoint (with respect to (, )) operators, their 
Schouten bracket is a tri-linear map [H,K]sb '■ J^ 1 x fi 1 x fi 1 — > f2°, given by (cf. equation 
(4.60)) 

(4.63) [H,K] SB (ui,u 2 ,U3) = (ui,H(L K ( U3 )U) 3 )) + (wi, K(L h ^ 2 )uj 3 )) + (cycl. perm.). 
We have the following corollary of Proposition 4.20(b): 

Corollary 4.23. A linear map H : Q 1 — > q 9 is a Hamiltonian operator if and only if it is 
skew-adjoint and [H,H]sb = 0. 

Proposition 4.24. Let H, K : S7 1 — ^ g a be Hamiltonian operators. Their graphs G{H) and 
Q{K) are compatible Dirac structures (see Definition 4-12) if and only if we have 

[H,K]sb(ui,W2,U3) = 0, 

for all u>i, LO2, 0J3 £ , such that K(u>2) £ Im (H), H(u>2) G I m (K). 
Proof. For the Dirac structures G(H) and G(K), we have 

(4.64) Mg(H),g(K) = {H(V) © K( V ) \ r, G ft 1 } C fl° Q d . 
The adjoint relation is, by the non-degeneracy of ( , ), 

^g(H),Q(K) = © I H(u) = K{J)} C ft 1 © n 1 . 
Then by Definition 2.1 G(H) and G{K) are compatible if for all u,u>',uj" £ ft 1 such that 

(4.65) = K(u/), = 
and for all 77, # £ fi , we have 

(4.66) 

- (v>, [K( V ), H(9)}) + (a/', [K(rj), K(0)]) = . 
We next use equation (3.16) to rewrite (4.66) as follows: 

{L H{v) u, H(9)) - {L H{ri) u f , K{9)) - {L K{rj) J, H(9)) + {L K{r]) J\ K(9)) 

(4.67) 

- L H(V) (u, H{9)) + L H{ri) {J, K{9)) + L K[V) (u/, H{9)) - L K[v) (u", K{9)) = . 

By (4.65) and skew-adjointness of H and K, the last four terms in the LHS of (4.67) cancel. 
Moreover, by assumption, [H, H]sb = [K,K]sb = 0, so that equation (4.67) becomes 

( v ,H(L H{uj) 9)) + {u,H(L H{e)V )) + {9,K{L H{J1) J)) 

+ (0,H(L K{v) u/)) + ( V ,K(L K(wl/) 9)) + (u",K(L m n)) =0. 

If we then use again (4.65), the above identity becomes [H, K]sb(Vi w 'i 9) = 0. □ 
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Remark 4.25. Let H,K : O — ► q" be Hamiltonian operators. It is clear from Definition 4.22 
that the Schouten bracket [H, K]sb is linear in each factor. Hence the following conditions are 
equivalent: 

(i) any linear combination aH + (3K is a Hamiltonian operator, 

(ii) the sum H + K is a Hamiltonian operator, 

(iii) [H,K]sb = 0. 

When the above equivalent conditions hold we say that (H, K) form a hi- Hamiltonian pair (cf. 
Definition 2.1). It follows from Proposition 4.24 that if (H, K) form a bi-Hamiltonian pair, then 
the graphs G(H) and Q{K) are compatible Dirac structures. 

Remark 4.26. Theorems 4.13 and 4.21 and Remark 4.25 imply Theorem 2.7. Indeed, due 
to Theorem 4.21 and Remark 4.25, if (H,K) form a bi-Hamiltonian pair (as in Definition 
2.1), the corresponding graphs G(H) and Q{K) are compatible Dirac structures. In order to 
apply Theorem 4.13, we need to check that assumptions (i) and (ii) of this theorem hold. For 
condition (i), by formula (4.64) we have that ^2{J^g<H),Q{K)) = Im(-ftT) C Q 9 . Hence condition 
(i) of Theorem 4.13 becomes the following: if F G V®^ is such that 

(4.68) f(K(d)G 2 ) ■ (D F (d) - D* F (d)) (K{8)G l ) = , for all G 1 , G 2 G V® e , 

then D F (d) — D F (d) = 0. For this we used the identifications VL 1 ~ V® e and q 9 ~ V e , and 
formula (3.45). Using the fact that K is skew-adjoint, equation (4.68) becomes 

JG 2 - (K(d)(D F (d)-D* F (d))K(d))G 1 = 0. 

Hence, recalling that the pairing (• , •) is non-degenerate, we see that condition (i) of Theorem 
4.13 is exactly the non-degeneracy assumption for K in Theorem 2.7. Next, suppose, as in 
Theorem 2.7, that {F n } n=0 ,.., N C V® 1 is an (H, ET)-sequence, and that F° = ^f, F l = 
Then we have uj F n X H{d)pn G Q{H) and, since H(d)F n = K(d)F n+1 , we also have 
ujp n +i © X H ^ F n G G(K), which is condition (ii) of Theorem 4.13. Then, by Theorem 4.13(b), 
all the elements uj F n , n = 0, • • • , N, are closed, thus proving Theorem 2.7. 

Remark 4.27. Let H(d) be a Hamiltonian differential operator, and let Jh G V/dV be a local 
functional, and consider the corresponding Hamiltonian evolution equation (4.12), where X = 
X H (Q\6h- This equation is integrable with respect to the Dirac structure G{H) if and only if it 
is integrable in the sense of Definition 1.27. 

4.8 The symplectic case. In Section 3.5 we gave a definition of a symplectic operator as 
a closed 2-form in the variational complex. As we shall see in Theorem 4.30, the following 
definition is equivalent to it in the case of differential operators. 

Definition 4.28. A symplectic operator is a linear map S : — > such that its graph 

G(s) = {S(x) © x i x g q 9 } c n 1 © q 9 , 

is a Dirac structure. 
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In the case of a Dirac structure given by the graph G(S) of a symplectic operator S, the 
corresponding space of Hamiltonian functionals is !F(G(S)) = 5~ 1 (S(g 9 )), while the space of 
Hamiltonian vector fields is TL{Q{S)) = S~ 1 (5fl°). Hence, by the discussion in Section 4.3, the 
space 5~ 1 (S(g 9 )) C V/dV acquires a structure of a Lie algebra with the bracket 

(4.69) {//, jg} s = JX(g) , where X G g 9 is s.t. S(X) = 5(jf) . 

Moreover, given a Hamiltonian vector field X G 7i{Q{S)) associated to a Hamiltonian functional 
Jh G J-(Q(S)), i.e. such that S(X) = 5(Jh), the corresponding evolution equation is (4.12). 

Proposition 4.29. Let S : g 9 ->• Q, 1 be a linear map, and let Q(S) = {S(X) © X \ X £ q 9 } C 
tt 1 (B Q 9 be its graph. 

(a) The space G(S) C © Q 9 is isotropic with respect to the bilinear form ( , ), i.e. G(S) C 
G(S)- 1 -, if and only if it is maximal isotropic, i.e. G{S) = G{S), which in turn is equivalent 
to saying that S is skew-adjoint: 

(4.70) (S(X), Y) = -(S(Y),Y) , for all X, Y e q 9 . 

(b) A skew-adjoint operator S : g 9 — > is a symplectic operator if and only if one of the 
following equivalent conditions holds: 

(i) for every Xi,X2,X% £ g 9 , we have 

(4.71) (L Xl S(X 2 ),X 3 ) + (L X2 S(X 3 ),X 1 ) + (Lx 3 S(X 1 ),X 2 ) =0, 

(ii) for every X,Y E g 9 we have 

(4.72) L X (S(Y)) - ly(8S(X)) = S([X,Y}) . 

Proof. The proof is the same as that of Proposition 4.20. □ 



Theorem 4.30. Let S : g 9 = V e ^ fl 1 = V®^ be a differential operator, S = (5^(5)). . 
Then Definitions 3.11 and 1^.28 of a symplectic operator are equivalent. 

Proof. In view of Proposition 4.29, we only have to check that, if S = (Sij(d))i : j^j : V — > V® 1 
is a skew-adjoint matrix valued differential operator, then equation (4.72) reduces to equation 
(3.48). Let X = X P ,Y = Xq, with P, Q G V*. We can use equations (1.15), (3.46) and (3.47) 
to rewrite equation (4.72) as follows: 

D S (d) Q (d)P + D* P (d){S(d)Q) - D s{d)P (d)Q + D* s(d)p (d)Q = S(d)(D Q (d)P - D P {d)Q) . 

Written out explicitly, the above equation reads (for k G /): 

E f%^(^) + (-^(^y(%(a Wi )) 



(4.73) -^PH^J + (-^[-^^Qi 



A' 



i,jGl,n£'Z+ \ j 
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We next use Lemma 1.2 to rewrite, in the LHS of (4.73), the first term as 

the third term as 
and the last term as 

E (- 9 ) n (^(f#^) + fly (Wft)) • 

ijeJ,nez + V V <K ^ 9 < / 

The first term in the RHS of (4.73) cancels with the second term in (4.74), and similarly the 
second term in the RHS of (4.73) cancels with the second term in (4.75). Furthermore, by 
the skew-adjointness assumption on 5*, the second term in the LHS of (4.73) cancels with the 
second term in (4.76). Hence, equation (4.73) can be rewritten as 

Since the above identity holds for every P, Q G V , we can replace d acting on P by A and d 
acting on Q by fx. We thus get that (4.72) is equivalent to the following equation (for k G I): 

(4.77) V f ^#A™ - + (-A - /, - «)» = • 

To conclude, we just notice that, by the Definition 1.23 of the Beltrami A-bracket, equation 

(4.77) is the same as equation (3.48). □ 

Next, we study how the notion of compatibility of Dirac structures is translated in the case 
of symplectic operators. Let S : V e —> V® e and T : V e — > V®^ be symplectic operators, and 
consider the corresponding Dirac structures G(S) and G(T). Recalling the definition (4.14), we 
have 

(4.78) Mg{S),Q{T) = {X © X' G Q 9 © Q 9 I S(X) = T(X') } , 
and, similarly, recalling (4.15), we have 

(4.79) Ng {S)t g(T) = {cu © w G ft 1 © ft 1 | (u, X) = {J ', X') , for all X © X' G Afg(3),g(T)} ■ 
We also consider the space 

(4.80) M g{S ),Q{T) = {S(X) © T(X) | X G Q 9 } Cll 1 ©!! 1 . 

It is immediate to check that -Mg(s),g(T) C Ngrs) Q(T) ^ or ever y P a i r °f symplectic operators 5 
and T, while the opposite inclusion is not necessarily true. 
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Definition 4.31. A pair (S,T) of symplectic operators is called strong, if 
( 4 - 81 ) ^g{S),g{T) = M 9{s),g(T) ■ 

Proposition 4.32. (cf. [DJ) Let (S,T) be a strong pair of symplectic operators. The following 
three conditions are equivalent: 

(i) G(S),Q(T) are compatible Dirac structures, 

(ii) for every X,X',Y,Y' £ g 9 satisfying S(X) = T(X'), S(Y) = T(Y'), there exists Z £ g 9 
such that 

(4.82) L X ,S{Y) - L X S(Y') = S(Z) , L X ,T(Y) - L X T(Y>) = T(Z) . 

(Hi) for X 1 ,X[,X 2 ,X^,X 3 ,X' 3 £ g 9 such that S(Xi) = T(X' i ), i = 1,2, 3, we have 

(4.83) (L Xl S(X 2 ),X' 3 ) + {L X ,S(X 2 ),X 3 ) + (cycl. perm.) = 0. 

Proof. Recall Definition 4.12. By (4.78) and assumption (4.81), G(S) and Q{T) are compatible 
Dirac structures if and only if we have 

(4.84) (S(Z'),[X,Y])-(S(Z),{X,Y'])-{S{Z),[X',Y]) + (T(Z),[X',Y'\) = 0, 
for every X, X' , Y, Y' , Z, Z' £ Q 9 such that 

(4.85) S(X) = T(X') , S(Y) = T(Y') , S(Z) = T(Z') . 
Using identity (3.16), we can rewrite equation (4.84) as follows: 

(4.86) -(L X S(Z'),Y) + (L X S(Z),Y') + {L X ,S(Z),Y) - (L X ,T(Z),Y') 
L X {S(Z'),Y) - L X {S{Z),Y') - L X ,{S(Z),Y) + L X ,{T(Z),Y>) = 0. 

By equations (4.85) and the skew-adjointness of the operators S and T, we have (S(Z'),Y) = 
(S(Z),Y') and (S(Z),Y) = (T{Z),Y'), hence the last four terms in the LHS of (4.86) cancel. 
Equation (4.86) then reads 

(4.87) {L X 'S(Z) - L X S(Z>), Y) = {L X ,T(Z) - L X T(Z'), Y') . 

Since the above equation holds for every Y,Y' £ g 9 such that S(Y) = T(Y'), by assumption 
(4.81), it is equivalent to saying that there exists an element W £ g 9 such that L X /S(Z) — 
L X S{Z') = S(W), L X ,T{Z) - L X T{Z') = T(W). If we replace Z by Y and W by Z, these 
equations are the same as (4.82). This proves the equivalence of conditions (i) and (ii). Next, 
we prove that equation (4.87) is equivalent to condition (iii). Since S and T are symplectic 
operators, we have, by the integrability condition (4.71), that 

(L X S(Z'),Y) = -{L Z ,S(Y),X)-{L Y S(X),Z'), 
(L X ,T(Z),Y') = -(L Z T(Y'),X')-(L Y ,T(X'),Z). 

Using the above identities, equation (4.87) becomes 

(L Z ,S(Y),X) + (L Y S{X),Z') + (L X >S(Z),Y) + (L X T{Z'),Y') 
+ (L Z T(Y>),X') + (L Y ,T(X'),Z) = 0, 

which is the same as equation (4.83) for X\ = X, X 2 = Z, X 3 = Y, and similarly for X[, % = 
1,2,3. □ 
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Definition 4.33. A strong pair (S, T) of symplectic operators which satisfies one of the equiv- 
alent conditions (i)-(iii) of Proposition 4.32 is called a hi- symplectic pair. 

In view of Proposition 4.32, we get the following corollaries of Theorem 4.13, Proposition 
4.16 and Corollary 4.17, which translate the Lenard scheme of integrability in the symplectic 
case. 

Theorem 4.34. Let V be an algebra of differential functions in the variables {ui}i<^i, and let 
S(d), T{&) be two symplectic operators given by differential operators: V — > V ■ Suppose 
moreover that {P n }o< n <N C V 1 is an (5, T)- sequence, namely 

S(d)P n+1 = T(d)P n , n = 0, . . . , N - 1 . 

(a) If the orthogonality condition 

(5pan c {P n }o<n<Ar) ± C Im S(8) 
hods, then we can extend the sequence {P n }o< n <N to an infinite (S, T)-sequence {P n } n£ x + ■ 

(b) Suppose in addition that (S,T) is a bi-symplectic pair, and that it satisfies the following 
non- degeneracy condition: if, for some F G V® £ , fQ ■ (Dp(d) — Dp(d))P = for all 
P,Q G V e such that T(d)P, T(d)Q G lm{S(d)), then D F {d) - D* F (d) = 0. Suppose, 
moreover, that 

s(d)p° = 6 -^ , nd)p° = 5 -^, 

ou du 
for some local functionals fho, Jhi G V/dV. Then all elements F n = S(d)P n , n G Z + , 
are closed, namely Dpn(d) = D* Fn (d). In particular, if the algebra V is normal, then they 
are all exact, namely there exist local functionals Jh n G V/ dV such that 

F n = S(d)P n = ^ , n G Z + . 
ou 

These local functionals are in involution with respect to both Lie brackets associated to S 
and T (cf. {4.11)): 

{Jh m ,Jh n } s = {Jh m ,Jh n } T = , for all m,n£Z + . 

(c) Suppose, in addition, that Ker (5) D Ker (T) = 0. Then the evolutionary vector fields 
Xpn, n G Z+, commute, hence the equations of the hierarchy 

du „„ 

— = P n , n G Z+ , 
at n 

are compatible, and the local functionals in involution Jh n are their integrals of mo- 
tion. Thus, this hierarchy is integrable, provided that the set {i ?n } ng z + spans an infinite- 
dimensional space. 
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4.9 The potential KdV hierarchy. Let V = C[u, u' , u", . . . ], and consider the following 
pair of differential operators on V (c 6 C): 

(4.88) S(d) = d , T(d) = u" + 2u'd + cd 3 = 2(u') 1/2 do( u ') 1/2 + cd 3 . 

Proposition 4.35. (a) (S,T) form a bi-symplectic pair (cf. Definition 4-33). 
(b) The pair (S,T) satisfies the following non- degeneracy condition (cf. Theorem 4-34(b)): if 
(4.89) f Q • M{d)P = , 

where M(d) is a differential operator, for all P,Q£V such that T(d)P, T(d)Q G Im S(d), 
then M(d) = 0. 

Proof. We already know, from Example 3.13, that S(d) and T(d) are symplectic operators. 
We want to prove that the pair (S,T) is strong (cf. Definition 4.31). Let F © G £ g(ry 
namely 

(4.90) JFP = JGQ , 
whenever 

(4.91) dP = cd 3 Q + 2d{u'Q) - u"Q . 

We need to prove that F © G G ■M.g(s),g(T)i namely that there exists R G V such that 

(4.92) F = dR , G = cd 3 R + 2u'dR + u"R . 

Equation (4.91) implies that u"Q G dV, namely Q = zjrdf for some / G V. Hence, all pairs 
(P, Q) solving (4.91) are of the form 

(4-93) P = cd 2 (%)+2^df-f , Q = %, 

\ u" J u" u" 

for some / G V such that 

(4.94) df G u"V . 

Equation (4.90), combined with (4.93), gives, after integration by parts, that 

(4 , 5) / /W ^) +8( ^ )+F _ 8( |) } .„. 

Note that the space of elements / G V satisfying condition (4.94) contains the ideal generated 
by (u") 2 . Hence, we can use Proposition 1.3(b) to conclude, from (4.95), that 

w 

It immediately follows from (4.96) that F = dR for some R G V. If we combine this fact with 
equation (4.96), we easily get that F and G have the form (4.92). 
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In order to prove part (a) we are left to prove that the pair (S, T) satisfies condition (ii) in 
Proposition 4.32. Let P , P\, Q , Q\ G V be such that 



(4.97) S(d)P = T(d)P 1 , S(d)Q = T(d)Q 1 . 

According to condition (ii), we need to find R G V such that 

Lx Pl {S(d)Q ) - L Xpq {S(d)Q 1 ) = S(d)R , 
L Xpi {T(d)Q ) - L Xpq {T{d)Q 1 ) = T(d)R . 

Using (3.47) the above equations can be rewritten as follows 

(a q R ) Ds(8)Q (d)Pi ~ D s{d)Ql (d)P + D* Pi (d)S(d)Q - D* P{) (d)S(d)Q 1 = S(d)R, 
1 " ] D T{d)Qo (d)P 1 -D nd)Ql (d)Po + D Pi (d)T(d)Q -D Po (d)T(d^ 

In order to check these equations we use the fact that the pairs (Pq,Pi) and (Qo,Qi) are of 
the form (4.93), and Lemma 1.8. A long but straightforward computation, which is left to the 
reader, shows that equations (4.98) hold for 

R = D Qo (d)P l - D Ql (d)P + PidQi - QidP t . 

Let us next prove part (b). Recall from the above discussion that Q G V is such that 
T(d)Q £ Im S(d) if and only if Q = ^rrdf for some / E V satisfying condition (4.94) (in 
particular, for every element of the ideal of V generated by (u") 2 ). By Proposition 1.3(b), 
condition (4.89) implies that M(d)^Jr = 0, for every / G V satisfying (4.94), and this of course 
implies that M{d) = 0, as we wanted. □ 

Let P° = 1, P 1 = u'. We clearly have S(d)P 1 = T(d)P°, namely {P ,^ 1 } is an (S, re- 
sequence. Moreover, we have (Span c {P°, P 1 })' L C {P )^ = dV = Im S(d), namely the or- 
thogonality condition in Theorem 4.34(a) holds. We also have S(d)P° = , S(d)P 1 = 
for fh(j = and fh± = —^J(u') 2 . Hence, all the assumptions of Theorem 4.34 hold, and we 
can extend {P ,^ 1 } to an infinite sequence {P n } n&+ , such that S(d)P n = T(d)P n ~ l = , 
for some local functionals Jh n G V/dV. This gives us an infinite hierarchy of Hamiltonian 
equations associated to both symplectic operators S(d) and T(d) {n G Z+): 

(4.99) 

for which all local functionals Jh m , m G Z + are integrals of motion. 

It is easy to compute the first few terms of the sequence {P n , fh n } n £z + . In fact the whole 
sequence P n , n G Z + , can be obtained inductively by P° = 1 and the recursive equation 

(4.100) dP n+1 = cd 3 P n + 2u'dP n + u" P n , 

while the local functionals jh n , n G Z+, are obtained by solving the variational problem 

— = dP n 

5u 
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We then get 



P 1 = u> , Jh = -\j{u'f , P 2 = cu>" + 2(«') 2 , Jh 2 = J (|(u") 2 - ^u'f) , 

P 3 = cV 5 ) + 5 cuV'' + ^(i/') 2 + ^u') 3 , fh 3 = j\- \c\u"'f + 5 -cu'(u"f - l(uT) ■ ■ ■ 

The corresponding Hamiltonian equations £1X6 clS follows: 

^=u', ^ = cu'" + ^(u') 2 , ^ = c^ + ^u'u''' + \c{u''f + \{u'f 1 ... 

The second evolution equation above is known as the potential KdV-equation. 

We claim that (4.99) is an integrable hierarchy of evolution equations and hence, in par- 
ticular, the potential KdV equation is integrable (see Definition 4.11). First, the elements 
F n = dP n , n > 1, are linearly independent, since, by the recurrence formula (4.100), the high- 
est degree term (in it, it', ... ) in P n is obtained by putting c = 0. In this case the recurrence 
equation (4.100) can be solved explicitly: 

p n _ (2n - 1)!! . 
and the corresponding integrable hierarchy 

dt n 

may be called the dispersionless pKdV hierarchy. We thus get P n = P^q-I- terms of lower 
degree. The linear independence of the P n 's follows immediately. Finally the evolutionary 
vector fields P n commute since Ker (S(d)) n Ker (T(d)) = 0. 

4.10 The KN hierarchy. Let V be an algebra of differential functions in one variable in 
which u' is invertible, for example V = C[u, (u') , it",...]. Consider the following pair of 
differential operators on V (c € C): 

s(d) = (u'y 1 do( u 'y 1 , T{d) = do( u 'y 1 do( u 'y 1 d. 

Recall from Example 3.14 that S(d) and T{d) are symplectic operators. 

Proposition 4.36. (S, T) form a bi- symplectic pair, and the non- degeneracy condition of 
Proposition 4- 35(b) holds. 

Proof. The proof is similar to the proof of Proposition 4.35 and it is left to the reader. □ 

As usual, in order to find an integrable hierarchy, it is convenient to start from the kernel 
of S(d). Let then P° = it', so that S(d)P° = 0. We have the following identities, which can be 
checked directly, 

/ 1 / 1 a\ _ S /l (i/') 2 \ _ 1 / it'" 3 (u") 2 \ _ tt W u"u"' (u"f 
W lit' J ) ~ 5u \2 {u 1 ) 2 ) ~ w ' \u> 2 (it') 2 J ~ {u'f (it') 3 + (it') 4 " 
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These identities can be rewritten in the form T(d)P° = S(d)P 1 = where P 1 = u'" — 

§£jj£ and J hi = kJ&-- In Particular, {P ,^ 1 } is an (S, T)-sequence. Moreover, we 

have (Span c {P°, P 1 })' L C (P )^ = -KdV = Im S(d), hence the orthogonality condition in 
Theorem 4.34(a) holds. Therefore all the assumptions of Theorem 4.34 hold, and we can 
extend {P ,^ 1 } to an infinite sequence {P n } neZ+ , such that S(d)P n = T(d)P n ~ 1 = ^f, 
for some local functionals jh n £ V/dV. This gives us an infinite hierarchy of Hamiltonian 
equations associated to both symplectic operators S(d) and T(d) for which all local functionals 
fh m , m £ Z + are integrals of motion. The 1-st equation of the hierarchy, associated to P , is 

du_ w 3 Q") 2 
~dT l ~ U ~ 2 u' ' 

which is known as the Krichever-Novikov (KN) equation. 

To prove integrability of the KN equation (and of the whole hierarchy 4r- = P n ), it suffices to 
show that the elements F n = S{d)P n , n > 1, are linearly independent. For this, we notice that, 
by the recurrence relation S(d)P n+1 = T(d)P n , the elements P n have the form P n = -u( 2n+1 ) + 
terms of lower differential order. The linear independence of P n , n > 0, and hence of F n , n > 1, 
follows immediately. Finally the evolutionary vector fields P n commute since Ker (S(d)} = Cu', 
and therefore Ker (S(d)) D Ker (T(d)) = 0. 
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